engo problemas con las operaciones de vectores o matrices.

Tengo tres matrices 3D y quiero encontrar el promedio de ellas. ¿Cómo puedo hacer eso? no podemos usarmean() ya que solo devuelve un único valor.

Lo más importante es que algunas de las celdas en las matrices son NA, lo que significa que si solo las agrego como

A = (B + C + D)/3 

Los resultados de mostrarán NA también.

¿Cómo puedo dejar que reconozca si la celda es NA y luego omitirla?

Me gust

 A = c(NA, 10, 15, 15, NA)
 B = c(10, 15, NA, 22, NA)
 C = c(NA, NA, 20, 26, NA)

Quiero que la salida del promedio de estos vectores sea

(10, (10+15)/2, (15+20)/2, (15+22+26)/3, NA)

Tampoco podemos usarna.omit, porque moverá el orden de los índices.

Este es el código correspondiente. desearía que fuera útil.

for (yr in 1950:2011) {
    temp_JFM <- sst5_sst2[,,year5_sst2==yr & (month5_sst2>=1 & month5_sst2<=3)]
       k = 0
       jfm=4*k+1
    for (i in 1:72) {
        for (j in 1:36) {
            iposst5_sst2[i,j,jfm] <- (temp_JFM[i,j,1]+temp_JFM[i,j,2]+temp_JFM[i,j,3])/3
        }
    }      
}

Gracias

Ya se resolvió.

a forma más fácil de corregirlo se puede mostrar a continuació

iposst5_sst2[i,j,jfm] <- mean(temp_JFM[i,j,],na.rm=TRUE)