Same results in Limma Toptable as in DESEQ2 results
Dear all,
I am trying to generate exactly the same results in Limma for my result table as I did in DESEQ2 (Differential Expression Analysis).
With the toptable function I am not getting lfcSE and basemean.
Can anyone help me how to extract this information via Limma?
I do get logFC, AveExpr, t, P-Value, adj.Pvalue and B in Limma.
Thank you, Bine
• 3,028 views
•
link
updated
by
jared.andrews07
•
written
by
Nona •
•
written
by
Nona •
1 answer
It would be superfluous to include a SE column in the limma output because it is not needed and because it is immediately computable from the other columns.
• 0 views
•
link
written
by
Gordon Smyth
Yes.
• 0 views
•
link
written
by
Gordon Smyth
Log in to answer this question.
More posts like this
-
Limma for DEG analysis
written by rahma.khoualdia •
-
microarray differential expression analyses
written by bioinformatics •
-
Limma Differential Expression contrast matrix with several variables
written by Nona •
-
understanding log2fold change of 0.31 and -0.25
written by arshad1292 •
-
How to work out Z score for heatmaps for RNA seq dataset
written by v.johnson •
-
DESeq2 normalised data
written by bioinfo_ga •
-
Explanation of Limma logFC calculation
written by John S. •
-
RNA-seq analysis: How to output normalized counts or normalized FPKM values for all genes in each c…
written by sclittlerain •
-
Meta output vs Toptable output
written by vn383 •
-
R: How to convert log2FC (Fold Change) obtained by limma's topTable() function to FC
written by komal.rathi
These are different tools with slightly different approaches. Why would you expect the exact same results?
AveExpris the limma equivalent of DESeq2'sbasemeancolumn, but it has nolfcSEequivalent.Thank you very much. I am replicating the analysis with Limma and for my result tables I would like it to be exactly the same output.
I get that. My question is why would you expect that to be plausible? These methods account for similar things and tend to return largely overlapping results, but they have differences in approach and thus results as well. You're not going to get the exact same results, and I don't know why you'd expect to. Robust findings should be reasonably replicable between the major DE methods (DESeq2, edgeR, limma).