r - 两个不同维度矩阵之间的相关性

Question

我对 R 非常陌生。我有两个不同维度的矩阵，C(3 行，79 列)和 T(3 行，215 列)。我希望我的代码计算 C 的第一列与 T 的所有列之间的 Spearman 相关性，并返回与索引和列的最大相关性。然后，C 的第二列和 T 的所有列，依此类推。事实上，我想找到两个矩阵之间最相关的列。希望它是清楚的。我所做的是一个嵌套的 for 循环，但结果不是我搜索的。

for (i in 1:79){
    for(j in 1:215){
        print(max(cor(C[,i],T[,j],method = c("spearman"))))
  }
}

Answer 1

您不必循环遍历各列。

x <- cor(C,T,method = c("spearman"))

out <- data.frame(MaxCorr = apply(x,1,max), T_ColIndex=apply(x,1,which.max),C_ColIndex=1:nrow(x))

head(out)

给予，

  MaxCorr T_ColIndex C_ColIndex
1       1          8          1
2       1          1          2
3       1          2          3
4       1          1          4
5       1         11          5
6       1          4          6

虚假数据:

C <- matrix(rnorm(3*79),nrow=3)
T <- matrix(rnorm(3*215),nrow=3)