Question: limma: Using log2 fold changes as input for a linear model fit
gravatar for Moritz E. Beber
3.5 years ago by
European Union
Moritz E. Beber10 wrote:

My data are log2 fold changes (FCs) between control and treatment usually coming from triple replicated affymetrix microarrays. I also have p-values for those data determined by t-test. As I currently have only those data available and not the replicates themselves, I wonder if it is sensible to use the FCs as input for the limma function lmFit?

The results I get seem reasonable but I have found no way to include the p-value information which makes me uneasy about the reliability. Is this possible in some way or is my only solution to obtain the raw replicate data?

Thank you for your insights.

ADD COMMENTlink modified 3.5 years ago by Aaron Lun21k • written 3.5 years ago by Moritz E. Beber10
gravatar for Aaron Lun
3.5 years ago by
Aaron Lun21k
Cambridge, United Kingdom
Aaron Lun21k wrote:

Your best option would be to obtain the raw data (i.e., the raw intensities from each of the replicated array) and use limma on that. The log-fold change between control and treatment for each gene/probe does not provide any information about the variability of expression across replicates. In fact, it's not obvious to me how you would do the limma analysis at all, if you were to give it only a single log-fold change for each gene.

ADD COMMENTlink modified 3.5 years ago • written 3.5 years ago by Aaron Lun21k

The data are from treatment experiments that have a dose and time dependency, so that I just fit a linear model to the FC as a function of time and dose. However, the lack of variability and significance information of the FCs is a big concern to me, too, so I'll try to obtain the raw data, I suppose.

ADD REPLYlink modified 3.5 years ago • written 3.5 years ago by Moritz E. Beber10
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: 369 users visited in the last hour