如何使用 R 获得数据框中每个列对的相关系数？ [复制]答案

【问题标题】：How to obtain correlation coefficients for every column pair in a dataframe with R? [duplicate]如何使用 R 获得数据框中每个列对的相关系数？ [复制]
【发布时间】：2021-04-04 21:37:51
【问题描述】：

这里有一些代码：

for(i in 1:colCount){
 for(j in colCount:i){
     cor(newCovidNoNaDF[,i],newCovidNoNaDF[,j], method = "pearson")
 }
}

但是当我点击进入时没有任何输出。我做错了什么？

【问题讨论】：

你试过print()吗？
cor(newCovidNoNaDF)
谢谢马丁。 print() 在我之前使用数据表时不起作用，所以我认为它不适用于数据框。我错了！

标签： r for-loop

【解决方案1】：

您可以使用combn。以mtcars 为例：

val <- combn(seq_along(mtcars), 2, function(x) 
             cor(mtcars[[x[1]]], mtcars[[x[2]]], method = 'pearson'))
name <- combn(names(mtcars), 2, paste0, collapse = '-')
data <- data.frame(name, val)

head(data)

#      name        val
#1  mpg-cyl -0.8521620
#2 mpg-disp -0.8475514
#3   mpg-hp -0.7761684
#4 mpg-drat  0.6811719
#5   mpg-wt -0.8676594
#6 mpg-qsec  0.4186840

【讨论】：

【解决方案2】：

corrr 包中的 correlate() 函数计算数据帧的成对列相关性：

library(corrr)

correlate(df)

输出：

Correlation method: 'pearson'
Missing treated using: 'pairwise.complete.obs'

# A tibble: 3 x 4
  term       a      b      c
  <chr>  <dbl>  <dbl>  <dbl>
1 a     NA      0.182  0.223
2 b      0.182 NA     -0.262
3 c      0.223 -0.262 NA

数据：

set.seed(123)
n <- 20
df <- tibble(a = rnorm(n), b = rnorm(n), c = rnorm(n))

【讨论】：