Assessing DMRs between treatment groups: subtracting control group DMRs to account for "baseline" methylation
1
0
Entering edit mode
Ellen O • 0
@ellen-o-12679
Last seen 4.1 years ago
Stellenbosch University, South Africa

I am currently awaiting Illumina methylationEPIC (850k) data and in the mean time figuring out how to assess DMRs with demo data. I am new to these analyses but I am making use of the Maksimovic et al. cross-package workflow (https://f1000research.com/articles/5-1281/v2), which I have found very useful thus far. 

I have three sample groups: healthy controls, poor treatment responders, and good treatment responders. I was wondering if it is possible to "subtract" the differential methylation seen within the controls before comparing the two treatment groups, so that I can account for a baseline/ "normal" level of variable methylation. Would this be possible in limma, minfi or any other packages? Or would it make more sense to just compare all three groups to each other?  

Thanks in advance for any assistance! 

illumina methylationepic 450k limma minfi dmr analysis dmr • 1.2k views
ADD COMMENT
0
Entering edit mode
Aaron Lun ★ 28k
@alun
Last seen 9 hours ago
The city by the bay

Assuming you're fitting a linear model to the M-values, you're saying that you want to correct for the control:

corrected_poor = poor - control
corrected_good = good - control

The null hypothesis would then be:

H0: corrected_poor = corrected_good

However, some arithmetic reveals that this is the same as:

H0: poor - control = good - control
              poor = good

So you might as well just compare the two treatment groups directly, if you're interested in their differences.

ADD COMMENT
0
Entering edit mode

Thanks Aaron!

ADD REPLY

Login before adding your answer.

Traffic: 761 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