My next suggestion (I don't have time to work out or test an example at the moment):
library(abind) tmparr <- abind(m1,m2,m3,...,along=3) OR tmparr <- do.call(c(matlist,list(along=3))) apply(tmparr,c(1,2),mean,na.rm=TRUE) or something along those lines. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
