# data was created as follows:
## Not run:
# library('MMDiffBamSubset')
# dataDir <- system.file("extdata", package="MMDiffBamSubset")
# library('DiffBind')
# olddir <- setwd(dataDir)
# DBA <- dba(sampleSheet="Cfp1.csv", minOverlap=3)
# Peaks <- dba.peakset(DBA, bRetrieve = TRUE)
# DBA <- dba.count(DBA, minOverlap=3)
# setwd(olddir)
# peaks <- dba.peakset(DBA, bRetrieve=TRUE)
# C <- Counts(MMD)
# idx <- which(C[,1]>150 & C[,3]>150&width(Peaks)>1000&width(Peaks)<5000)
# Peaks <- Peaks[idx[1:500]]
# ## End(Not run)
Run the code above in your browser using DataLab