Предупреждение: нецелые #successes в биномиальном glm! (обзорные пакеты)

Я использую пакет twang для создания оценок склонности, которые используются в качестве веса в биномиальном glm с использованием survey::svyglm. Код выглядит примерно так:

pscore <- ps(ppci ~ var1+var2+.........., data=dt....)

dt$w <- get.weights(pscore, stop.method="es.mean")

design.ps <- svydesign(ids=~1, weights=~w, data=dt,)

glm1 <- svyglm(m30 ~ ppci, design=design.ps,family=binomial)

Это вызывает следующее предупреждение:

Warning message:
   In eval(expr, envir, enclos) : non-integer #successes in a binomial glm!

Кто-нибудь знает, что я могу делать неправильно?

Я не был уверен, будет ли это сообщение лучше на stats.SE, но на балансе я подумал, что сначала попробую.

Ответ 1

Нет ничего плохого, glm просто придирчиво, когда дело доходит до определения биномиальных (и пуассоновских) моделей. Он предупреждает, обнаруживает ли он, что нет. испытаний или успехов не является неотъемлемой частью, но он идет вперед и все равно соответствует модели. Если вы хотите подавить предупреждение (и вы уверены, что это не проблема), используйте family=quasibinomial вместо этого.

Ответ 2

Нет ничего плохого, вычислительно, но статистически вы не можете делать что-то, что имеет смысл. В таком случае, вероятно, лучше использовать метод надежной регрессии, который, как правило, является хорошей идеей для данных пропорционального ответа, если ваши данные включают единицы с ровно 1 или ровно 0.

Ответ 3

Извините, но он более устойчив в том смысле, что если базовый механизм является переопределенной биномиальной моделью, переопределенный бином будет учитывать его при оценке стандартного erorr. Следовательно, вы получите лучший охват, хотя оценки точек одинаковы.