问题描述
我正在使用混合模型进行重复测量方差分析.我想运行一个事后测试来查看交互 TREAT*TIME 的 p 值，但我只设法使用了以下 ghlt Tukey 测试，它没有给我我正在寻找的交互.
I am doing a repeated measures anova with a mixed model. I would like to run a post hoc test to see the p-values of the interaction TREAT*TIME, but I only managed to use the following ghlt Tukey test which do not give me the interaction I am looking for.
library(multcomp)
library(nlme)
oi<-lme(total ~ TREAT * TIME, data=TURN, random = ~1|NO_UNIT)
anova(oi)
summary(glht(oi, linfct=mcp(TIME="Tukey", TREAT="Tukey")))
我要找的是这样的:
summary(glht(oi, linfct=mcp(TIME="Tukey",TREAT="Tukey",TREAT*TIME="Tukey")))
推荐答案
使用 snk.test(model, term="TREAT*TIME", between="TREAT", inside="TIME") 来自 GAD 包，如果您有一个平衡模型和 summary( lsmeans( oi, pairwise ~ TIME*TREAT), infer=TRUE) 来自 lsmeans 如果您的模型不平衡
Use snk.test(model, term="TREAT*TIME", among="TREAT", within="TIME") from the package GAD if you have a balanced model and summary( lsmeans( oi, pairwise ~ TIME*TREAT), infer=TRUE) from lsmeans if your model is unbalanced
这篇关于双向混合模型方差分析的事后检验的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持！