我是新来的。我要求它拟合一个模型，该模型考虑了以两种花粉处理方式喂养的蜜蜂，放置在不同的笼子中，并每天记录它们的重量。我试图解决的问题是蜜蜂的体重是否随饮食而变化。笼子，复制品嵌套在复制品中。我的响应变量是体重。我的解释变量是：饮食、日和复制/笼子/蜜蜂。需要两个模型，一个具有混合效应，包括随机变量，第二个。我正在运行的代码如下所示：
使用带有功能 glmmNP 的 gamlss.mx 包
库（gamlss.mx） 库（gamlss.dist）
m1<- glmmNP(Weight.g.~ Diet+ (1|Replicate/Cage/Bee), family="weibull", data=beewt)
m1<- glmmNP(Weight.g.~ Diet+ (1|Replicate/Cage/Bee), family="Gamma", data=beewt)
m1<- glmmNP(Weight.g.~ Diet+ (1|Replicate/Cage/Bee), family="Gumbell", data=beewt)
错误信息是：找不到函数 glmmNP#
我已经尝试过其他软件包，如 MASS 的功能：glmm(PQL) library(MASS) model<- glmmPQL(Weight.g.~ Diet+ (1|Replicate/Cage/Bee), family="Gamma", data=beewt)
错误消息是：缺少参数“random”，没有默认值
请告知我如何适当地调整我的代码。