Question: limma for metabolite data
gravatar for mahes.muniandy
5 months ago by
mahes.muniandy0 wrote:


My name is Mahes Muniandy and I am a doctoral student working on twin data. I have used limma previously for gene expression as well as methylation data and would now like to use it on my metabolomics data - to detect within twin-pair differences. I have two quesions:

1) Since metabolite data has high multicolinearity, can I safely use my metabolite data as is or should I reduce the data to PCAs and then use the PCA values to represent my metabolites in my limma model?

Here is my design matrix:
design <- model.matrix(~Pair+Smoking+subcutaneousfat+bmi)
fit <- lmFit(metabolites, design)

2) In my paired samples design, I am using a continuous variable (subcutaneousfat) which might or might not be discordant within the twin-pairs - I think this is not a problem but should my data be sorted in a way that the higher subcutaneousfat is always compared against the lower one - or does the model care?

Many Thanks,

Mahes Muniandy,

Finnish Institute of Molecular Medicine,

Helsinki, Finland

ADD COMMENTlink modified 5 months ago by Aaron Lun16k • written 5 months ago by mahes.muniandy0
gravatar for Aaron Lun
5 months ago by
Aaron Lun16k
Cambridge, United Kingdom
Aaron Lun16k wrote:

Are you referring to collinearity between metabolites? This does not matter to lmFit, as a linear model is fitted to each metabolite separately. Correlations between metabolites do not affect the individual model fits. (It might reduce the benefit of information sharing between metabolites via empirical Bayes, but that's another matter.) Collinearity between covariates/factors is more concerning, e.g., if all smokers had higher BMIs; this would reduce power to detect DE caused by smoking or BMI, as the effects of one cannot be separated from the other.

As for your second question - how you order the samples doesn't matter to lmFit, as long as the order is the same in the design matrix. Instead, the interpretation of the coefficients depends on the order of levels in each factor (for categorical variables). See ?factor for more details.

ADD COMMENTlink modified 5 months ago • written 5 months ago by Aaron Lun16k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 104 users visited in the last hour