Is it correct to plot Log2FC vs adj p-value from Likelihood ratio test?
1
0
Entering edit mode
Laia ▴ 10
@239caad3
Last seen 6 weeks ago
Belgium

Hi,

I am trying to analyze what genes are affected by the interaction factor in my model.

When I do 2 by 2 contrasts using Wald, the volcano plots look ok (V shape). When I want to plot the genes that were significantly affected by the interaction -reducing it in the LRT- (adjusted p-value this time), then I get noise in the central area of the volcano plot, getting kind of an inverted V shape. I am using -log10(adjp) in LRT, and -log10(pval) in Wald.

Why is that happening? How is that genes with no much fold change still get high adjusted p-values in LRT? If I put it in words, genes that have no significant fold change when compared with the reference condition (log2FC close to 0), still get a significant impact from the interaction factor.. (or whatever factor I test with LRT).

Am I plotting this wrongly?

Thanks,

Laia

DESeq2 adjp Volcano LRT • 185 views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 32 minutes ago
United States

See the information in vignette and ?results about LRT. The coefficient/LFC is just one of many that are represented in the LRT.

ADD COMMENT
0
Entering edit mode

Hey! How did you manage to get the beta values then? word wipe

ADD REPLY

Login before adding your answer.

Traffic: 407 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