Я пытаюсь вписать GLMM в R и получаю два предупреждения.Пробовал страницу ?help
для пакета / функции, а также ?convergence
, но на самом деле не решил.=> Я начал работать с распределением Пуассона, затем диагностика показала избыточную дисперсию и нулевое раздувание.Затем я перешел на отрицательное биномиальное распределение с glmer.nb
.Затем страница ?help
говорит, что если мы ЗНАЕМ, что есть чрезмерная дисперсия, мы должны использовать *glmer(family=MASS::negative.binomial(theta=1.75))*
вместо этого.
ВОПРОСЫ: 1) это значение тета (1.75) фиксируется формулой, или я должен найти тета для моей модели?2) Пробовал страницу ?convergence
, и не нашел ответа.Итак, кто-нибудь знает возможный способ решения этого предупреждения: "1: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model failed to converge with max|grad| = 0.172241 (tol = 0.001, component 1)
2: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model is nearly unidentifiable: very large eigenvalue
- Rescale variables?"
- Исходная переменная: количество «мертвых»;К независимым переменным относятся категориальные и три случайных эффекта.* Иерархические данные, поэтому каждый «орд» (единица) имеет уровни «источник» и «неделя».(250 или 7 недель * 3 источника = 5250 строк).
Если у кого-нибудь есть совет, я был бы очень признателен.
Мой код:
nbm <- glmer(dead ~ week + year + wean.month + classB.Q + euth.p.cat + source + protec + wean.cat.Q+ arr.euth.p.cat + (1|source:ord) + (1|ord) + (1|barngrp), offset=log(survival), data=final, family=MASS::negative.binomial(theta=1.75),control=glmerControl(optimizer = c("bobyqa", "Nelder_Mead"), optCtrl=list(maxfun=15000000)))