Opções adicionais consiste em usar funções de pacotes para manipulação de dados, veja da <- expand.grid(A=1:14, B=1:4, r=1:4) da <- cbind(da, as.data.frame(matrix(rnorm(nrow(da)*3), nrow(da)))) str(da) require(reshape) require(plyr) ddply(da[,-c(1,2)], .(da$A, da$B), mean) ddply(da[,-c(1,2)], .(da$A, da$B), sd) db <- melt(da, id=1:3) str(db) ddply(db[,-3], .(A, B, variable), summarise, m=mean(value), sd=sd(value)) library(doBy) summaryBy(V1+V2+V3~A+B, data=da, FUN=function(x){ c(mean(x), sd(x)) }) À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================