r - glmer 的 emmeans 错误 - 在参考网格中无法识别变量

Question

我在使用带有 glmer 对象的 emmeans 时遇到问题。我有一个不满足比例优势假设的序数响应，我想在每个切点处生成一个二元模型作为替代方法。我最终想绘制每个响应级别的预测概率（对于变量组合）（相关问题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 

有任何想法吗？