3

我一直在使用 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

谢谢！

桑德

4

2 回答 2

2

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

谢谢你的想法，伙计们！

于 2015-03-25T11:20:46.510 回答
0

这是一个相当统计的问题，但我尝试回答希望对您有所帮助。如果你有一个在原假设下不可能的样本，你肯定会得到零值。要根据 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

于 2015-03-24T12:48:35.273 回答