0

我正在尝试建立一个泊松混合模型,但残差有些异方差。所以,我想为模型定义一个方差结构。据我所知,我可以通过使用包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)')

有任何想法吗?谢谢

4

0 回答 0