56

パッケージを使用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 で良いかどうかはわかりませんでしたが、バランスを考えて、最初にここで試してみようと思いました。

4

4 に答える 4

77

二項（およびポアソン）モデルの指定に関してglmは、うるさいだけです。いいえを検出すると警告します。試行回数または成功回数は整数ではありませんが、先に進んでモデルに適合します。警告を抑制したい場合 (そしてそれが問題ではないと確信している場合) は、family=quasibinomial代わりに を使用してください。

于 2012-10-18T12:00:53.100 に答える