我有两个相同维度的矩阵(p 和 e),我想在同名的列之间建立一个 spearman 相关性。我想在矩阵(M)中输出对相关性
我使用了corr.test()
Psych 库中的函数,这就是我所做的:
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
你能告诉我有什么问题吗?或建议另一种方法?