emmeans可以对重复测量的对比度应用球度校正吗？

Question

我正在分析一个双向重复测量数据集，建模如下：

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棱镜的输出进行比较，我确认没有G-G校正。

有没有一种简单的方法来实现对比度的球形性校正分析？

Answer 1

感谢评论员确定了这个解决方案。afex()包是专门为重复测量析因设计而设计的，并允许进行适当的校正。

emmeans软件包(在afex软件包中)应自动应用模型校正应指定模型为multivariate

• contrast (在
• 软件包中)应指定对"p"

的适当调整

因此，解决方案应该是这样的……

library(afex)
library(emmeans)

model= aov_ez("subject","result",df,within=c("treatment","time"),type="III")
emm= emmeans(model,c("treatment","time"),model="multivariate")
contrast(emm,"trt.vs.ctrl",ref="Control",by="time",adjust="dunn")

希望这对其他有类似问题的人有帮助！