bamFiles <- system.file("exdata", c("rep1.bam", "rep2.bam"), package="csaw")
data <- windowCounts(bamFiles, width=100, filter=1)
data$param
data <- windowCounts(bamFiles, width=100, filter=1, param=list(readParam(), readParam(minq=100)))
data$param
data[,1]$param
data[,2]$param
plist <- data$param
reformList(plist, restrict="chrA")
reformList(plist, minq=10)
checkList(plist)
checkList(reformList(plist, minq=10))
Run the code above in your browser using DataLab