Differential methilation analysis batch effect removal
Entering edit mode
Mageiros • 0
Last seen 15 days ago

Hi all,

I am trying to analyze 43 RRBS samples using edgeR and I am having some issues as I am very new to this. I mainly follow section 4.7 of the manual (https://www.bioconductor.org/packages/release/bioc/vignettes/edgeR/inst/doc/edgeRUsersGuide.pdf) to perform my analysis.

My issue is that my samples have been sequenced in 4 batches so I would like to remove batch effects.
Applying ComBat to matrix M (see section 4.7.4 of the manual) indeed seems to work n the MDS plots.

How can I incorporate this change though, to my DGEList object before I proceed to differential methylation analysis is unclear to me. Can I just apply ComBat to my DGEList counts matrix?

Should I normalize my data any other way apart from what is described in section 4.7.3 of the manual?

Any help will be truly appreciated.

Kind regards


BatchEffect combat sva edgeR RRBSdata • 115 views
Entering edit mode
Last seen 10 minutes ago
WEHI, Melbourne, Australia

ComBat can't be run on a DGEList, but it is not necessary to do so.

The best method to correct for batch effects is to include batch as a factor in the edgeR linear model. Then edgeR will correct for any batch effects as part of the differential methylation analysis.

In my opinion, it is always better to adjust for batch effects as part of the linear model (see Section 3.4.3 of the edgeR User's Guide) rather than as a separate step. Running batch correction as a separate step tends to change the data in risky ways and can very often result in spurious significant results.


Login before adding your answer.

Traffic: 172 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6