w*****r 发帖数: 165 | 1 求助这里的大牛:
我的数据有3个group,5个时间点。用 R 里面的 lme 来做repeated measure ANOVA:
imgdata.lme <- lme(measdat ~ Group * Time, random = ~1 | Subject,
data=imgdata, method="REML")
我现在想比较每个时间点的group间(A,B,C) 的pair-wise的区别 (A vs. B, B vs. C
, A vs. C)。不需要跨时间点比较。所以会是:
T1: A-B
T1: A-C
T1: B-C
T2: A-B
T2: A-C
T2: B-C
T3: A-B
T3: A-C
T3: B-C
T4: A-B
T4: A-C
T4: B-C
T5: A-B
T5: A-C
T5: B-C
但需要correct for multiple comparisons,比如Tukey方法。我试了用glht,但对于
该怎么设contrast matrix感觉很晕。求指点。
谢谢大家。 |
|