小编use*_*522的帖子

glmer - 用二项式数据预测(cbind计数数据)

我试图预测在我的二项式数据上运行的glmer模型随时间推移的值(x轴中的天数).Total Alive和Total Dead是计数数据.这是我的模型,以及下面的相应步骤.

full.model.dredge<-glmer(cbind(Total.Alive,Total.Dead)~(CO2.Treatment+Lime.Treatment+Day)^3+(Day|Container)+(1|index),
                         data=Survival.data,family="binomial")
Run Code Online (Sandbox Code Playgroud)

正如您在代码中看到的那样,我们已经考虑了过度离散(1:索引).

然后我们使用dredge命令确定具有主效应的最佳拟合模型(CO2.Treatment,Lime.Treatment,Day)及其相应的相互作用.

dredge.models<-dredge(full.model.dredge,trace=FALSE,rank="AICc")
Run Code Online (Sandbox Code Playgroud)

然后为它们创建了一个工作区变量

my.dredge.models<-get.models(dredge.models)
Run Code Online (Sandbox Code Playgroud)

然后,我们进行了模型平均,以对最佳拟合模型的系数求平均值

silly<-model.avg(my.dredge.models,subset=delta<10)
Run Code Online (Sandbox Code Playgroud)

但是现在我想要创建一个图形,Y轴上的Total Alive和X轴上的Days,以及取决于模型输出的拟合线.我知道这很棘手,因为模型连接了Total.Alive和Total.Dead(参见cbind(Total.Alive,Total.Dead)模型).

当我尝试运行预测命令时,我收到错误

# 9: In UseMethod("predict") :
#   no applicable method for 'predict' applied to an object of class "mer"
Run Code Online (Sandbox Code Playgroud)

r predict lme4 glm

6
推荐指数
1
解决办法
2万
查看次数

标签 统计

glm ×1

lme4 ×1

predict ×1

r ×1