gcapc.R
In gcapc: GC Aware Peak Caller

## ----para, echo = FALSE, results='hide'---------------------------------------
BiocStyle::markdown()
knitr::opts_chunk$set(dev="png",fig.show="hold",
               fig.width=8,fig.height=4.5,fig.align="center",
               message=FALSE,collapse=TRUE)
set.seed(1)

## ----library------------------------------------------------------------------
library(gcapc)

## ----data---------------------------------------------------------------------
bam <- system.file("extdata", "chipseq.bam", package="gcapc")

## ----rc5end-------------------------------------------------------------------
cov <- read5endCoverage(bam)
cov

## ----bdwidth1-----------------------------------------------------------------
bdw <- bindWidth(cov, range=c(50L,300L), step=10L)
bdw

## ----gcbias-------------------------------------------------------------------
layout(matrix(1:2,1,2))
gcb <- gcEffects(cov, bdw, sampling=c(0.25,1), plot=TRUE, model='poisson')

## ----gcapc, results='hide'----------------------------------------------------
layout(matrix(1:2,1,2))
peaks <- gcapcPeaks(cov, gcb, bdw, plot=TRUE, permute=100L)
peaks <- gcapcPeaks(cov, gcb, bdw, plot=TRUE, permute=50L)

## ----gcapcpeaks---------------------------------------------------------------
peaks

## ----refine, results='hide', fig.width=4.5------------------------------------
newpeaks <- refinePeaks(cov, gcb, bdw, peaks=peaks, permute=50L)
plot(newpeaks$es,newpeaks$newes,xlab='old score',ylab='new score')

## ----refinepeaks--------------------------------------------------------------
newpeaks