MAnorm output as DiffBind input.
4
0
Entering edit mode
ajeet • 0
@ajeet-14838
Last seen 2.8 years ago

Hello Everyone/Rory,

I'm working on ChIP and ATAC data sets analyses, however, I need help at this point. My goal is to use pre-normalised data (e.g. using MAnorm) as DiffBind input so that later I can create contrast list (dba.contrast), reports (dba.report) and plot (dba.plotXYZ).

Question: Is there a way by which I could use the output (MAnorm_result_commonPeak_merged.xls) obtained after MAnorm normalization as an input for DiffBind?

Any help or recommendation is appreciated. 

Thanks,

Ajeet 

FYI: I'm aware of these posts e.g.

 question about DiffBind 

 https://rdrr.io/bioc/DiffBind/man/dba.peakset.html

Raw Count Matrix from DiffBind

MAnorm diffbind normalization chipseq • 956 views
ADD COMMENT
1
Entering edit mode
Rory Stark ★ 4.1k
@rory-stark-5741
Last seen 2 days ago
CRUK, Cambridge, UK

I'm looking now adding features to DiffBind to give much more fine-grained control over the normalization, including adding pre-normalized data.

In the meantime, there is a feature that can be used for this. See the man page for dba.peakset(), specifically the Counts parameter. You can input count values explicitly using this parameter, or by writing them to files and using a Counts column in your samplesheet; see man page for dba(). Note that you should not call dba.count() if you are using the Counts parameter.

If you do not want these counts to be further normalized, you should execute the following line of code on your DBA object (myDBA):

myDBA$class[8,] <- 1e+06

This will set all the library sizes to be the same and prevent additional normalization, provided you use the DBA_DESEQ2 method with bFullLibrarySize=TRUE in your call to dba.analyze().

 

ADD COMMENT
0
Entering edit mode
ajeet • 0
@ajeet-14838
Last seen 2.8 years ago

Thanks! Rory,

Sounds promising. Let me try this out and will get back to you.

Ajeet 

ADD COMMENT
0
Entering edit mode
ajeet • 0
@ajeet-14838
Last seen 2.8 years ago

Thanks! Rory,

It worked nicely for me.

Best,

Ajeet 

ADD COMMENT
0
Entering edit mode
Rory Stark ★ 4.1k
@rory-stark-5741
Last seen 2 days ago
CRUK, Cambridge, UK

The latest version of DiffBind has a new interface function, dba.normalize(), that allows vectors of normalization factors and library sizes (or a matrix of normalization offsets) to be specified instead of using the internal normalization calculations. Examples in the help page and vignette.

ADD COMMENT

Login before adding your answer.

Traffic: 317 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6