GlmmNP包给我错误

Question

我是glmm的新手。我要求它符合一个模型，该模型考虑到用两种花粉处理方法喂养的蜜蜂，放在不同的笼子里，并每天记录它们的体重。我试图解决的问题是，蜜蜂的体重是否随饮食而变化。笼子，复制体被嵌套在复制中。我的响应变量是权重。我的解释变量是:饮食、日期和复制/笼子/蜜蜂。需要两个模型，一个包含随机变量的混合效应，第二个。我运行的代码如下所示：

使用带有函数glmmNP的gamlss.mx包

库(gamlss.mx)库(gamlss.dist)

m1<- glmmNP(权重.~ Diet+ (1|复制/笼子/蜜蜂)，family=“威布尔”，data=beewt)

m1<- glmmNP(权重.~ Diet+ (1|复制/笼/蜜蜂)，family="Gamma"，data=beewt)

m1<- glmmNP(权重.~ Diet+ (1|复制/笼子/蜜蜂)，family="Gumbell"，data=beewt)

错误消息是:找不到函数glmmNP#

我尝试过其他软件包，比如MASS及其功能: glmm(PQL)库(MASS)模型<- glmmPQL(Weight.g.~ Diet+ (1|Replicate/Cage/Bee)，family="Gamma"，data=beewt)

错误消息是:参数"random“缺失，没有默认值

请告诉我如何适当地调整我的代码。

Answer 1

怎么样

bee.pql <- lmer(Weight.g. ~ Diet+Cage+Day+(1|Bee), data=beewt)

这将体重视为一个连续的结果，一个蜜蜂的多个记录是相关的。