Проблемы хранения корреляции в матрице с помощью rcorr из пакета Hmisc

Question

Я хочу вычислить p-значения для корреляции двух матриц, используя Hmisc. У меня проблема с сохранением результирующих значений 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])
    }
}

Возвращает эту ошибку:

Error in P[i, j] <- rcorr(x[, i], y[, j]) : 
  number of items to replace is not a multiple of replacement length