0

我运行了以下混合模型:

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 包的一部分,他们提到应该支持混合对象。

有什么建议么?

4

1 回答 1

0

我无法重现您的问题。

根据包中支持的机型列表,机型emmeans 直接通过包支持。mixedafexafex

mixed我们可以使用模型拟合中包含的样本数据集之一来验证边际均值的计算afex

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
于 2019-08-07T08:13:23.343 回答