使用 purrr、broom 的许多单变量模型的整齐输出答案

【问题标题】：Tidy output from many single-variable models using purrr, broom使用 purrr、broom 的许多单变量模型的整齐输出
【发布时间】：2017-02-09 14:49:58
【问题描述】：

我有一个数据框，其中包含一个二元结果列 (y) 和多个独立预测变量列 (x1、x2、x3...)。

我想运行许多单变量逻辑回归模型（例如y ~ x1、y ~ x2、y ~ x3），并为每个模型提取指数系数（优势比）、95% 置信区间和 p 值模型成数据框/小标题的行。在我看来，应该可以使用 purrr 和 broom 的组合来解决。

This 问题类似，但我无法确定以下步骤：

仅提取我需要的值并
正在整理数据框/tibble。

从引用问题中的示例开始：

library(tidyverse)
library(broom)

df <- mtcars

df %>%
 names() %>%
 paste('am~',.) %>%
 map(~glm(as.formula(.x), data= df, family = "binomial"))

【问题讨论】：

标签： r purrr broom

【解决方案1】：

睡在上面之后，我想到了解决方案。需要使用map_df 运行每个模型，并使用tidy 从每个模型中提取值。

希望这对其他人有用：

library(tidyverse)
library(broom)

df <- mtcars

output <- df %>%
  select(-am) %>%
  names() %>%
  paste('am~',.) %>%
  map_df(~tidy(glm(as.formula(.x), 
               data= df, 
               family = "binomial"), 
               conf.int=TRUE, 
               exponentiate=TRUE)) %>%
  filter(term !="(Intercept)")

【讨论】：