对 class 'mixed' 的模型与 emmeans 进行 post-hoc 比较不起作用

Question

我运行以下混合模型：

fit_mixed_POST <- mixed(correct~probability_simple*letter_position+
(1|PP),data=all_data_cleaned_POST_only, method = c("LRT"),
family=(binomial(link = "logit")))

现在，我想使用 emmeans 对 probability_simple 和 letter_position 执行所有成对比较。但是，当尝试运行（在引入库（emmeans）之后）：

emmeans(fit_mixed_POST, "probability_simple")

我收到以下错误：

Error in ref_grid(object, ...) : Can't handle an object of class  “mixed”

但是，根据文档，emmeans 应该支持执行的混合模型。混合模型是 afex 包的一部分，他们提到应该支持混合对象。

有什么建议吗？

Answer 1

我无法重现您的问题。

根据 models supported by emmeans 的列表，afex 包中的 mixed 模型直接通过 afex 包得到支持。

我们可以使用 afex

中包含的样本数据集之一验证 mixed 模型拟合的边际均值计算

library(afex)
library(emmeans)

data(md_15.1)
model <- mixed(iq ~ timecat + (1 + time|id), data = md_15.1)

emmeans(model, "timecat")
# timecat emmean   SE   df lower.CL upper.CL
# 30         103 4.14 12.3     94.0      112
# 36         107 3.88 14.8     98.7      115
# 42         110 3.91 14.8    101.6      118
# 48         112 4.23 12.2    102.8      121
#
#Degrees-of-freedom method: kenward-roger
#Confidence level used: 0.95

对 class 'mixed' 的模型与 emmeans 进行 post-hoc 比较不起作用

Performing post-hoc comparisons with emmeans for model of class 'mixed' does not work

r

mixed-models

pairwise

emmeans