如何解决predict.glmmTMB() 具有随机效应和 beta 误差族 模型
我在使用 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 举报,一经查实,本站将立刻删除。