每行乘以每行

Question

我有一个20k行和100列的矩阵。我想在所有可能的行对之间执行除法操作。这会创建向量的组合。如何在R编程语言中最有效地做到这一点？

Answer 1

n <- 20 
x <- matrix(rnorm(2 * n), nrow = n) 
apply(
    combn(nrow(x), 2), 
    2, 
    function(y){ 
    x[y[1], ]/x[y[2], ] 
    } 
)