Arrays 数组r的矩阵加法

我有一个包含272个矩阵的数组，每个矩阵是2乘2。现在我想用矩阵加法把这些矩阵求和。所以我希望返回是一个2乘2的矩阵。下面是我使用过的一些代码

 y <- as.matrix(faithful) 
 B <- matrix(c(0,0,0,0),nrow = 2)


 sigma <- function(n = 272,u_new) {
     vec = replicate(272,B)

     for (i in 1:n) {

w <- (y-u_new)[i,]
x <- ptilde1[i]*(w%*%t(w))
vec[,,i][1,1] <- x[1,1]
vec[,,i][1,2] <- x[1,2]
vec[,,i][2,1] <- x[2,1]
vec[,,i][2,2] <- x[2,2]}

  vec
  }

---

y这里是循环多次（272次）并向同一列表添加矩阵的代码

B <- matrix(c(0,0,0,0),nrow = 2)
list <- list(B)
for (i in 2:272) {
    list[[i]] <- B
}

这是一个人为的例子，因为所有的矩阵都是零矩阵。我将把它作为家庭作业留给你，让你使用你真正想要求和的矩阵。
可能重复@TimBiegeleisen，我试过了。它没用，只给了我一个号码。我猜那是因为我的向量是一个数组。不是列表，我不知道如何建立矩阵列表。
sum <- Reduce('+', list)
> sum
     [,1] [,2]
[1,]    0    0
[2,]    0    0