R emmeans能否对重复测量的对比度进行球度校正？

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

model= lmer(result ~ treatment * time + (1|subject), data=df)

。。。每个受试者每次接受治疗，每次接受测试。然而，在分析对比度时，似乎没有对球度进行校正。在这里，我使用emmeans测试每种治疗和对照治疗在每种“时间”水平上的差异

---

emm感谢评论者确定此解决方案。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","time"),model="multivariate")
contrast(emm,"trt.vs.ctrl",ref="Control",by="time",adjust="dunn")

希望这对其他有类似问题的人有所帮助
 
emmeans
不分析您的数据；它解释了你的模型。您拟合的模型包含非常简单的同质性假设。如果你认为你的模型是正确的，他们的EMM和对比度是正确的。但如果你认为这些都是错误的，你需要一个不同的模型来解释你所关心的任何违规行为。也许是多元模型？我建议在stats.stackexchange.com上询问。在多水平回归的背景下，异方差残差是一个困难的问题，你可能会在另一个网站上获得更好的见解。我同意@Simon的观点，关于建模问题的更好建议可以在CV上找到。如果我在评论中过于轻蔑，我会补充说，你可以看看afex软件包。它的
aov_ez
函数（或类似名称）将适用于单变量和多变量模型，提供关于哪个更好的指导，并支持通过emmeans对您决定使用的任何模型进行事后测试。我非常感谢您的评论，谢谢。我将重新发布并检查afex——听起来它是为这类实验设计的。干杯
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")