Fernando, > long_shape <- reshape2::melt(dados, id.vars = 'Trat') > model <- function(df) lm(value ~ Trat, data = df) > fits <- plyr::ldply(long_shape, .(variable), model)
Deste modo, "fits" é uma lista que armazena os ajustes para todas as variáveis. Minha sugestão e usar outro **ply para tratar tais ajustes e obter os testes de média e demais analíses necessárias. Att, FH 2016-12-06 13:29 GMT-06: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. >
_______________________________________________ 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.