使用 rcorr 的相关矩阵中的 p=0 是什么意思？答案

【问题标题】：What do the p=0 mean in a correlation matrix using rcorr?使用 rcorr 的相关矩阵中的 p=0 是什么意思？
【发布时间】：2015-05-27 15:35:22
【问题描述】：

我一直在使用 R 中 Hmisc 包中的 rcorr 函数。我弄清楚了它是如何工作的，我可以提取结果，将它们放在一个漂亮的表中。我看了这篇文章：p-values of correlation coefficients。但是......剩下的一件事困扰着我：如果 p 值 = 0.0000 怎么办？即使您使用打印语句添加更多数字，p-value = 0。

print(M_rcorr$P, digits = 20)

这是否意味着这是一个非常显着的相关性？我的意思是，R 的计算能力是否已经达到极限，这就是为什么它只会返回“0”？

例如下面的一些结果。

P
              Calcification Collagen Atheroma IPH    Macrophages Mastcells Neutrophils SMCs   Vessels
Calcification               0.0000   0.3390   0.0000 0.6488      0.0000    0.0063      0.0000 0.0000 
Collagen      0.0000                 0.0000   0.8866 0.7528      0.0000    0.0044      0.0000 0.0056 
Atheroma      0.3390        0.0000            0.0000 0.0000      0.7003    0.0000      0.0000 0.0063 
IPH           0.0000        0.8866   0.0000          0.0000      0.0474    0.0000      0.0000 0.0000 
Macrophages   0.6488        0.7528   0.0000   0.0000             0.0000    0.5536      0.0000 0.0000 
Mastcells     0.0000        0.0000   0.7003   0.0474 0.0000                0.0000      0.0001 0.0000 
Neutrophils   0.0063        0.0044   0.0000   0.0000 0.5536      0.0000                0.0000 0.0000 
SMCs          0.0000        0.0000   0.0000   0.0000 0.0000      0.0001    0.0000             0.0000 
Vessels       0.0000        0.0056   0.0063   0.0000 0.0000      0.0000    0.0000      0.0000

当我用“test$P”进一步检查时，我会得到这个：

              Calcification     Collagen     Atheroma          IPH  Macrophages    Mastcells
Calcification            NA 3.552714e-15 3.389508e-01 6.436717e-06 6.487711e-01 8.083783e-10
Collagen       3.552714e-15           NA 0.000000e+00 8.865855e-01 7.527929e-01 3.092686e-05
Atheroma       3.389508e-01 0.000000e+00           NA 0.000000e+00 3.623466e-09 7.003210e-01
IPH            6.436717e-06 8.865855e-01 0.000000e+00           NA 1.728324e-05 4.742716e-02
Macrophages    6.487711e-01 7.527929e-01 3.623466e-09 1.728324e-05           NA 1.545497e-11
Mastcells      8.083783e-10 3.092686e-05 7.003210e-01 4.742716e-02 1.545497e-11           NA
Neutrophils    6.319943e-03 4.367848e-03 1.532108e-14 0.000000e+00 5.535871e-01 0.000000e+00
SMCs           3.185481e-05 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 1.081798e-04
Vessels        1.372520e-05 5.565543e-03 6.339438e-03 1.666998e-06 0.000000e+00 0.000000e+00

谢谢！

砂光机

【问题讨论】：

help("format.pval")

标签： r correlation p-value hmisc

【解决方案1】：

好的。我想到了。当 p 值非常、非常、非常、非常小，小到 R 无法再计算时，p 值将为零。这是答案的实际部分，它只是意味着它在计算上不再可能。那就是说。要得到一个完全为零的 p 值可能在统计上/数学上是不可能的。

谢谢你的想法，伙计们！

【讨论】：

【解决方案2】：

这是一个相当统计的问题，但我尝试回答希望对您有所帮助。如果你有一个在原假设下不可能的样本，你肯定会得到一个 p 值为零。要根据 p 值解释结果，请看这里

http://www.graphpad.com/guides/prism/6/statistics/index.htm?stat_interpreting_results_correlati.htm

或者这里

http://www.eecs.qmul.ac.uk/~norman/blog_articles/p_values.pdf

【讨论】：

相关性永远不会得到完全为零的 p 值。