Uma alternativa

colMeans(cbind(m1, m2, m3), na.rm=TRUE)

b