Média e desvio de muitas variáveis de um fatorial!
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! \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}
um exemplo reproduzivel... set.seed(1) da = data.frame(Trat=sample(6, 1000, rep=T), x=rnorm(1000), y=rexp(1000, 1)) aggregate(da[,-1], by=list(Trat=da[['Trat']]), function(x) c(media=mean(x, na.rm=T), stderr=sd(x, na.rm=T))) b
Perfeito! Obrigado! \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature} ________________________________ De: Benilton Carvalho <beniltoncarvalho@gmail.com> Para: r-br@listas.c3sl.ufpr.br; Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> Enviadas: Terça-feira, 18 de Outubro de 2011 8:26 Assunto: Re: [R-br] Média e desvio de muitas variáveis de um fatorial! um exemplo reproduzivel... set.seed(1) da = data.frame(Trat=sample(6, 1000, rep=T), x=rnorm(1000), y=rexp(1000, 1)) aggregate(da[,-1], by=list(Trat=da[['Trat']]), function(x) c(media=mean(x, na.rm=T), stderr=sd(x, na.rm=T))) b
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 ==========================================================================
participantes (3)
-
Benilton Carvalho -
Ivan Bezerra Allaman -
Walmes Zeviani