Hi guys,
I have a doubt regarding the logFC calculation in limma. I searched for the answer in Biostars but I have not found anything that help me to understand my doubt. Supposedly logFC is: LOG2(mean (Group A) / mean (GroupB)) Meaning that if a gene logFC has an expression of 50 in Group A and 100 in Group B, the log2 FoldChange would be -1. Thus, the gene would be two times more expressed in Group B. However, I have checked the logFC that appears in the topTable obtained with the functions in limma and the logFC seems to be something else. Here, the logFC represents the difference between the mean of the groups. The value is then elevated to 2 and log transform, to address correctly the difference between the groups. Thus, the value indicates directly the difference between the two groups, it is not the ratio between the groups. So, it is the value of the topTable a proper logFC? Or should I calculate the logFC manually? Maybe I am confused with the terminology and the value of limma is called something else or the "supposedly logFC" that I addressed is named differently.
Thank you for your precious time.
0 answers
No answers yet.
Log in to answer this question.
-
Create the correct contrast in limma when want to combine 2 or more groups
written by Lila M
-
Same results in Limma Toptable as in DESEQ2 results
written by Nona •
-
Can we get the results for Student t test and Welch t test by using limma?.
written by Minushi •
-
Why do I get a big log fold change but small mean change in b value when plotting differential meth…
written by Christine •
-
logFoldChange computed by EdgeR different from naive log2FCcomputed based on reads
written by ZheFrench
-
Calculation of Fold change for already log transformed gene expression data
written by Bioinformatic •
-
sign of logfc in limma
written by im •
-
[limma] LogFC interpretation
written by flappix •
-
Different logFC (log2foldchange) values for genes from limma-voom and other tools (edgeR and DESeq…
written by SMILE
-
Differential Gene Expression Analysis using data_RNA_Seq_v2_expression_median RSEM.Normalized
written by gokce.ouz
Can you give references or links where you found
With only words it is difficult to understand what you mean.
I just realized my mistake. I was dividing values that were log transformed, thus, I should have done a subtraction instead. By doing so, I have the same results that limma does. Meaning that here the logFC is directly the difference between the groups rather the ratio (since the values are log transform). Thank you for your quick answer.
Hello John S.!
Solved by OP.
For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.
If you disagree please tell us why in a reply below, we'll be happy to talk about it.
Cheers!