Differene between the Wald and LRT test for time-series gene expression data
1
0
Entering edit mode
Denise • 0
@860ad694
Last seen 2 days ago
Spain

Hi,

I am using DESeq2 to analyze a (short) time series gene expression dataset. Does anybody know what is the advantage, if any, of using an LRT (likelihood ratio test) for differential expression analysis across all time points versus a Wald test separately for each time point?

I have tried both but I do not see much qualitative differences.

Thank you in advance.

Denise

DESeq2 • 233 views
ADD COMMENT
1
Entering edit mode
@james-w-macdonald-5106
Last seen 4 days ago
United States

The advantage to using the LRT is that it's an omnibus test that tells you if there are changes at any time point. If you have lots of time points you could have any number of Wald tests you would have to do. The alternative is to do the LRT, get the significant genes, and then use e.g., k-means clustering to identify groups of genes with different patterns of expression. You can then make high level statements about those groups, and maybe do gene set testing of some sort to identify likely pathways that are being affected in a time-dependent manner.

Think of it as inspecting the forest versus looking at individual trees. In general the higher level view provides (IMO) a more identifiable signal.

0
Entering edit mode

This makes a lot of sense, thank you very much!

ADD REPLY

Login before adding your answer.

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