GlmmNP 包给了我错误

Question

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

使用 gamlss.mx 包及其函数 glmmNP

库（gamlss.mx） 库（gamlss.dist）

m1

m1

m1

错误信息是：找不到函数glmmNP#

我已经尝试过其他软件包，例如 MASS 的功能：glmm(PQL) 图书馆（大众） 模型

错误信息是：参数“random”丢失，没有默认值

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

Answer 1

怎么样

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

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