Мои данные заполнены нулями, поэтому я использую модель с нулевым раздувом, используя glmmamdb:
Model3z <- glmmadmb (Count3 ~ Light3 + (1 | Site3), zeroInflation = T, семейство= "Пуассон", данные = dframe3) </strong>
Однако, когда я пытаюсь выполнить попарные сравнения различных типов освещения в этой модели pwcs3 <- lsmeans (Model3z, "Light") </strong> Я получаю сообщение об ошибке:
Ошибка в ref_grid (объект, ...): не удается обработать объект класса «glmmadmb». Использовать справку ("models", package = "emmeans") дляинформация о поддерживаемых моделях.
Когда я захожу на сайт пакета emmeans, он говорит, что glmmadmb больше не поддерживается.
Я переключился на pscl и функцию zeroinfl, но не уверен, какреструктурировать мой код, чтобы соответствовать формату PSCL. При наборе P <- zeroinfl (Count3 ~ Light3 + (1 | Site3), family = poisson, data = dframe3) </strong> появляется сообщение об ошибке:
Ошибка в contrasts<- (*tmp*, value = contr.funs [1 + isOF [nn]]): контрасты могут применяться только к факторам с 2 или более уровнями. Дополнительно: Предупреждение: In Ops.factor (1, Site3): '|'не имеет значения для факторов
Есть ли другой способ использования glmmadmb с lsmeans? Если нет, кто-нибудь знает, как должен выглядеть код модели с нулевым раздувом в pscl? Спасибо.