r - r中GLMER多级的模拟“zelig风格”

Question

我用 r 运行逻辑混合效应回归。回归在某种程度上是这样的：

glmer ( Y~ X1 + X2 + X1:X2 + (1 | country), data = hdp, family = binomial)

现在，使用固定效应，我想绘制 Y 的预测概率。我尝试使用 Zelig，因为这是我学到的最简单的模拟和获得预测概率的方法，但我看到新版本不包括多级模型和以前的 Zelig Multilevel 非常“不稳定”。有没有简单的替代方案？如何进行可以绘制的模拟？

提前致谢！

Answer 1

您可以使用merTools包。