R中的3D数组列表的平均值

Question

如何计算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) 

[...] 

应该产生包含平均值的维数4,5,5的数组。

我可以做它矩阵列表：

apply(simplify2array(myList), 1:2, mean) 

，但它并不适用于我的目的工作...

感谢您的任何线索！

Answer 1

可能是我们可以使用Reduce

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

Answer 2

你可以先改变你数组列表变为四维阵列，然后使用apply：

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