在R中的矩阵的列上应用成对函数

Question

我想在大矩阵的所有列上应用一些成对相关估计。由于我正在处理大尺寸，我正在寻找一个可能对此有帮助的功能。我一直在试验apply函数，但我还没有走得太远所以非常感谢所有的帮助。

require(ccaPP)
require(mvtnorm)
d<- 10
Sigma <- matrix(0.2, nrow = d, ncol = d)
diag(Sigma) <- 1
#Data generation
X <- rmvnorm(100, sigma = Sigma) # 100 x d matrix
Q <- apply(X, 2, FUN = corQuadrant, consistent = TRUE)

如果之前已经问过这个问题，我很抱歉，但搜索没有透露我可以使用的内容。谢谢。

Answer 1

嵌套for循环是一种方式：

require(ccaPP)
require(mvtnorm)
d<- 10
Sigma <- matrix(0.2, nrow = d, ncol = d)
diag(Sigma) <- 1
#Data generation
X <- rmvnorm(100, sigma = Sigma) # 100 x d matrix
Q <- matrix(nrow = d,ncol = d)
for (i in 1: (d-1)) { 
  for (j in (i+1):d) {
    Q[i,j] <- corQuadrant( X[,i], X[,j] ) 
  }
}