r - R 中带扫帚的 glm 的快速 Wald 置信区间

Question

我想在一个有点大的数据集上计算 glm 系数的 Wald 置信区间，并使用 broom为了一个整洁的输出。

mydata <- data.frame(y = rbinom(1e5,1,0.8), 
                 x1 = rnorm(1e5), 
                 x2 = rnorm(1e5))
glm.1 <- glm(y ~ x1 + x2, data = mydata, family = "binomial")

使用 broom::tidy在大数据上花费大量时间，因为它使用 confint.glm ，它根据分析的对数似然函数计算置信区间。

tidy(glm.1, conf.int = TRUE) # can take literally hours

Answer 1

confint和 confint.glm分别不为用于计算置信区间的方法提供参数。如果要使用另一种方法，则需要使用不同的函数，例如confint.default对于沃尔德。
broom::tidy反过来，使用的函数没有参数(或者我错过了什么？)，它总是调用 confint.glm为glm。

要使用不同的函数计算置信区间，broom有 confint_tidy ，您可以在其中指定要使用的函数:

confint_tidy(glm.1, func = stats::confint.default)

将其与估计值放在一起:

cbind(tidy(glm.1), confint_tidy(glm.1, func = stats::confint.default))