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 模型和置信区间，那么上面的命令将帮助你。

此外，如果您glm()仅使用一个参数运行 , 然后使用该confint命令，您将不会收到任何警告。

于 2016-12-14T15:24:19.510 回答