Ошибка: (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
не является целым числом