Question: Extract RPKM values from Diffbind
0
gravatar for francesca3
4 weeks ago by
francesca320
francesca320 wrote:

Hi everyone, does anyone know if there is a way to extrapolate rpkm values from Diffbind?? Thanks Francesca

chip-seq diffbind rpkm • 104 views
ADD COMMENTlink modified 4 weeks ago by Rory Stark730 • written 4 weeks ago by francesca320
1
gravatar for Rory Stark
4 weeks ago by
Rory Stark730
University of Cambridge, Cancer Research UK - Cambridge Institute
Rory Stark730 wrote:

You can set the read score to DBA_SCORE_RPKM in dba.count(). You can do this when you originally count, or anytime after by calling

DBA <-  dba.count(DBA,peaks=NULL,score=DBA_SCORE_RPKM)

Then you can retrievethe RPKM values in a GRanges object (or data.frame):

rpkm <- dba.peakset(DBA, bRetrieve=TRUE)
ADD COMMENTlink written 4 weeks ago by Rory Stark730

Thank you Rory. I have another question. When I retrieve these RPKM values I obtain very low numbers (ranging from around 0.5 to 2). Which values should I expect for trusted regions? If I retrieve TMM values (DBA_SCORE_TMM_MINUS_FULL) these values range from around 10 to 200. Usually I was counting using DBA_SCORE_TMM_MINUS_FULL. Thank you.

ADD REPLYlink written 4 weeks ago by francesca320
1
gravatar for Rory Stark
4 weeks ago by
Rory Stark730
University of Cambridge, Cancer Research UK - Cambridge Institute
Rory Stark730 wrote:

RPKM values are much lower in general than read counts. The first step is to divide by how many million sequence reads you have in your library. So if you sequenced to a depth of 20M reads, all the counts would be divided by at least 20.

Here's how it's implemented in DiffBind:

rpkm    <- (counts/(width(intervals)/1000))/(libsize/1e+06)
ADD COMMENTlink written 4 weeks ago by Rory Stark730
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 2.3.0
Traffic: 1661 users visited in the last hour