如果您需要反复进行,请考虑以下事项:
require(tidyverse)
require(Hmisc)
构建功能:
flattenCorrMatrix
DF % as.matrix() %>% Hmisc::rcorr()
ut
flat
cor = (DF$r)[ut], p = DF$P[ut], n = DF$n[ut])
return(flat)
}
运行数据(我使用mtcars数据集):
mtcars %>%
flattenCorrMatrix() %>%
filter(cor < abs(0.5))
输出:
row column cor p n
1 mpg cyl -0.85216196 6.112688e-10 32
2 mpg disp -0.84755138 9.380328e-10 32
3 mpg hp -0.77616837 1.787835e-07 32
4 cyl drat -0.69993811 8.244636e-06 32
5 disp drat -0.71021393 5.282022e-06 32
6 hp drat -0.44875912 9.988772e-03 32
7 mpg wt -0.86765938 1.293958e-10 32
8 drat wt -0.71244065 4.784260e-06 32
9 mpg qsec 0.41868403 1.708199e-02 32
10 cyl qsec -0.59124207 3.660533e-04 32
11 disp qsec -0.43369788 1.314404e-02 32
12 hp qsec -0.70822339 5.766253e-06 32
13 drat qsec 0.09120476 6.195826e-01 32
14 wt qsec -0.17471588 3.388683e-01 32
15 cyl vs -0.81081180 1.843018e-08 32
16 disp vs -0.71041589 5.235012e-06 32
17 hp vs -0.72309674 2.940896e-06 32
18 drat vs 0.44027846 1.167553e-02 32
19 wt vs -0.55491568 9.798492e-04 32
20 cyl am -0.52260705 2.151207e-03 32
21 disp am -0.59122704 3.662114e-04 32
22 hp am -0.24320426 1.798309e-01 32
23 wt am -0.69249526 1.125440e-05 32
24 qsec am -0.22986086 2.056621e-01 32
25 vs am 0.16834512 3.570439e-01 32
26 mpg gear 0.48028476 5.400948e-03 32
27 cyl gear -0.49268660 4.173297e-03 32
28 disp gear -0.55556920 9.635921e-04 32
29 hp gear -0.12570426 4.930119e-01 32
30 wt gear -0.58328700 4.586601e-04 32
31 qsec gear -0.21268223 2.425344e-01 32
32 vs gear 0.20602335 2.579439e-01 32
33 mpg carb -0.55092507 1.084446e-03 32
34 disp carb 0.39497686 2.526789e-02 32
35 drat carb -0.09078980 6.211834e-01 32
36 wt carb 0.42760594 1.463861e-02 32
37 qsec carb -0.65624923 4.536949e-05 32
38 vs carb -0.56960714 6.670496e-04 32
39 am carb 0.05753435 7.544526e-01 32
40 gear carb 0.27407284 1.290291e-01 32