[R-br] Loop script

[Jose Claudio Faria]

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 em 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 em 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161206/4a1b151d/attachment.html>
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: pasted1
Tipo: image/png
Tamanho: 56704 bytes
Descrição: não disponível
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161206/4a1b151d/attachment.png>