Entering edit mode
shen_
•
0
@shen_-22158
Last seen 5.1 years ago
I tried to count reads around peak set based on paired-end ChIP-seq data. But I don't want to get a consensus peak set. So I add "peaks=NULL" according to the document of DiffBind.
Here is my code:
> data$config$singleEnd <- F
> data <- dba.count(data,peaks=NULL,bUseSummarizeOverlaps=T)
Error in dba.count(data, peaks = NULL, bUseSummarizeOverlaps = T) :
DBA object must contains only counts