小编beb*_*bru的帖子

我想在一个有点大的数据集上计算 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")

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

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