警告:非整数#成功的二项式glm！ (调查包) [英] Warning: non-integer #successes in a binomial glm! (survey packages)

问题描述

我正在使用twang包来创建倾向得分，该得分用作使用survey::svyglm的二项式glm中的权重.代码看起来像这样:

I am using the twang package to create propensity scores, which are used as weights in a binomial glm using survey::svyglm. The code looks something like this:

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!

有人知道我在做什么错吗?

Does anyone know what I could be doing wrong ?

我不确定stats.SE上的信息是否会更好，但总的来说，我认为我会先在这里尝试.

I wasn't sure if this message would be better on stats.SE, but on balance I thought I would try here first.

推荐答案

没有错，glm在指定二项式(和Poisson)模型时只是挑剔的.如果它检测到否，它会发出警告.尝试或成功的过程不是完整的，但无论如何它还是可以适应模型的.如果您想取消该警告(并且确定这不是问题)，请改用family=quasibinomial.

There's nothing wrong, glm is just picky when it comes to specifying binomial (and Poisson) models. It warns if it detects that the no. of trials or successes is non-integral, but it goes ahead and fits the model anyway. If you want to suppress the warning (and you're sure it's not a problem), use family=quasibinomial instead.

这篇关于警告:非整数#成功的二项式glm！ (调查包)的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！