Caro Fernado,
O loop é muito simples. Segue um script em que analisei 11 variáveis de uma só vez e ainda apliquei dois testes de comparações de medias. As duas primeiras colunas sao respectivamente
Trat e Repetiçoes, as demais são as variavies.
Att
--
José Wellingthon dos Santos
Pesquisador - Métodos Quantitativos
Embrapa Algodão
Empresa Brasileira de Pesquisa Agropecuária (Embrapa)
Campina Grande, PB

jose-wellingthon.santos@embrapa.br
Telefone: +55 (83) 3182-4401 | Fax: +55 (83) 3182-4367
www.embrapa.br/algodao | https://twitter.com/embrapa_algodao
Confira também: www.facebook.com/agrosustentavel



De: "Fernando Rodrigo Bortolozo via R-br" <r-br@listas.c3sl.ufpr.br>
Para: r-br@listas.c3sl.ufpr.br
Enviadas: Terça-feira, 6 de dezembro de 2016 16:29:24
Assunto: [R-br] Loop script

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.  

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.