有沒有辦法用系數、se、p、優勢比和 CI 從 R 編寫回歸表？

Question

我有相當可靠的數據，並希望將我的回歸輸出從 R 寫入“期刊就緒”格式。

目前，我已經使用 outreg 包來獲取系數、se 和 p 值（下面的代碼）。 但是，我似乎無法找到如何還包括優勢比或置信區間。 有沒有人對這樣的包有建議，或者我如何使用 outreg 實現這一點？

上面提到的outreg表代碼：

dm1.output <- list(dm1a, dm1b, dm1c, dm1d, dm1e, dm1f, dm1g, dm1h)
dm1.output2 <- as.data.frame(outreg(dm1.output, pv=T, starred = c("pv")))

Answer 1

stargazer包非常適合您的任務。 它的設計目標是生成“期刊就緒”表格，因為開發人員是學者。

Answer 2

如果它們是標准模型，您可以在此處查看模型modelsummary ：

library(modelsummary)
models <- list("Model 1" <- glm(am ~ mpg, data = mtcars, family = binomial),
               "Model 2" <- glm(am ~ cyl, data = mtcars, family = binomial))
modelsummary(models, exponentiate = TRUE, stars = T, statistic = 'conf.int', conf_level = .95)

---

或stargazer ：

library(stargazer)   
model_1 <- glm(am ~ mpg, data = mtcars, family = binomial)
stargazer(model_1, coef = list(exp(model_1$coefficients)), type = "text")

---

或texreg

library(texreg)
model_trexreg <- texreg::extract(model_1)
screenreg(model_1, override.coef = exp(model_trexreg@coef), override.se = exp(model_trexreg@se))