如何获得rcorr()中p小于0.05, r大于一定数值的数据

最新推荐文章于 2024-10-03 08:59:23 发布

梁韵750

最新推荐文章于 2024-10-03 08:59:23 发布

阅读量2.4k

点赞数

文章标签： r语言开发语言

本文链接：https://blog.csdn.net/weixin_60409733/article/details/121951327

版权

该篇博客介绍了如何在R语言中使用rcorr函数计算相关矩阵，并通过代码展示如何筛选出p值小于0.05以及相关系数大于0.75的变量对，以便于分析显著且高度相关的数据。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

如果一个数据集包含了较多的变量，在计算相关性之后会产生庞大的相关矩阵，而我们关心的可能只是其中显著相关（p<0.05）并且相关度比较高的数据，为了方便找到这些数据，以下code仅作参考：

library(Hmisc)

#flattenCorrMatrix在之后会用到，将相关矩阵转换成为一个每行由两两变量对应的“cor”和“p”组成的data.frame，参考网址Correlation matrix : Formatting and visualization - Easy Guides - Wiki - STHDA

flattenCorrMatrix <- function(cormat, pmat) {
ut <- upper.tri(cormat)
data.frame(
row = rownames(cormat)[row(cormat)[ut]],
column = rownames(cormat)[col(cormat)[ut]],
cor =(cormat)[ut],
p = pmat[ut]
)
}
result<-rcorr(as.matrix(data))
cor_p<-flattenCorrMatrix(result$r, result$P)

#首先找出 p<0.05的数据
cor_p.005 <