r - R中矩阵列与向量元素相乘的最快方法

Question

我有一个矩阵m和一个向量v。我想将矩阵 m 的第一列乘以向量 v 的第一个元素，并将矩阵 m 的第二列乘以第二个元素向量v的元素，依此类推。我可以使用以下代码来完成此操作，但我正在寻找一种不需要两个转置调用的方法。如何在 R 中更快地完成此操作？

m <- matrix(rnorm(120000), ncol=6)
v <- c(1.5, 3.5, 4.5, 5.5, 6.5, 7.5)

system.time(t(t(m) * v))

#   user  system elapsed 
#   0.02    0.00    0.02 

Answer 1

使用一些线性代数并执行矩阵乘法，这在 R 中相当快。

例如

m %*% diag(v)

一些基准测试

m = matrix(rnorm(1200000), ncol=6)

 v=c(1.5, 3.5, 4.5, 5.5, 6.5, 7.5)
library(microbenchmark)
microbenchmark(m %*% diag(v), t(t(m) * v))
##    Unit: milliseconds
##            expr      min       lq   median       uq      max neval
##   m %*% diag(v) 16.57174 16.78104 16.86427 23.13121 109.9006   100
##     t(t(m) * v) 26.21470 26.59049 32.40829 35.38097 122.9351   100