Question: Integrating the computed spike-in coefficients for the normalisation before differential peaks analysis
gravatar for Konstantin Okonechnikov
4 weeks ago by

Hi! Have a question about using DiffBind for ChIP-seq data with drosophila spike-ins.

I already have computed peaks from usage of control data and spike-in normalisation coefficients for a cohort of target samples across two conditions. The spike-in coefficients computation is based on drosophila alignment results as described in ActiveMotif documentation for down-sampling (adjustment based on minimum).

Is there an easy way to use these computed coefficients for the DiffBind analysis? I found this post, but there is a link to pipeline which starts from reads, while I would like only to correct the computed scores matrix before differential analysis. The matrix is computed from original samples without subsampling, but with the usage of control.

chipseq diffbind spike-ins • 68 views
ADD COMMENTlink modified 24 days ago by Rory Stark3.0k • written 4 weeks ago by Konstantin Okonechnikov20
Answer: Integrating the computed spike-in coefficients for the normalisation before diff
gravatar for Rory Stark
24 days ago by
Rory Stark3.0k
CRUK, Cambridge, UK
Rory Stark3.0k wrote:

You could retrieve the computed read counts using dba.peakset(), then correct them and use them as in the reference.

ADD COMMENTlink written 24 days ago by Rory Stark3.0k

Thanks a lot for the reply! My further question: how to create the object with adjusted coefficients for further processing? Is there a specific way to set the reference? Here's my code example:

# standard object creation
dtCounts <- dba.count(dt)
# retrieve peaks
peaksRes <- dba.peakset(dtCounts,1:numSamples,bRetrieve=TRUE)
# use spike-in coefficients to adjust results, here's example for one sample
peaksRes$sample1 = peaksRes$sample1 * k1
# how to re-write the target analysis object with adjusted values?
dtCountsAdj <- <???> 
# standard analysis continues....
dtCountsAdj <- dba.contrast(dtCountsAdj, categories=DBA_CONDITION)
dtCountsAdj <- dba.analyze(dtCountsAdj)

Would be grateful for comments

ADD REPLYlink modified 17 days ago • written 17 days ago by Konstantin Okonechnikov20
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: 145 users visited in the last hour