我有一个相关矩阵,我试图测试使用rcorr()得到相关矩阵的p值。这可以在Hmisc
包中找到。我有一个矩阵,你可以在这里看到KCN
,你可以看到这里的前6个记录:
0 0 0 0 1 0
0 0 0 0 1 0
0 0 1 0 2 0
2 1 5 3 5 0
0 0 1 0 0 0
0 2 1 0 2 0
然后我运行rcorr:
rKCR<-rcorr(as.matrix(KCN[,2:7]),type="pearson")
如果我然后调用rKCR,我会继续得到p值矩阵的所有0,并且我认为这些是非常弱的相关性,正在被舍入。我在KCN的两列上运行了corr.test()函数:
cor.test(KCN[,2],KCN[,3],type="pearson")
我的p值为2.2e-16,显然非常小。我在rcorr()函数中使用了整个数据框,在查看文档之后,我想知道是否有任何方法可以获得以科学计数法表示的p值矩阵?谢谢你的帮助。
答案 0 :(得分:0)
您可以使用:
k <- cor.test(KCN[,2],KCN[,3],type="pearson")
pVal <- k$p.value