将每一行乘以每一行

Question

我有一个20k行和100列的矩阵。我想在所有可能的行对之间执行除法运算。这将创建＆lt; 20k选择2＆gt;矢量组合。如何在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], ]
  }
)