Question: Diffbind: getting the raw reads
gravatar for veronique.storme
14 months 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 14 months 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 14 months 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 14 months ago by Rory Stark2.5k
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: 140 users visited in the last hour