在“corrr”package中，如何为统计显着性添加星号

Question

我有一个相关矩阵，其中包括 14 个变量之间的双变量相关性。 我如何用 append 星号来表示统计显着性？ 我正在使用以下代码：

pretty.matrix<-a %>%  
  correlate() %>%  
  shave() %>% 
  fashion() %>% 
  print()

Answer 1

您可以使用 function 打印统计显着性。 使用colpair_map可以很容易地制作一个漂亮的矩阵。 我以 mtcars 数据集为例。 您可以使用以下代码：

library(corrr)
  library(tidyverse)
  
  # Function
  calc_p_value <- function(vec_a, vec_b, sig_level){
    test_res <- cor.test(vec_a, vec_b)
    sig <- if_else(test_res$p.value < sig_level, "*", "")
    paste0(round(cor.test(vec_a, vec_b)$estimate, 2), sig)
  }
  
  # Matrix with p = 0.05
  colpair_map(mtcars, calc_p_value, 0.05) %>%
    shave()

Output 看起来像这样：

# A tibble: 11 × 12
   term  mpg    cyl    disp   hp     drat   wt     qsec   vs     am    gear  carb 
   <chr> <chr>  <chr>  <chr>  <chr>  <chr>  <chr>  <chr>  <chr>  <chr> <chr> <chr>
 1 mpg   NA     NA     NA     NA     NA     NA     NA     NA     NA    NA    NA   
 2 cyl   -0.85* NA     NA     NA     NA     NA     NA     NA     NA    NA    NA   
 3 disp  -0.85* 0.9*   NA     NA     NA     NA     NA     NA     NA    NA    NA   
 4 hp    -0.78* 0.83*  0.79*  NA     NA     NA     NA     NA     NA    NA    NA   
 5 drat  0.68*  -0.7*  -0.71* -0.45* NA     NA     NA     NA     NA    NA    NA   
 6 wt    -0.87* 0.78*  0.89*  0.66*  -0.71* NA     NA     NA     NA    NA    NA   
 7 qsec  0.42*  -0.59* -0.43* -0.71* 0.09   -0.17  NA     NA     NA    NA    NA   
 8 vs    0.66*  -0.81* -0.71* -0.72* 0.44*  -0.55* 0.74*  NA     NA    NA    NA   
 9 am    0.6*   -0.52* -0.59* -0.24  0.71*  -0.69* -0.23  0.17   NA    NA    NA   
10 gear  0.48*  -0.49* -0.56* -0.13  0.7*   -0.58* -0.21  0.21   0.79* NA    NA   
11 carb  -0.55* 0.53*  0.39*  0.75*  -0.09  0.43*  -0.66* -0.57* 0.06  0.27  NA