微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

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

如何解决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 举报,一经查实,本站将立刻删除。