相同维度的两个矩阵之间的Spearman相关性

时间:2015-03-16 03:34:58

标签: r loops matrix correlation psych

我有两个相同维度的矩阵(p和e),我想在同名的列之间建立一个spearman相关。我想在矩阵(M)

中得到对相关的输出

我使用了图书馆Psych的corr.test()功能,这就是我所做的:

library(psych)
M <- data.frame(matrix(ncol=3,nrow=ncol(p)))
M[,1] <- as.character()
G <- colnames(p)
for(rs in 1:ncol(p){
      M[rs,1] <- G[rs]     
      cor <- corr.test(p[,rs],e[,rs],method="spearman",adjust="none")
      M[rs,2] <- cor$r
      M[rs,3] <- cor$p
}

但是我收到一条错误消息:

Error in 1:ncol(y) : argument of length 0
你能告诉我出什么问题吗?或建议另一种方法?

2 个答案:

答案 0 :(得分:4)

不需要所有这些循环和索引等:

# test data
p <- matrix(data = rnorm(100),nrow = 10)
e <- matrix(data = rnorm(100),nrow = 10)

cor <- corr.test(p, e, method="spearman", adjust="none")
data.frame(name=colnames(p), r=diag(cor$r), p=diag(cor$p))

#  name           r         p
#a    a  0.36969697 0.2930501
#b    b  0.16363636 0.6514773
#c    c -0.15151515 0.6760652
# etc etc

如果矩阵的名称不匹配,则match

cor <- corr.test(p, e[,match(colnames(p),colnames(e))], method="spearman", adjust="none")

答案 1 :(得分:0)

由于两个矩阵很大,所以在所有可能的对上执行函数corr.test()需要很长的system.time,但最终工作的循环如下:

    library(psych)
    M <- data.frame(matrix(ncol=3,nrow=ncol(p)))
    M[,1] <- as.character()
    G <- colnames(p)
    for(rs in 1:ncol(p){
          M[rs,1] <- G[rs]     
          cor <- corr.test(as.data.frame(p[,rs]),as.data.frame(e[,rs]),
method="spearman",adjust="none")
          M[rs,2] <- cor$r
          M[rs,3] <- cor$p
    }