predict.glmmTMB() 具有随机效应和 beta 误差族 模型

编程问答 2022-05-12

问题描述

我在使用 predict.glmmTMB() 时遇到了一些麻烦，当我从我的模型中去掉随机效应时，我可以让它正常工作，但是一旦我有了它们，它就拒绝工作.

模型

MOD10_temp<-glmmTMB (Adj_Prop ~  AvgGorg + TOD
                + (1|Site_ID),family=beta_family(link = "logit"),data = temp)

然后我设置了预测组件：

new_data_temp = data.frame(AvgGorg = rep(seq(0,50,1),2),TOD = c(rep(c('midday'),51),rep(c('Crepuscular'),51)))

# then the prediction
predict(MOD10_temp,type = 'response',newdata = new_data_temp)

但是我总是收到这个错误：

Error in eval(predvars,data,env) : object 'Site_ID' not found

所以在这个模型中，Site_ID 是随机效应，但我不明白为什么 predict() 函数不能用具有随机效应的模型进行预测？使用 glmer() 进行预测时我没有这个问题，但在这种情况下我需要一个 beta 分布，因此 glmmTMB。

非常感谢任何想法/帮助！

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

glmmtmbpredictrrrandom-effects