如何在R中的列表中插入p值低于5%的回归器的变量名

Question

我试图在R中构建一个列表，其中包含所有的回归者名称，这些名称的p值低于5%的阈值。例如：

第一回归

#gender (male - female)
regr1 <- lm(salary ~ female, data = test)
summary(regr1)

输出第一回归：

 Coefficients:
         Estimate Std. Error t value Pr(>|t|)    
 (Intercept)  0.855618   0.001888  453.24   <2e-16 ***
 female      -0.054514   0.003088  -17.65   <2e-16 ***

第二次回归：

#education (PhD - Master - Bachelor - HighSchool - None)
regr2 <- lm(salary ~ Master + Bachelor + HighSchool + None, data = test)
summary(regr2)

输出二次回归：

Coefficients:
        Estimate Std. Error t value Pr(>|t|)    
(Intercept)  0.91008    0.02829  32.164  < 2e-16 ***
Master      -0.05446    0.02836  -1.920 0.054811 .  
Bachelor    -0.10291    0.02848  -3.613 0.000303 ***
HighSchool  -0.10173    0.02911  -3.495 0.000475 ***
None        -0.12590    0.02864  -4.396 1.11e-05 ***

变量母版没有意义，所以我不希望它出现在列表中。这是我想要的清单：

varnames <- c("female", "Bachelor", "HighSchool", "None")  

Answer 1

您可以使用broom::tidy，然后操作表，如下所示：

library(tidyverse)

tab <- lm(data = mtcars, mpg ~ cyl + disp + hp) %>% summary() %>% broom::tidy()
tab
# A tibble: 4 x 5
  term        estimate std.error statistic  p.value
  <chr>          <dbl>     <dbl>     <dbl>    <dbl>
1 (Intercept)  34.2       2.59       13.2  1.54e-13
2 cyl          -1.23      0.797      -1.54 1.35e- 1
3 disp         -0.0188    0.0104     -1.81 8.09e- 2
4 hp           -0.0147    0.0147     -1.00 3.25e- 1

然后过滤p.value列：

tab %>% filter(p.value < 0.05)
# A tibble: 1 x 5
  term        estimate std.error statistic  p.value
  <chr>          <dbl>     <dbl>     <dbl>    <dbl>
1 (Intercept)     34.2      2.59      13.2 1.54e-13

因此，现在您可以使用回归者的名称：

tab %>% filter(p.value < 0.05) %>% select(term) %>% as.character()
[1] "(Intercept)"