Question: DiffBind dba.peakset + bRetrieve=T gives peak scores? or raw counts?
0
gravatar for jchap14
20 months ago by
jchap140
jchap140 wrote:

Hello,

I've got a set of binding site intervals in a DBA object and use dba.count() to count the reads in those peaks. According to the manual, the default score=DBA_SCORE_TMM_MINUS_FULL. But it also mentions "Note that all raw read counts are maintained for use by dba.analyze".

So, the question is, if I then use dba.peakset(bRetrieve=T, DBA_DATA_FRAME) to get a dataframe, will that dataframe contain the TMM_MINUS_FULL scores or the raw read counts?

Thanks,

James

diffbind • 328 views
ADD COMMENTlink modified 20 months ago by Rory Stark2.8k • written 20 months ago by jchap140
Answer: DiffBind dba.peakset + bRetrieve=T gives peak scores? or raw counts?
0
gravatar for Rory Stark
20 months ago by
Rory Stark2.8k
CRUK, Cambridge, UK
Rory Stark2.8k wrote:

Hi James-

When you call dba.peakset() with bRetrieve=TRUE, you will get back the binding matrix using the score set in dba.count(). So in this case, you will get the TMM normalized scores (with control reads subtracted, and using the full library sizes instead of just the reads in peaks).

However you can quickly and easily change the read scores using dba.count() with peaks=NULL. So if you have already done the counting using the default TMM score, you can get the raw unaltered read counts as follows:

> myDBA <- dba.count(myDBA, peaks=NULL, score=DBA_SCORE_READS)
> rawreads <- dba.peakset(myDBA, bRetrieve=TRUE, DataType=DBA_DATA_FRAME)

Likewise, to get RPKM scores:

> myDBA <- dba.count(myDBA, peaks=NULL, score=DBA_SCORE_RPKM)
> rpkm <- dba.peakset(myDBA, bRetrieve=TRUE, DataType=DBA_DATA_FRAME)

Cheers-

Rory

ADD COMMENTlink written 20 months ago by Rory Stark2.8k
Please log in to add an answer.

Help
Access

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