Question: Diffbind: getting the raw reads
gravatar for veronique.storme
2.7 years ago by
veronique.storme0 wrote:

Dear Rory,

I found the solution to my previous questions. However, now I am comparing the results from the available edgeR analysis with the robust edgeR pipeline with QL method based on raw read counts. This is the pipeline that I want to use. I found the raw read counts by

reads.5 = dba.count(prol.5, peaks=NULL, score=DBA_SCORE_READS)
bindingMatrix.5 = dba.peakset(reads.5, bRetrieve=TRUE, DataType=DBA_DATA_FRAME)

However, these counts are different from:

de.5.edger = dba.analyze(contrast.5, method=DBA_EDGER, bFullLibrarySize=FALSE, bTagwise=TRUE, bReduceObjects=FALSE)

and as a result, I also get totally different DE genes.

How do I get the same counts with the dba.count function?

kind regards,



edger diffbind raw reads • 551 views
ADD COMMENTlink written 2.7 years ago by veronique.storme0

Dear all,

I found the answer myself. When I ask for 

dba.count(prol.5, peaks=NULL, score=DBA_SCORE_READS_MINUS)

I can see that when all negative values are put to 1, I obtain the same reads as in 



ADD REPLYlink written 2.7 years ago by veronique.storme0

You can turn off the subtraction of control reads by setting bSubControl=FALSE in the call to dba.analyze().


ADD REPLYlink written 2.6 years ago by Rory Stark3.0k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 315 users visited in the last hour