R中三维阵列列表的平均值

如何计算R中列表中存储的数组的平均值？结果应该是一个大小相同的数组，其中包含。。。例如：

a1 <- array(runif(100), dim=c(4, 5, 5))
a2 <- array(runif(100), dim=c(4, 5, 5))
a3 <- array(runif(100), dim=c(4, 5, 5))
a4 <- array(runif(100), dim=c(4, 5, 5))
a5 <- array(runif(100), dim=c(4, 5, 5))

l <- list(a1, a2, a3, a4, a5)

[...]

但这不符合我的目的

---

谢谢你的提示

也许我们可以使用

Reduce

Reduce(`+`, l)/length(l)

---

您可以先将数组列表转换为4D数组，然后使用

apply

：

a <- array(unlist(l), dim = c(dim(l[[1]]), length(l)))
apply(a, 1:3, mean)

---

a感谢您的选择！
a <- array(unlist(l), dim = c(dim(l[[1]]), length(l)))
apply(a, 1:3, mean)