Симуляция "zelig style" для многоуровневого GLMER в р
Я запускаю логистическую регрессию со смешанными эффектами с помощью r. Регрессия как-то так:
glmer (Y ~ X1 + X2 + X1: X2 + (1 | страна), data = hdp, family = binomial)
Теперь, с фиксированными эффектами, я хотел бы построить прогнозируемые вероятности Y. Я попробовал с Zelig, поскольку это то, что я узнал, как самый простой способ делать симуляции и получать прогнозируемые вероятности, но я видел, что новая версия не включает многоуровневую Модели Zelig и бывшего Multilevel очень "нестабильны". Есть ли легкая альтернатива? Как я могу делать симуляции, которые можно построить?
Заранее спасибо!