how the Log2 FCs is calculated in case of single factor or condition lets say three groups, in LRT it would run as reduced model and the it compares gene across all the 3 condition it is explained that the compares the full vs reduced model to report the significant genes. So in case of LRT when i run the reduced model which sample or group becomes my reference?
Interpreting likelihood ratio tests using DESeq2
Hello,
I am hoping to receive some guidance on a longitudinal RNA-sequencing project I've conducted. I have sequences at multiple ages across 3 different disease models in addition to control. I am running a likelihood ratio test in DESeq2 with the full and reduced models as specified:
full: ~ genotype + age
reduced: ~ age
I realize that I am testing for the effect of genotype, controlling for age, one gene expression with this test. However, I'm unsure how to interpret the DESeq results and particularly the log fold change result. The padj values are simply telling me whether my full model is a better model to capture deviance of gene expression, without taking into account log2fold change differences. Does this mean I could have significant DEGs for which the log fold change in expression is actually 0?
• 3,367 views
•
link
updated
by
1769mkc
•
written
by
jdhindsa1999 •
•
written
by
jdhindsa1999 •
1 answer
The Log2 FCs are not actually directly related to the LRT test p-values, but rather the result of a comparison between two groups (controllable by the names argument). For your purposes just ignore them for now.
The p-values in your case will show genes where the additional coefficient (genotype) produces an increase in log-likelihood that is larger than would be expected if the coefficient was 0 (functionally having no effect).
• 0 views
•
link
written
by
rpolicastro
Log in to answer this question.
More posts like this
-
DESeq2 LRT test with multifactorial design
written by lessismore
-
Wald's test in DESeq2
written by tommy •
-
How to interpret different reduced models in time series analysis with DESeq2?
written by arshad1292 •
-
scanpy differential expression log fold change
written by Skamboj •
-
Determining up-regulation vs down-regulation DESeq2 likelihood ratio test
written by jdhindsa1999 •
-
DESeq multiple factor question
written by jdhindsa1999 •
-
DESeq2 LRT: Does the order of the factor levels matter?
written by basch •
-
DESeq2 likelihood ratio test log fold change interpretation/calculation
written by Argonaut •
-
High logfold change but padj > 0.1 in DESEQ2
written by thjnant
-
edgeR : How to test differential exon usage on case/control study across all time points when accou…
written by sarin.heikki.v •
