Question: Diffbind: getting the raw reads
gravatar for veronique.storme
6 days 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,



ADD COMMENTlink written 6 days 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 5 days 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 4 days ago by Rory Stark1.7k
Please log in to add an answer.


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