<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>
Só completando, o pacote ExpDes é a maneira mais fácil, mas caso você queira fazer da forma que esta fazendo é só utilizar os seguintes comandos:<div><br></div><div><font face="Courier New">rend<-read.table("C:\\Program Files\\R\\RDADOS\\UFPR1.txt",h=T)</font></div><div><span style="font-family: 'Courier New'; font-size: 10pt;">attach(rend)</span></div><div><div><font face="Courier New">A<-factor(A)</font></div><div><font face="Courier New">K<-factor(K)</font></div></div><div><font face="Courier New">bloc<-factor(bloc)</font></div><div><span style="font-family: 'Courier New'; font-size: 10pt;">m1 <- aov(rg~bloc+A/K)</span></div><div><span lang="EN-US" style="font-size:10.0pt;line-height:
115%;font-family:"Courier New";mso-fareast-font-family:"Times New Roman";
mso-ansi-language:EN-US;mso-fareast-language:PT-BR;mso-bidi-language:AR-SA">KinA<-sapply(paste("A",levels(A),sep=""),simplify=F,grep,x=names(coef(m1)[m1$assign==3]))</span></div><div><span lang="EN-US" style="font-size:10.0pt;line-height:
115%;font-family:"Courier New";mso-fareast-font-family:"Times New Roman";
mso-ansi-language:EN-US;mso-fareast-language:PT-BR;mso-bidi-language:AR-SA"><p class="MsoNormal" style="margin-bottom:0cm;margin-bottom:.0001pt;text-align:
justify;line-height:150%;tab-stops:45.8pt 91.6pt 137.4pt 183.2pt 229.0pt 274.8pt 320.6pt 366.4pt 412.2pt 458.0pt 503.8pt 549.6pt 595.4pt 641.2pt 687.0pt 732.8pt"><span lang="EN-US" style="font-size: 10pt; line-height: 150%;">DesKinA<-summary(m1, split=list("A:K"=KinA))<o:p></o:p></span></p><p class="MsoNormal" style="margin-bottom:0cm;margin-bottom:.0001pt;text-align:
justify;line-height:150%;tab-stops:45.8pt 91.6pt 137.4pt 183.2pt 229.0pt 274.8pt 320.6pt 366.4pt 412.2pt 458.0pt 503.8pt 549.6pt 595.4pt 641.2pt 687.0pt 732.8pt"><span lang="EN-US" style="font-size: 10pt; line-height: 150%;">DesKinA</span></p><p class="MsoNormal" style="margin-bottom:0cm;margin-bottom:.0001pt;text-align:
justify;line-height:150%;tab-stops:45.8pt 91.6pt 137.4pt 183.2pt 229.0pt 274.8pt 320.6pt 366.4pt 412.2pt 458.0pt 503.8pt 549.6pt 595.4pt 641.2pt 687.0pt 732.8pt"><span lang="EN-US" style="font-size: 10pt; line-height: 150%;"><br></span></p><p class="MsoNormal" style="margin-bottom:0cm;margin-bottom:.0001pt;text-align:
justify;line-height:150%;tab-stops:45.8pt 91.6pt 137.4pt 183.2pt 229.0pt 274.8pt 320.6pt 366.4pt 412.2pt 458.0pt 503.8pt 549.6pt 595.4pt 641.2pt 687.0pt 732.8pt"><span lang="EN-US" style="font-size: 10pt; line-height: 150%;"><div style="font-family: Tahoma; line-height: normal; text-align: start;"><font face="Courier New">m2 <- aov(rg~bloc+K/A)</font></div><div style="font-family: Tahoma; line-height: normal; text-align: start;"><span lang="EN-US" style="font-size: 10pt; line-height: 14px; font-family: 'Courier New';">AinK<-sapply(paste("K",levels(K),sep=""),simplify=F,grep,x=names(coef(m2)[m2$assign==3]))</span></div><div style="font-family: Tahoma; line-height: normal; text-align: start;"><span lang="EN-US" style="font-size: 10pt; line-height: 14px; font-family: 'Courier New';"><p class="MsoNormal" style="margin-bottom: 0.0001pt; text-align: justify; line-height: 19px;"><span lang="EN-US" style="font-size: 10pt;">DesAinK<-summary(m2, split=list("A:K"=AinK))<o:p></o:p></span></p><p class="MsoNormal" style="margin-bottom: 0.0001pt; text-align: justify; line-height: 19px;"><span lang="EN-US" style="font-size: 10pt;">DesAinK</span></p></span></div></span></p><p class="MsoNormal" style="margin-bottom:0cm;margin-bottom:.0001pt;text-align:
justify;line-height:150%;tab-stops:45.8pt 91.6pt 137.4pt 183.2pt 229.0pt 274.8pt 320.6pt 366.4pt 412.2pt 458.0pt 503.8pt 549.6pt 595.4pt 641.2pt 687.0pt 732.8pt"><span lang="EN-US" style="font-size: 10pt; line-height: 150%;"><br></span></p><p class="MsoNormal" style="margin-bottom:0cm;margin-bottom:.0001pt;text-align:
justify;line-height:150%;tab-stops:45.8pt 91.6pt 137.4pt 183.2pt 229.0pt 274.8pt 320.6pt 366.4pt 412.2pt 458.0pt 503.8pt 549.6pt 595.4pt 641.2pt 687.0pt 732.8pt"><span lang="EN-US" style="font-size: 10pt; line-height: 150%;">Outra forma seria pelo método dos efeitos aninhados</span></p><p class="MsoNormal" style="margin-bottom:0cm;margin-bottom:.0001pt;text-align:
justify;line-height:150%;tab-stops:45.8pt 91.6pt 137.4pt 183.2pt 229.0pt 274.8pt 320.6pt 366.4pt 412.2pt 458.0pt 503.8pt 549.6pt 595.4pt 641.2pt 687.0pt 732.8pt"><span lang="EN-US" style="font-size: 10pt; line-height: 150%;"><br></span></p><p class="MsoNormal" style="margin-bottom:0cm;margin-bottom:.0001pt;text-align:
justify;line-height:150%;tab-stops:45.8pt 91.6pt 137.4pt 183.2pt 229.0pt 274.8pt 320.6pt 366.4pt 412.2pt 458.0pt 503.8pt 549.6pt 595.4pt 641.2pt 687.0pt 732.8pt"><span lang="EN-US" style="font-size: 10pt; line-height: 150%;"><div style="font-family: Tahoma; line-height: normal; text-align: start;"><font face="Courier New">rend<-read.table("C:\\Program Files\\R\\RDADOS\\UFPR1.txt",h=T)</font></div><div style="font-family: Tahoma; line-height: normal; text-align: start;"><span style="font-size: 10pt;"><font face="Courier New">attach(rend)</font></span></div><div style="font-family: Tahoma; line-height: normal; text-align: start;"><div><font face="Courier New">A<-factor(A)</font></div><div><font face="Courier New">K<-factor(K)</font></div><div><span style="font-family: 'Courier New';">bloc<-factor(bloc)</span></div></div><div style="font-family: Tahoma; line-height: normal; text-align: start;"><font face="Courier New">m1 <- aov(rg~bloc+A/K)</font></div><div style="font-family: Tahoma; line-height: normal; text-align: start;"><font face="Courier New">effects(m1) #A partir da matriz de efeitos você seleciona os efeitos de interesse.</font></div><div style="font-family: Tahoma; line-height: normal; text-align: start;"><span style="font-family: 'Courier New'; line-height: 19px; text-align: justify;">DesKinA<-summary(m1, split=list("A:K"=list(A37.5=c(1,4,7,10),A50=c(2,5,8,11),A62.5=c(3,6,9,12))))</span></div><div style="font-family: Tahoma; line-height: normal; text-align: start;"><span style="font-family: 'Courier New'; line-height: 19px; text-align: justify;">DesKinA</span></div><div style="font-family: Tahoma; line-height: normal; text-align: start;"><span style="font-family: 'Courier New';">m2 <- aov(rg~bloc+K/A)</span></div><div style="font-family: Tahoma; line-height: normal; text-align: start;"><span style="font-family: 'Courier New';">effects(m2)</span></div><div style="font-family: Tahoma; line-height: normal; text-align: start;"><span style="font-family: 'Courier New'; line-height: 19px; text-align: justify;">DesAinK<-summary(m2, split=list("K:A"=list(K0=c(1,6),K30=c(2,7),K60=c(3,8),K120=c(4,9),K180=c(5,10))))</span></div><div style="font-family: Tahoma; line-height: normal; text-align: start;"><span style="font-family: 'Courier New'; line-height: 19px; text-align: justify;">DesAinK</span></div><div style="font-family: Tahoma; line-height: normal; text-align: start;"><span style="font-family: 'Courier New'; line-height: 19px; text-align: justify;"><br></span></div><div style="font-family: Tahoma; line-height: normal; text-align: start;">Espero que tenhamos sanado sua duvida.</div><div style="font-family: Tahoma; line-height: normal; text-align: start;"><br></div><div style="font-family: Tahoma; line-height: normal; text-align: start;">Att.</div><div style="font-family: Tahoma; line-height: normal; text-align: start;"><br></div><div style="font-family: Tahoma; line-height: normal; text-align: start;">Tiago.</div></span></p></span></div><div><br><div><div id="SkyDrivePlaceholder"></div><hr id="stopSpelling">Date: Sat, 17 Nov 2012 09:21:23 -0200<br>From: jhguilhen@gmail.com<br>To: r-br@listas.c3sl.ufpr.br; camposagro@yahoo.com.br<br>Subject: Re: [R-br] Experimento fatorial. Desdobramento<br><br>Uma outra maneira de fazer isto e bem mais fácil é utilizar o pacote ExpDes, ele já faz o desdobramento automaticamente quando dá interação e bem fácil de utilizar.<br><br>Atte..<br><br><br><br><div class="ecxgmail_quote">Em 16 de novembro de 2012 21:41, Odirley Campos <span dir="ltr"><<a href="mailto:camposagro@yahoo.com.br">camposagro@yahoo.com.br</a>></span> escreveu:<br>
<blockquote class="ecxgmail_quote" style="border-left:1px #ccc solid;padding-left:1ex"><div><div style="font-size:10pt;font-family:tahoma,new york,times,serif"><div><span>#Pessoal! Boa noite ! <br></span></div>
<div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:tahoma,new york,times,serif">#Estou tentando sem sucesso fazer o desdobramento de interações com R. Cheguei a baixar alguns arquivos da net no endereço "<span style="text-decoration:underline"><a href="http://www.leg.ufpr.br/doku.php/ridiculas#desdobramento_de_interacao_em_experimento_fatorial" target="_blank">http://www.leg.ufpr.br/doku.php/ridiculas#desdobramento_de_interacao_em_experimento_fatorial</a></span>" . Porêm tive o mesmo problema que tive com meus dados. <br>
</div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:tahoma,new york,times,serif">#Usei os mesmos comandos retirados da net para os dados tirados de lá. A anova "geral" saiu normalmente:</div>
<div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:tahoma,new york,times,serif"><br></div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:tahoma,new york,times,serif">
>rend<-read.table("C:\\Program Files\\R\\RDADOS\\UFPR1.txt",h=T) #após armezenar arquivo no meu PC#</div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:tahoma,new york,times,serif">
>rend<br></div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:tahoma,new york,times,serif">>m1 <- aov(rg~bloc+A/K, data=rend)</div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:tahoma,new york,times,serif">
>summary(m1) #Anova geral saiu norma sem problemas#<br></div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:tahoma,new york,times,serif"> Df Sum Sq Mean Sq F value Pr(>F) <br>
bloc <a target="_blank">1 51.3 51.3 3.585</a> 0.0624 . <br>A 1 495.6 495.6 34.633 1.2e-07 ***<br>
A:K 1 2213.7 2213.7 154.712 < 2e-16 ***<br>Residuals 71 1015.9 14.3 <br>---<br>Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 <br></div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:tahoma,new york,times,serif">
<br></div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:tahoma,new york,times,serif">>summary(m1,split=list("A:K"=list("A-37.5"=c(1,4,7,10),"A-50.0"=c(2,5,8,11),"A-62.5"=c(3,6,9,12)))) #<span style="background-color:rgb(255,255,0)"><span>Mas o desdobramento não mostra Quadrado médio nem mesmo, F-value ou mesmo o P-value. O que eu posso ter feito de errado??</span></span><br>
</div><div style="font-style:normal;font-size:13.3333px;background-color:transparent;font-family:tahoma,new york,times,serif"> Df Sum Sq Mean Sq F value Pr(>F) <br>bloc 1 51.3
51.3 3.585 0.0624 . <br>A 1 495.6 495.6 34.633 1.2e-07 ***<br>A:K 1 2213.7 2213.7 154.712 < 2e-16 ***<br> A:K: A-37.5 4 <br> A:K: A-50.0 4 <br>
A:K: A-62.5
4 <br>Residuals 71 1015.9 14.3 <br>---<br>Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 </div><span class="ecxHOEnZb"><font color="#888888"><div style="text-align:center;font-style:normal;font-size:13px;background-color:transparent;font-family:arial,helvetica,clean,sans-serif">
Odirley Campos<br>Engenheiro Agrônomo UFV/MG<br></div></font></span></div></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>
<br>_______________________________________________
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.</div></div> </div></body>
</html>