r - 乘以不同维度的数组

Question

我想乘(或除)两个不同维度的数组(a1=[48,38,31] 和 b1[48,38])。到目前为止，我在第三维上使用了 for 循环。但是，我想知道如何使用(如果可能的话)apply。假设我有以下示例:

a1<- array(rnorm(20), dim=c(2,3,3))
b1<- array(rnorm(20), dim=c(2,3))

如果我尝试直接执行 a1/b1(或 *)我不能，因为它们需要具有相同的尺寸。所以我使用了一个 for 循环:

for(i in 1:3){
  m1[,,i] <- a1[,,i]/b1
}

但我想避免使用循环。

Answer 1

一个选项是

array(c(a1)/rep(b1, dim(a1)[3]), dim= dim(a1))

或者我们可以使用apply

apply(a1, 3, FUN=function(x) x/b1)