“模拟”；“泽利格风格”；对于r中的GLMER多级

我用r。回归是这样的：

glmer（Y~X1+X2+X1:X2+（1 |国家），数据=hdp，家庭=二项式）

现在，在固定效应下，我想绘制Y的预测概率。我尝试了Zelig，因为这是我学到的最简单的模拟和预测概率的方法，但我看到新版本不包括多层次模型，以前的Zelig多层次模型非常“不稳定”。有什么简单的选择吗？我如何进行可以绘制的模拟

提前谢谢

您可以使用该软件包。

您是否尝试过

模拟（model，newdata=hdp[，-1]）

或

预测

？