我有一个鸟类数据框架 . 我有参与者的身份证号码,他们计算的鸟类数量,他们计算的年份,他们的纬度和长度坐标以及他们的努力 . 我做了这个模型:
model = lmer(count~year+lat+long+effort+(1|participant), data = df)
我现在希望模型从同一数据集中绘制预测值 . 因此,该数据是针对1997 - 2017年的,我希望该模型能够给出每年的预测值 . 我想绘制这些图,因此最终的图将在y轴上具有预测计数,在x轴上具有年(分类) . 每年将有一个数据点w /置信区间 .
我试图搞清楚 predict()
,但我有一个新的数据集来运行模型来预测未来的数量 . 我希望模型返回并根据 summary(model)
输出中的Beta值处理我之前已经放入的数据 .
我找到了这个线程,它似乎基本上是我 sjPlot
依赖下载, sjlabelled
每次都抛出一个错误:How to plot predicted values with standard errors for lmer model results?
2 回答
您可以尝试ggeffects-package,它将在即将到来的 sjPlot -update中用于绘制预测值 .
如果您缺少依赖项,请尝试:
您甚至可能希望安装ggeffects from GitHub,因为当前的dev-version对混合模型有一些修复和改进 .
我找到了我正在寻找的包,它被称为
predictedmeans
并且有一个函数,你可以放入模型和模型术语,你想要预测predictmeans(model, model term)
. 它完美无缺!