У меня есть линейная смешанная модель lme (из пакета nlme) с расстоянием (log 10 преобразовано) в качестве переменной отклика, факторной переменной до или после перемещения в качестве предиктора, время log10 преобразовано как фиксированный эффект смещения, а затем два вложенные случайные эффекты ID (имена), вложенные в популяцию (фактор с двумя уровнями), а затем пол (мужской или женский), вложенный в сезон (размножение или разведение)
m.lme2 <- lme(l.S_Dists ~ PrePost_Trans + offset(l.St), random=list(Population=pdDiag(~ ID), Season= pdDiag(~Sex)), data = Q1a)
Я выполняю этот штраф - и проверяю AI C, что составляет 3562,669
Пока все в порядке.
I затем go, чтобы проверить этот AI C с нулевой моделью, чтобы увидеть, лучше ли тронется m.lme2. Используя код:
m.lme0 <- lme(l.S_Dists ~ 1 + offset(l.St), random=list(Population=pdDiag(~ ID), Season= pdDiag(~Sex)), data = Q1a)
Однако, когда я пытаюсь запустить эту модель, я получаю сообщение об ошибке:
Ошибка в lme.formula (l.S_Dists ~ 1 + смещение (l.St), random = список (Population = pdDiag (~ ID),: проблема nlminb, код ошибки конвергенции = 1 сообщение = единственная конвергенция (7)
Я уже запускал эту модель для очень похожий набор данных и не получает ошибок.
Может кто-нибудь сказать, почему я получаю это сообщение об ошибке при попытке запустить нулевую модель?
Вот как выглядит сводка для набора данных Q1a. Загрузите ссылки на изображения этого. Сводка набора данных Q1a . Сводка набора данных Q1a