R：如何在没有此警告消息的情况下获得置信区间？答案

【问题标题】：R: How to get confidence interval without this warning message?R：如何在没有此警告消息的情况下获得置信区间？
【发布时间】：2017-04-29 22:12:39
【问题描述】：

型号是：

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

对于这个特定的模型：

是用confint还是其他功能更好？
如何解决该警告？
我能否获得可靠的置信区间？如何获得？

提前致谢

【问题讨论】：

置信区间是多少？模型中的每个参数都需要一个吗？
是的，我忘了指定
Here you go。读一读
@Sotos 谢谢。我真的不知道该怎么办。对于这种特定情况，您有什么解决方案吗？
查看使用arm 包的帖子的第三个答案。很容易理解

标签： r regression logistic-regression confidence-interval interpretation

【解决方案1】：

您始终可以在glm 中计算置信区间，而无需依赖任何类型的命令：

exp(confint.default(model))

您始终可以使用 Sotos 推荐的贝叶斯方法。取决于你想做什么。如果这就像一个硬件问题，告诉你只做一个 glm 模型和置信区间，那么上面的命令会帮助你。

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

【讨论】：

我注意到confint.default 不会产生警告。 confint 和 confint.default 有什么区别？
默认方法假定渐近正态性。 lm 对象中常见的东西。我强烈推荐使用?confint 并从文档中获取所有好东西。