Fernando, Você poderia pensar em fazer uma manova e em seguida uma aov do objeto que armazena a manova. Por exemplo, no código fonte da documentação da manova você pode fazer uma aov do objeto npk2.aov npk2 <- within(npk, foo <- rnorm(24)) (npk2.aov <- manova(cbind(yield, foo) ~ block + N*P*K, npk2)) summary(aov(npk2.aov)) HTH, ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)99966.9100 - VIVO 55(73)99100.7351 - TIM 55(73)98817.6159 - OI 55(73)98129.9942 - CLARO ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ 2016-12-06 16:29 GMT-03:00 Fernando Rodrigo Bortolozo via R-br < r-br@listas.c3sl.ufpr.br>:
Eu estou fazendo anova e com posterior teste de médias de um experimento com muitas variáveis. Então eu tenho que repetir o script muitas vezes o mesmo script para cada variável.
Os colegas sabem me dizer se existe alguma maneira de fazer um Loop para um script, ou mesmo criar um elemento com todas variáveis para analisar de uma só vez?
Aqui vai uma pequena parte das variáveis que vou analisar, tenho 7 tratamentos, 4 repetições.
[image: pasted1]
Aqui esta o script (Agricolae Package)
dados <- read.csv("~/Documents/R/Data306090DAPv2.csv", header = TRUE) head(dados) model<-aov(Cas30 ~Trat, data=dados) out <- HSD.test(model,"Trat", group=TRUE,console=TRUE, main="Calcio") par(mfrow=c(1,2)) out$means par (mar = c (3,3,2,0),cex=0.9) bar.group(out$groups,ylim=c(0, 360),density=50,border="blue", col="blue", las=1) bar.err(out$means,variation="SE",horiz=FALSE, ylim = c(0, 360), col=colors()[15],space=0.3,bar=TRUE, las=1) out<-HSD.test(model,"Trat", group=FALSE) means<-out$means df<-df.residual(model) MSerror<-deviance(model)/df with(dados,HSD.test(Cas30,Trat,df,MSerror, group=TRUE,console=TRUE, main= "Cálcio"))
Muito obrigado desde já.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
