首页 文章

来自lmer模型的预测(?)值

提问于
浏览
0

我有一个鸟类数据框架 . 我有参与者的身份证号码,他们计算的鸟类数量,他们计算的年份,他们的纬度和长度坐标以及他们的努力 . 我做了这个模型:

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 回答

  • 0

    您可以尝试ggeffects-package,它将在即将到来的 sjPlot -update中用于绘制预测值 .

    library(ggeffects)
    dat <- ggpredict(model, terms = "dat")
    plot(dat)
    

    如果您缺少依赖项,请尝试:

    install.packages(
      c("sjlabelled", "sjmisc", "sjstats", "ggeffects", "sjPlot"),
      dependencies = TRUE
    )
    

    您甚至可能希望安装ggeffects from GitHub,因为当前的dev-version对混合模型有一些修复和改进 .

    devtools::install_github("strengejacke/ggeffects")
    
  • 0

    我找到了我正在寻找的包,它被称为 predictedmeans 并且有一个函数,你可以放入模型和模型术语,你想要预测 predictmeans(model, model term) . 它完美无缺!

相关问题