如何解决emmeans 能否将球度校正应用于重复测量对比?
我正在分析一个 2 路重复测量数据集,建模如下:
model= lmer(result ~ treatment * time + (1|subject),data=df)
...每个受试者都接受每次治疗并每次都接受测试。然而,在分析对比时,似乎没有对球形度进行校正。在这里,我使用 emmeans 来测试每个处理和对照处理之间的差异,在每个“时间”级别...
emm <- emmeans(model,c("treatment","time")
contrast(emm,"trt.vs.ctrl",ref="Control",by="time")
当我查看contrast() 的输出时,通过与GraphPad Prism 对同一数据集的输出进行比较,我确认没有G-G 校正。
是否有一种简单的方法可以对对比度进行球度校正分析?
解决方法
感谢评论者确定此解决方案。 afex() 包专为重复测量因子设计而设计,并允许进行适当的修正。
- aov_ez(在 afex 包中)自动应用非球形校正
- emmeans 应该指定模型是多变量的
- 对比(在 emmeans 包中)应指定对“p”的适当调整
因此解决方案看起来像这样......
library(afex)
library(emmeans)
model= aov_ez("subject","result",df,within=c("treatment","time"),type="III")
emm= emmeans(model,c("treatment",model="multivariate")
contrast(emm,"trt.vs.ctrl",ref="Control",by="time",adjust="dunn")
希望这对有类似问题的其他人有所帮助!
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。