Interpreting likelihood ratio tests using DESeq2
1
1
Entering edit mode
7 weeks ago
jdhindsa1999 ▴ 10

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?

RNA-sequencing expression differential DESeq2 R gene • 184 views
ADD COMMENT
1
Entering edit mode
7 weeks ago

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).

ADD COMMENT

Login before adding your answer.

Traffic: 2747 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6