一个矩阵(n×m)与另一矩阵(n×n)的所有列相乘的总和

时间:2018-11-13 10:13:05

标签: r matrix linear-algebra

我希望标题是不言自明的。我希望制作这样的标题,以便其他人也可以找到它。我知道如何使用循环执行该操作,但是使用某种矩阵乘法必须更快,这是我感兴趣的。

带有循环的代码看起来像这样

x <- matrix(rexp(300, rate=.1), nrow=20)
y <- matrix(rexp(400, rate=.1), nrow=20)

res <- as.data.frame(matrix(0,ncol = 15, nrow = 20))

for (i in 1:20){

  res <- res + x*y[,i]

}

0 个答案:

没有答案