ANOVA-like test via treat() in limma
1
0
Entering edit mode
@kvittingseerup-7956
Last seen 8 months ago
European Union

Is it possible to do an anova style test (testing multiple contrast jointly) when using treat() instead of eBays() ?

Would be very useful to test against the log2FC cutoff across all conditions.

limma anova • 1.1k views
ADD COMMENT
3
Entering edit mode
@gordon-smyth
Last seen 3 hours ago
WEHI, Melbourne, Australia

No, you can't use treat() to do an F-test. I'm not even sure what a fold-change threshold would mean for an F-test. Fold changes are only defined for individual contrasts, not for an F-test, which in effect tests all possible contrasts.

Of course, if you just want a gene to satisfy a treat() test for several different contrasts, then you could just run treat() for each of those contrasts. The decideTests() function provides an easy way to select genes satisfying several contrasts at once.

ADD COMMENT

Login before adding your answer.

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