数组 – 从n,k维矩阵的矩阵数中减去n,k维的矩阵

如果我有阵列A.

A <- array(0, c(4, 3, 5))
for(i in 1:5) {
  set.seed(i)
  A[, , i] <- matrix(rnorm(12), 4, 3)
}

如果我有矩阵B

set.seed(6)
B <- matrix(rnorm(12), 4, 3)

从数组A的每个矩阵中减去B的代码是:

d<-array(0, c(4,3,5))
for(i in 1:5){
  d[,,i]<-A[,,i]-B
}

但是,使用“apply”系列函数执行相同计算的代码是什么?

最佳答案
这就是扫描的目的.

sweep(A, 1:2, B)

转载注明原文:数组 – 从n,k维矩阵的矩阵数中减去n,k维的矩阵 - 代码日志