Ошибка: (maxstephalfit) разделение на шаги PIRLS не смогло уменьшить отклонение в pwrssUpdate в lme4 для данных отношения

Мои данные - данные соотношения. Поэтому я пытаюсь использовать lme4 с binomail для анализа.

Вот мой код:

fate.reP = glmer(predated~type+(1|island),data=fate.rate,family="binomial")

Вот примерный набор данных:

type    cluster  tree   predated
 B        B7-1    1       0.48  
 B        B7-1    2       0.66
 B        B7-2    3       0.18
 M         I63    8       0.55
 M         I63    9       0.6
 M         I63   20       0.41
 M         I63   21       0.42
 S         I14    5       0.75
 S         I14   17       0.53
 S         I15    6       0.23
 S         I15    7       0.03

Когда я запускаю модель, она показывает, что:

Error: (maxstephalfit) PIRLS step-halvings failed to reduce deviance in pwrssUpdate
In addition: Warning message:
In eval(expr, envir, enclos) : non-integer #successes in a binomial glm!

Есть ли ошибки для моих данных или другие ошибки? спасибо большое я использую r 3.0.3 для windows

1 ответ

Решение

Предупреждение non-integer #successes in a binomial glm! дает вам подсказку, что биномиальные ответы в lme4 должны быть целыми числами. Я нигде не вижу знаменателей (т.е. общего числа людей, подверженных хищничеству): если они есть в ваших данных как (например) total_exposed, ты можешь использовать

fate.reP <- glmer(predated~type+(1|island),
     data=fate.rate,family="binomial",
     weights=total_exposed)

Поведение lme4 немного отличается от glm в базе R, которая предупредит вас, но все равно даст результаты, если glm не является целым числом

Другие вопросы по тегам