我正在尝试建立一个泊松混合模型,但残差有些异方差。所以,我想为模型定义一个方差结构。据我所知,我可以通过使用包glmmPQL
中的命令来做到这一点MASS
。但这由于某种原因不起作用。这是代码:
library(nlme)
library(MASS)
VS=varIdent(form=~1|Experience.group)
fit1<-glmmPQL(Total.pecks~Age+Sex+Trial.no+Experience.group,
random=~1|as.factor(Name),family=poisson,weights=VS)
我认为一切都很酷,但我收到了以下错误消息:
Error in model.frame.default(weights = VS, formula = ~Total.pecks + Age + :
variable lengths differ (found for '(weights)')
有任何想法吗?谢谢