R:计算两个向量的 "product"的最佳方法

Question

假设我有一个向量

r <- rnorm(4)

和一个矩阵 W尺寸为 20000*200 例如:

W <- matrix(rnorm(20000*200),20000,200)

我想计算一个新矩阵 M尺寸为 5000*200，使得 m11 <- r%*%W[1:4,1] , m21 <- r%*%W[5:8,1] , m12 <- r%*%W[1:4,2]等(即按 4×4 对行进行分组并计算乘积)。

这样做的最佳(速度，内存)方式是什么？

提前致谢。

Answer 1

这对我来说似乎运行得最快:

array(r %*% array(W, c(4, 20000 * 200 / 4)), c(5000, 200))