Limma vs. lm for methylation data
1
0
Entering edit mode
fluentin44 • 0
@45f45212
Last seen 20 hours ago
United Kingdom

Hi,

Im currently using limma to interrogate methylation data using multiple regression models, typically in the format of methyaltion (beta) values ~ phenotypic variable + covariate 1 +...

However I was wondering, what is / is there a benefit or difference for an OLS regression using the lmFit, eBayes and toptable functions of limma to generate results vs. using multiple testing adjusted lm function for each cg?

We've conducted both limma-based and lm-based (parallelised) regressions using the same data, same regression model and seen ~10% difference in the top 100k hits (when arranged by P.Value).

Thanks, Matt

limma MethylationArrayData MethylationArray • 199 views
ADD COMMENT
2
Entering edit mode
@gordon-smyth
Last seen 40 minutes ago
WEHI, Melbourne, Australia

Yes, there's a big benefit from limma's empirical Bayes approach in terms of improved statistical power and better ROC curves, especially if your sample sizes are small.

There are plenty of published papers that show this. Just follow the references given in the limma documentation or do a Google Scholar search for published papers.

ADD COMMENT
0
Entering edit mode

Ok much appreciated, I’ll have a look around.

ADD REPLY

Login before adding your answer.

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