使用Hmisc包中的rcorr在矩阵中存储相关性时出现的问题

我想使用

Hmisc

计算两个矩阵相关性的p值。我在矩阵中存储生成的p值时遇到问题：

library("Hmisc")
x <- as.matrix(read.table("tmp01"))
y <- as.matrix(read.table("tmp02"))
P <- matrix(0, ncol(x), ncol(y))

for (i in 1:ncol(x)) {
    for(j in 1:ncol(y)) {
        P[i,j] <- rcorr(x[,i],y[,j])
    }
}

库（“Hmisc”）
x
rcorr
返回一个列表，因此您需要仅对p值进行子集运算，例如，使用类似于
rcorr（x[，i]，y[，j]）$p[1,2]
非常感谢您的帮助！
Error in P[i, j] <- rcorr(x[, i], y[, j]) : 
  number of items to replace is not a multiple of replacement length