我无法将emmeans与glmer对象一起使用。我有一个不满足比例赔率假设的序数响应,我认为可以在每个切割点生成一个二进制模型作为替代方法。我最终想绘制出在每个响应级别(related question here)中的预测概率(对于变量组合)。
如果我跑步:
# Final clmm model
summary(mmod <- clmm(sc ~ visit + group + lang + (1|id), data = dat, nAGQ = 10))
# Emmeans
mmod_em <- emmeans(mmod, ~ sc + visit + group + lang, mode = "prob")
这很好用。但是,如果我再尝试使用简化的二进制模型之一,则emmeans无法识别因变量:
# Final glmer models
summary(mmod_2 <- glmer(sc_2 ~ visit + group + lang + (1|id), data = dat, family = "binomial", nAGQ = 10))
# Emmeans
mmod_2_em <- emmeans(mmod_2, ~ sc_2 + visit + group + lang)
Error in emmeans(mmod_2, ~sc_2 + visit + group + :
No variable named sc_2 in the reference grid
有什么想法吗?