Пакет GlmmNP выдает ошибки

Question

Я новичок в glmm. Я требую, чтобы он соответствовал модели, которая учитывала бы пчел, которых кормили двумя способами обработки пыльцы, помещенных в разные клетки, и их вес регистрировался ежедневно. Я хочу ответить на вопрос, меняется ли вес пчел в зависимости от рациона. Клетки, реплики были вложены в Replicate. Моя ответная переменная - вес. Мои объясняющие переменные: Диета, День и Репликация / Клетка / Пчела. Требуются две модели: одна со смешанными эффектами, включая случайную величину, и вторая. Коды, которые я использую, показаны ниже:

ИСПОЛЬЗОВАНИЕ пакета gamlss.mx с его функцией glmmNP

библиотека (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)

Сообщение об ошибке: аргумент "случайный" отсутствует, значение по умолчанию отсутствует.

Посоветуйте, пожалуйста, как я могу соответствующим образом изменить свои коды.

Answer 1

Как насчет

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

это рассматривает вес как непрерывный результат, и несколько записей об одной пчеле связаны.