1

モデルは次のとおりです。

model <- glm(DW ~ P + DV_1, family = "binomial")

変数は次のとおりです。

DW <- c(1, 1, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 1, 1)
P <- c(18.584898, 8.177430, -7.392020, -13.123626, 11.742363, 35.836419, 8.177430, 8.177430, 7.209096, 10.398933, -23.382043, -8.177430, 7.392020, 17.607980, -37.631207, -8.177430, 12.202439, -29.602930, -8.177430, 14.709837, 8.194932, 8.177430, -5.222738, 1.185302, 12.049662, 6.193046)
DV_1 <- c(45.49215, 55.40000, 51.63815, 36.12306, 34.78324, 41.17867, 59.14783, 62.45898, 55.04072, 53.76998, 52.31764, 44.71056, 42.23566, 50.08676, 61.34397, 49.59538, 38.21099, 51.05214, 44.69676, 40.83045, 46.09846, 53.45508, 54.73643, 50.26476, 48.75601, 53.68885)

パラメータごとに信頼区間を取得しようとすると、confint次の警告メッセージが表示されます。

glm.fit: 0 または 1 が発生した確率を数値的に当てはめました

この特定のモデルの場合:

  1. confintまたは他の機能を使用する方が良いですか？
  2. その警告を修正するにはどうすればよいですか?
  3. 信頼できる信頼区間を取得できますか?またその方法は?

前もって感謝します

4

1 に答える 1

0

glmどのような種類のコマンドにも依存することなく、いつでも次のように信頼区間を計算できます。

exp(confint.default(model))

Sotos が推奨するベイジアン アプローチをいつでも使用できます。何をしたいかによって異なります。これが glm モデルと信頼区間だけを実行するように指示する HW の質問のようなものである場合は、上記のコマンドが役立ちます。

glm()また、パラメータを 1 つだけ指定してを実行してからconfintコマンドを使用しても、警告は表示されません。

于 2016-12-14T15:24:19.510 に答える