У меня проблемы с использованием emmeans с блеском объекта. У меня есть порядковый ответ, который не удовлетворяет предположению о пропорциональных коэффициентах, и я подумал создать двоичную модель в каждой точке отсечения в качестве альтернативного подхода. В конечном итоге я хочу построить прогнозируемые вероятности (для комбинации переменных) для каждого уровня ответа ( связанный вопрос здесь ).
Если я бегу:
# 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
Есть идеи?