Question

我有13个不同维度的矩阵，我想在成对矩阵相关中使用自定义函数（计算Rv系数）。该函数接受两个参数（matrix1，matrix2）并生成一个标量（基本上是一个多变量r值）。我想在所有可能的矩阵对上运行该函数（总共78个相关），并生成一个13乘13的结果Rv值矩阵，其中包含行和列中13个矩阵的名称。我想通过将matricies放在一个列表中并使用double for循环来遍历列表的元素来尝试这样做，但这似乎非常复杂。我在下面给出了一个带有虚拟数据的示例。有没有人对如何处理这个有任何建议？提前谢谢。

# Rv function  
Rv <- function(M1, M2) {  
    tr <- function(x) sum( diag(x) )   
    psd <- function(x) x %*% t(x)   
    AA <- psd(M1)  
    BB <- psd(M2)  
    num <- tr(AA %*% BB)  
    den <- sqrt( tr(AA %*% AA) * tr(BB %*% BB) )  
    Rv <- num / den  
    list(Rv=Rv, "Rv^2"=Rv^2)  
}  

# data in separate matricies  
matrix1 <- matrix(rnorm(100), 10, 10)  
matrix2 <- matrix(rnorm(100), 10, 10)  
# ... etc. up to matrix 13  

# or, in a list  
matrix1 <- list( matrix(rnorm(100), 10, 10) )  
rep(matrix1, 13) # note, the matrices are identical in this example   

# call Rv function  
Rv1 <- Rv(matrix1, matrix2)  
Rv1$Rv^2  

# loop through all 78 combinations?  
# store results in 13 by 13 matrix with matrix rownames and colnames?

Answer 1

我过去使用的是expand.grid()后跟apply()。这是一个更简单的例子，只使用1：3而不是1:13。

R> work <- expand.grid(1:3,1:3)
R> work
  Var1 Var2
1    1    1
2    2    1
3    3    1
4    1    2
5    2    2
6    3    2
7    1    3
8    2    3
9    3    3
R> apply(work, 1, function(z) prod(z))
[1] 1 2 3 2 4 6 3 6 9
R>

你显然想要一个不同的工人职能。

R中的成对矩阵相关 - 如何迭代所有对？

1 个答案: