Bom dia senhores!

Tenho um experimento fatorial 14x4, com 12 variáveis, e gostaria de saber se tem como colocar de uma só vez em um único quadro a média e o desvio padrão dos tratamentos e de todas as variáveis. Fiz algo como:

medias  <- apply(da[,-c(1:2,15)],2,function(x) tapply(x,da$Trat,mean,na.rm=TRUE))

desvios <- apply(da[,-c(1:2,15)],2,function(x) tapply(x,da$Trat,sd,na.rm=TRUE)) 

descrit <- list()

for(i in 1:12){

  descrit[[i]] <- cbind(medias[,i],desvios[,i])

}                                                  

Mais sai como lista não sendo o ideal que procuro. 

Desde já agradeço pela atenção!