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

Question

我用r运行逻辑混合效应回归。回归是这样的：

glmer（Y~X1 + X2 + X1：X2 +（1 | country），data = hdp，family = binomial）

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

提前致谢！

Answer 1

您可以使用merTools包。