使用 Broom 包获得许多线性回归的 Anova 表中固定值的 P 值答案

【问题标题】：Obtain P-Value of Fixed Value in Anova Table of many Linear Regressions with Broom Package使用 Broom 包获得许多线性回归的 Anova 表中固定值的 P 值
【发布时间】：2021-11-24 06:08:46
【问题描述】：

在多元线性回归lm(FE_FCE2 ~ Trial + .x, data = DF_FCE3) 中有一个固定变量（试验）和许多 x 变量。我正在针对 FE_FCE2 分析每个 x 变量，并将试验作为固定效果。我比将繁荣包用于许多回归并将结果绘制在一张表中。我已经获得了回归结果的结果。但是无法将 ANOVA Table 中的数据添加到具有地图功能的 Broom 包中。

有可能吗？是的，如何？

我使用以下公式从回归结果中获取数据：

DF_FCE3 %>% 
  select(-FE_FCE2, -Trial) %>%  # exclude outcome, leave only predictors 
  map( ~lm(FE_FCE2 ~ Trial + .x, data = DF_FCE3)) %>% 
  map(summary) %>% 
  map_df(glance) %>% 
  round(3) -> rsme

但是，我想从 ANOVA Table of Trial 中获取 P 值 (**4.26e-08 *****)。到看看 Trial 是否对 x 变量有显着影响。

**$x1
Analysis of Variance Table

**Response: FE_FCE2
          Df  Sum Sq  Mean Sq F value   Pr(>F)    
Trial      3 0.84601 0.282002 15.0653 **4.26e-08 *****
.x         1 0.00716 0.007161  0.3826   0.5377    
Residuals 95 1.77827 0.018719**                     
---**

是否可以使用带有 map 功能的 broom 包来获得一个包含 anova 回归的所有许多 p 值的表？

【问题讨论】：

请创建一个reproducible example 或minimal reproducible example，其中包含可用于测试和验证可能解决方案的示例输入和所需输出。

标签： r regression purrr anova broom

【解决方案1】：

像这样（使用 mpg）？这将返回一个数据框，其中包含原始列和包含除结果和目标之外的 p 值的一行（在此示例中为 hwy 和 cyl，在您的情况下为 FE_FCE2 和 Trial）。

mpg %>%
  select(-hwy, -cyl) %>%  # exclude outcome, leave only predictors 
  map( ~lm(hwy ~ cyl + .x, data = mpg)) %>% 
  map(anova) %>% 
  map(broom::tidy) %>%
  map_df(~.$p.value[1])

【讨论】：