我正在使用 glmmTMB 包运行广义线性混合效应模型来分析野生捕食者的咬伤数量。
我的模型:
Mod<-glmmTMB(all~Sex+class+sseason+timeR+class:timeR+(1|Microchip)+(1|olre), data = mydata_select, family = nbinom1, REML = FALSE, na.action = na.fail)
我的数据过于分散,因此我添加了观察级别的随机效应。我想使用预测功能根据模型预测咬伤的数量。我知道对于 glmer 模型,我可以使用 re.form=~(1|Microchip) 来确保它不包含观察级别的随机效应。但是,当为 glmmTMB 使用相同的公式时,我收到此消息:
mydata_select$predictAll<-predict(Mod, type= "response",re.form=~(1|Microchip))
Error in predict.glmmTMB(SelMod2, type = "response", re.form = ~(1 | Microchip)) :
re.form must equal NULL, NA, or ~0
我如何指定只考虑微芯片随机效应而不考虑观察水平随机效应?
欢迎任何建议。我确信必须有一个简单的方法来解决这个问题,但我一直无法找到它。
干杯,苏菲