Question: Does limma calculate pvalues using log2 fold changes or with linear values?
1
gravatar for serpalma.v
3.2 years ago by
serpalma.v40
Germany
serpalma.v40 wrote:

Dear community,

I am currently trying to figure out this question. I have tried different keywords but I do not seem to find the right thread. 

What I want to know is:

- Does limma arrives to a pvalue by using the log2 scale of the normalized array data?

- Does limma arrives to a pvalue using the normalized array data directly and then shows the differences as log2 fold changes in the top table?

If you know the proper thered to refer to, that will be of great help, else, I am grateful for your responses.

 

Thanks!

limma log2 pvalue log2fc • 2.3k views
ADD COMMENTlink modified 3.2 years ago by Aaron Lun24k • written 3.2 years ago by serpalma.v40
Answer: Does limma calculate pvalues using log2 fold changes or with linear values?
3
gravatar for Aaron Lun
3.2 years ago by
Aaron Lun24k
Cambridge, United Kingdom
Aaron Lun24k wrote:

limma - or specifically, lmFit and its downstream functions - will operate on whatever matrix of values you give it. If you're working with microarray data, then - as svlachavas has mentioned - you should be giving lmFit the log-transformed intensities, for various reasons. lmFit will not log-transform the values for you, as it assumes that's already been done. Thus, it's worth checking what was done upstream to ensure that limma is running on the correct values. For example, if you're using normalizeBetweenArrays before running lmFit, then the former function will log-transform the data if your input is an EListRaw object, but not if your input is a matrix. (Specifically, the function returns normalized log-transformed intensities, so there's no distinction between normalization and log-transformation in all of the downstream functions.) In summary, assuming you've done the analysis correctly, then the p-values from limma will be computed from the log-intensities.

ADD COMMENTlink modified 3.2 years ago • written 3.2 years ago by Aaron Lun24k

Thank you very much Aaron,

I normalized the array data with the RMA algorithm.

According to this thread, RMA log transforms the data:

log transform in RMA normalization

ADD REPLYlink written 3.2 years ago by serpalma.v40

Yes, that's correct, the RMA algorithm returns values on a log scale.

ADD REPLYlink written 3.2 years ago by Ryan C. Thompson7.3k
Answer: Does limma calculate pvalues using log2 fold changes or with linear values?
0
gravatar for svlachavas
3.2 years ago by
svlachavas660
Greece/Athens/National Hellenic Research Foundation
svlachavas660 wrote:

Dear Serpalma,

i don't know where you heard or read anything about "linear values" ? Or also to which technology you refer to ? For instance about microarray gene expression data, topTable() from limma provides log2 fold changes ("effect sizes") about your comparisons of interest, along with the relative p-values . But for more specific details and not generalizations, you should definately read the 2 papers below:

1)  http://www.statsci.org/smyth/pubs/ebayes.pdf

2) http://www.ncbi.nlm.nih.gov/pubmed/25605792

Hope that solves any general questions

Best,

Efstathios

ADD COMMENTlink modified 3.2 years ago • written 3.2 years ago by svlachavas660

Dear Efstathios

using the "linear" term was probably the wrong wording, I meant to say "not log scaled". 

Best

ADD REPLYlink written 3.2 years ago by serpalma.v40

Actually, if you don't specify your exact design i can't make any general assumptions-on the other hand, for instance about microarrays, after normalization in the vast majority of existed methodologies,the expression data are already log2-transformed, which is no related with limma. If you are also interested in RNA-seq, you can read the specific part below in the second paper:

"Traditionally,RNA-seq data require specialized software based on the negative binomial or similar distributions .limma however is able to analyse RNA-seq read counts with high precision by converting counts to the log-scale and estimating the mean-variance relationship empirically.The mean-variance trend is converted by the voom function into precision weights, which are incorporated into the analysis of log-transformed RNA-seq counts using the same linear modelling commands as for microarrays".

ADD REPLYlink modified 3.2 years ago • written 3.2 years ago by svlachavas660
Please log in to add an answer.

Help
Access

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