在R中使用对多维数组的应用

Question

我想知道如何在多维数组上使用apply.我有以下内容:

A <- array(0, c(2, 2, 5))
for(i in 1:5) {
  A[, , i] <- matrix(rnorm(4), 2, 2)
}

我想取这些切片的平均值来得到一个2乘2的矩阵.我提出的任何方式都非常糟糕.

我希望能够使用apply,就像我想要矩阵的列的平均说法一样:

B <- matrix(rnorm(10), 5, 2)
B.mean <- apply(B, 2, mean)

但这似乎与我认为3D阵列的方式不同:

A.mean <- apply(A, 3, mean)

我感谢你的建议.

Answer 1

A.mean <- apply(A, c(1,2), mean)