<font face="trebuchet ms,sans-serif">É isso?<br><br><span style="font-family:courier new,monospace">resp <- c(109.30,124.00,100.00,87.45,115.12,93.10,95.60,101.00,98.50,<br> 68.00,47.00,78.00,118.00,141.60,117.70,112.30,91.30,125.65,115.00,<br>
109.00,97.50,104.00,78.00,87.00,117.70,127.20,116.50,108.00,156.70,<br> 93.10,83.00,111.50,114.00,128.00,65.00,84.00,122.50,114.60,129.00,<br> 108.00,99.75,135.40,111.00,140.00,83.50,89.00,88.50,60.00,126.60,<br>
112.00,123.60,103.20,113.00,158.50,130.00,119.00,87.50,74.00,82.00,67.00)<br>trat <- gl(5,12,label=paste("N",1:5,sep=""))<br>exp <- rep(gl(4,3,label=paste("E",1:4,sep="")),5)<br>
bloco <- rep(gl(3,1,label=paste("bloco",1:3,sep="")),20)<br>dados <- data.frame(trat=trat,exp=exp,bloco=bloco,resp=resp)<br><br>xtabs(~trat+bloco+exp, data=dados)<br><br>mod <- aov(resp ~ trat * exp + exp/bloco, data = dados)<br>
<br># análise individuais<br>mi <- lapply(split(dados, f=dados$exp), aov, formula=resp~bloco+trat)<br>lapply(mi, summary)<br><br># análise conjunta<br>mod <- aov(resp~exp/bloco+trat+Error(exp:trat), data=dados)<br>summary(mod)<br>
</span><br>À disposição.<br>Walmes.<br><br clear="all"></font><div><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Walmes Marques Zeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">fone: (+55) 41 3361 3573</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">VoIP: (3361 3600) 1053 1173</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a><br>
skype: walmeszeviani<br style="font-family:trebuchet ms,sans-serif"></span><span style="font-family:trebuchet ms,sans-serif">twitter: @walmeszeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">linux user number: 531218</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span></div>
<br><br><div class="gmail_quote">2013/6/14 Gilenio Borges Fernandes <span dir="ltr"><<a href="mailto:gilenio@ufba.br" target="_blank">gilenio@ufba.br</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"><div>Gostaria de contar com ajuda para os seguinte problema - analise conjunta</div><div>de 4 experimentos em blocos casualizados completos com 3 repetições</div><div><br></div><div>resp=c(109.30,124.00,100.00,87.45,115.12,93.10,95.60,101.00,98.50,</div>
<div> 68.00,47.00,78.00,118.00,141.60,117.70,112.30,91.30,125.65,115.00,</div><div> 109.00,97.50,104.00,78.00,87.00,117.70,127.20,116.50,108.00,156.70,</div><div> 93.10,83.00,111.50,114.00,128.00,65.00,84.00,122.50,114.60,129.00,</div>
<div> 108.00,99.75,135.40,111.00,140.00,83.50,89.00,88.50,60.00,126.60,</div><div> 112.00,123.60,103.20,113.00,158.50,130.00,119.00,87.50,74.00,82.00,67.00)</div><div>trat=gl(5,12,label=paste("N",1:5,sep="")) </div>
<div>exp<-rep(gl(4,3,label=paste("E",1:4,sep="")),5)</div><div>bloco<-rep(gl(3,1,label=paste("bloco",1:3,sep="")),20)</div><div>dados=data.frame(trat=trat,exp=exp,bloco=bloco,resp=resp)</div>
<div>mod = aov(resp ~ trat * exp + exp/bloco, data = dados)</div><div>summary(mod)</div><div><br></div><div><br></div><div>A declaração do modelo (mod) com a função aov testa tudo com o resíduo</div><div>(32 gl). Quero testar tratamento (trat) e experimento (exp) com a a interação trat*exp e</div>
<div>esta com o resíduo (32 gl). Como declarar os efeitos no modelo?</div><span class="HOEnZb"><font color="#888888"><div><br></div><div><br></div>-- <br>Gilenio Borges Fernandes<div>Professor Associado III<br>Universidade Federal da Bahia<br>
Instituto de Matemática<br>
Departamento de Estatística<br>Av. Adhemar de Barros, s/n – Ondina.<br>40.170-110 - Salvador - BA, Brasil<br>Tel.: (071)3283-6280/6336 Fax: (071)3283-6276<br>URL: <a href="http://lattes.cnpq.br/6764860618464860" target="_blank">http://lattes.cnpq.br/6764860618464860</a><br>
</div>
</font></span></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br>