数组 – R对于每个第三维,通过向量乘以3D数组的第二维

当尝试将数组的第一个维度乘以向量的每个索引乘以第二个维度时,我的数组将转换为矩阵,并且事情变得松散.我只能长时间做适当的乘法.

多么满口……

使用代码更容易解​​释……

Arr <- array(runif(10*5*3), dim = c(10,5,3))
dim(Arr)
Vect <- c(1:5)

Arr[,1,1] <- Arr[,1,1]*Vect[1]
Arr[,1,2] <- Arr[,1,2]*Vect[1]
Arr[,1,3] <- Arr[,1,3]*Vect[1]

Arr[,2,1] <- Arr[,2,1]*Vect[2]
Arr[,2,2] <- Arr[,2,2]*Vect[2]
Arr[,2,3] <- Arr[,2,3]*Vect[2]

Arr[,3,1] <- Arr[,3,1]*Vect[3]
Arr[,3,2] <- Arr[,3,2]*Vect[3]
Arr[,3,3] <- Arr[,3,3]*Vect[3]

Arr[,4,1] <- Arr[,4,1]*Vect[4]
Arr[,4,2] <- Arr[,4,2]*Vect[4]
Arr[,4,3] <- Arr[,4,3]*Vect[4]

Arr[,5,1] <- Arr[,5,1]*Vect[5]
Arr[,5,2] <- Arr[,5,2]*Vect[5]
Arr[,5,3] <- Arr[,5,3]*Vect[5]

如何将其清理为一个命令?

最佳答案
首先将Vect转换为数组,然后将元素乘以:

varr <- aperm(array(Vect, dim = c(5L, 10L, 3L)), perm = c(2L, 1L, 3L))

Arr <- varr * Arr

(当然,如果你想在一个命令中使用它,我们不需要存储varr)

(另外,事实证明这基本上是扫描引擎盖下的……)

转载注明原文:数组 – R对于每个第三维,通过向量乘以3D数组的第二维 - 代码日志