1. ホーム
  2. r

[解決済み】警告:二項式glmで非整数の#成功例! (調査パッケージ)

2022-02-22 23:08:50

質問内容

を使用しています。 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の方がいいのか迷いましたが、バランス的にまずこちらで試してみようと思いました。

どのように解決するのですか?

何も問題がない glm は、二項(およびポアソン)モデルを指定するときに、ただうるさいだけです。試行回数や成功回数が積分でないことを検出すると警告が出ますが、とにかくそのままモデルを当てはめます。この警告を消したい場合(そしてそれが問題でないと確信している場合)、次のようにします。 family=quasibinomial の代わりに