模型摘要中有和没有优势比的逻辑回归结果

Question

我喜欢将带和不带指数系数的逻辑回归与

modelsummary

包并排展示。该包产生了很好的 html 输出。它带有一个简单的选项来打开/关闭

exponentiate = TRUE

。但该选项适用于列表中的所有模型。

log_model = glm(vs ~ am, family = binomial(link = 'logit'), data = mtcars)

library(modelsummary)

modelsummary(list("No Odds Ratio" = log_model), exponentiate = FALSE)
modelsummary(list("Odds Ratio" = log_model), exponentiate = TRUE)

Answer 1

每个

modelsummary

函数调用返回一个字符串，其中包含表示 HTML/CSS 的“kableExtra”类。我们可以添加一个适当的 div 来获得你想要的东西。

log_model <- glm(vs ~ am, family = binomial(link = "logit"), data = mtcars)

library(modelsummary)

structure(
  paste(
    "<div style='display:flex; flex-wrap: wrap;'>",
    "<div style='flex-basis: 50%;'>",
    modelsummary(list("No Odds Ratio" = log_model), exponentiate = FALSE),
    "</div>",
    "<div style='flex-basis: 50%;'>",
    modelsummary(list("Odds Ratio" = log_model), exponentiate = TRUE),
    "</div>",
    "</div>",
    collapse = "\n"
  ),
  class = "kableExtra"
)

Answer 2

exponentiate

1.3.0 中的

modelsummary

参数接受一个逻辑向量来打开和关闭不同模型的开关。这在

?modelsummary

中有详细记录。例如：

library(modelsummary)
mod = glm(am ~ mpg, data = mtcars, family = binomial)
mod = list("OR" = mod, "Coef" = mod)
modelsummary(mod, exponentiate = c(TRUE, FALSE))

	或	系数
（拦截）	0.001	-6.604
	(0.003)	(2.351)
每加仑	1.359	0.307
	(0.156)	(0.115)
Num.Obs.	32	32
工商总局	33.7	33.7
比克	36.6	36.6
Log.Lik.	-14.838	-14.838
F	7.148	7.148
均方根误差	0.39	0.39