在R中的corrplot中如何计算显着性水平?

时间:2019-06-25 15:44:19

标签: r correlation significance

我正在调查13个类别变量的相关因子。

为了能够建立相关矩阵,我根据变量的级别将1或0赋予每个变量。然后,我构造了一个矩阵并计算了相关因子。 我已经使用了corrplot软件包。我已经使用以下代码计算了显着性水平。

cor.mtest <- function(mat, ...) {
  mat <- as.matrix(mat)
  n <- ncol(mat)
  p.mat<- matrix(NA, n, n)
  diag(p.mat) <- 0
  for (i in 1:(n - 1)) {
    for (j in (i + 1):n) {
      tmp <- cor.test(mat[, i], mat[, j], ...)
      p.mat[i, j] <- p.mat[j, i] <- tmp$p.value
    }
  }
  colnames(p.mat) <- rownames(p.mat) <- colnames(mat)
  p.mat
}
p.mat <- cor.mtest(matrice)

我不确定我是否确切了解每个相关性的显着性水平是如何计算的。我期望最高的相关因素是最重要的因素,但事实并非如此。例如,B和C之间的相关性是5个最高相关因子之一,但没有被标记为具有统计学意义(p值0.05)。

使用上述代码获得的具有统计意义的相关性为:C-A,D-N,D-J,E-N,E-J,I-D,M-N,H-I,J-A,E-D,I-A,I-J和K-D。

有人可以帮助我理解为什么吗?

这是我的矩阵:

           A          B           C           D          E          F          G          H          I           J          K          L          M          N
A          NA         -0.204124145-0.327253796-0.17393886-0.11047384 0.17770466 0.07747516-0.04936241 0.042170673 0.27964508-0.03084243 0.45412438 0.27273783 0.26841809
B          -0.20412415NA          -0.240481445-0.03873286-0.08118136-0.03627381 0.05693239-0.03627381 0.009390603-0.05708231 0.08426531 0.09004924 0.26674678 0.07927153
C           -0.3272538-0.240481445NA            0.1248002 0.06544512-0.05815453 0.12407634  0.1508383 0.004234248-0.10467018-0.04236275 -0.2065999 0.09535462  0.0679832
D          -0.17393886 -0.03873286 0.124800202NA          0.26119846 0.13874261 0.08866557-0.06322448-0.171860016-0.28065388 0.26335782 -0.2854407-0.17731715-0.22791043
E          -0.11047384-0.081181365 0.065445119 0.26119846NA         -0.06067985 0.15966387 0.14456083 0.074155081-0.26666679 0.28563978-0.08913118 -0.1248569-0.22093625
F           0.17770466-0.036273813-0.058154526 0.13874261-0.06067985NA           0.2248724-0.00877193-0.023844308-0.04555305 0.19815355-0.05187858 0.06450615 0.07104128
G           0.07747516 0.056932385  0.12407634 0.08866557 0.15966387  0.2248724NA           0.2248724 0.303847721-0.14028442 0.12021134 0.05830386 0.12951876-0.04026408
H          -0.04936241-0.036273813 0.150838301-0.06322448 0.14456083-0.00877193  0.2248724NA          0.367883604-0.04555305-0.04426835 0.16908576-0.13598594-0.12347651
I           0.04217067 0.009390603 0.004234248-0.17186002 0.07415508-0.02384431 0.30384772  0.3678836NA          -0.03135818-0.02619483 0.11639665  0.0196331-0.03349833
J           0.27964508-0.057082312-0.104670181-0.28065388-0.26666679-0.04555305-0.14028442-0.04555305-0.031358183NA          -0.1154414 0.35648931-0.09094761 0.09342798
K          -0.03084243 0.084265305-0.042362748 0.26335782 0.28563978 0.19815355 0.12021134-0.04426835-0.026194825 -0.1154414NA         -0.04940683-0.01173103-0.15568226
L           0.45412438 0.090049245-0.206599899 -0.2854407-0.08913118-0.05187858 0.05830386 0.16908576  0.11639665 0.35648931-0.04940683NA          0.11800138 0.16450294
M           0.27273783 0.266746783 0.095354617-0.17731715 -0.1248569 0.06450615 0.12951876-0.13598594 0.019633103-0.09094761-0.01173103 0.11800138NA          0.52140784
N           0.26841809 0.079271526 0.067983203-0.22791043-0.22093625 0.07104128-0.04026408-0.12347651-0.033498334 0.09342798-0.15568226 0.16450294 0.52140784NA

0 个答案:

没有答案