Mauricio,<div><br></div><div>Que bom que minha suspeita e sugestão tenham sido proveitosas...</div><div><br></div><div>Quanto ao "erro" fico sem ter como ajudar sem maiores detalhes do que, ao seu ver, esteja "errado"! Você está vizualizando tal erro por meio de quê? Quantos níveis existem em cada fator? Qual a resposta do seu str(couto)? Existe desbalanceamento nos seus dados? Você chegou a ler a documentação da função terms('formula', keep = TRUE)? Você pretende fazer os testes F todos pelo termo do Residuals?</div>
<div><br></div><div>att,</div><div>FH<br><br><div class="gmail_quote">2012/6/25 Mauricio Couto <span dir="ltr"><<a href="mailto:couto_estatistica@hotmail.com" target="_blank">couto_estatistica@hotmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div dir="ltr">
Toledo,<br>segui suas sugestões e em parte deu certo, contudo as intereções dentro de outro fator, tipo AMB*Familias/SET ou REP/AMB*SET, tá errado como vc pode ver abaixo:<br><br>couto<-read.table("C:/Users/Mauricio/Desktop/SR.txt",h=T)<br>
> couto<-transform(couto,AMB=factor(AMB),SET=factor(SET),Famílias=factor(Famílias),REP=factor(REP))<br>> couto<-aov(CE~AMB*SET+REP/AMB*SET+Famílias/SET+ AMB*Famílias/SET,data=couto)<br>> summary(couto)<br><br>
<font face="Courier New"> Df Sum Sq Mean Sq F value Pr(>F) </font><font face="Courier New"><br></font><font face="Courier New">AMB 1 10 10.46 1.287 0.25726 </font><font face="Courier New"><br>
</font><font face="Courier New">SET 7 403 57.52 7.077 5.85e-08 ***</font><font face="Courier New"><br></font><font face="Courier New">REP 1 57 56.56 6.959 0.00868 ** </font><font face="Courier New"><br>
</font><font face="Courier New">Famílias 192 13551 70.58 8.684 < 2e-16 ***</font><font face="Courier New"><br></font><font face="Courier New">AMB:SET 7 98 13.97 1.719 0.10298 </font><font face="Courier New"><br>
</font><font face="Courier New">AMB:REP 1 8 7.84 0.965 0.32661 </font><font face="Courier New"><br></font><font face="Courier New">SET:REP 7 53 7.58 0.933 0.48060 </font><font face="Courier New"><br>
</font><font face="Courier New">AMB:Famílias 192 2012 10.48 1.289 0.01911 * </font><font face="Courier New"><br></font><font face="Courier New">AMB:SET:REP 7 43 6.11 0.751 0.62845 </font><font face="Courier New"><br>
</font><font face="Courier New">Residuals <a href="tel:384%C2%A0%C2%A0%203121%C2%A0%C2%A0%C2%A0%208.13" value="+553843121813" target="_blank">384 3121 8.13</a> </font><font face="Courier New"><br>
</font><font face="Courier New">---</font><font face="Courier New"><br></font><font face="Courier New">Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 </font><br><br><br>Ps: desculpe o incomodo.<div class="im">
<br><br><p align="center"><font color="#008000"><strong>Mauricio Farias Couto</strong></font></p>
<p align="center"><strong><font color="#008000">Doutorando em Genética e Melhoramento de Plantas (UENF-RJ)</font></strong></p>
<p align="center"><strong><font color="#008000">Mestre em Estatística Aplicada e Biometria (UFV-MG)</font></strong></p>
<p align="center"><font color="#008000"><strong>Engº Agronomo (UESC-BA)</strong></font></p>
<p align="center"><font color="#008000"><strong></strong></font> </p>
<p align="center"><font color="#008000"><strong>Cel: <a href="tel:%2822%29%208117-3725" value="+552281173725" target="_blank">(22) 8117-3725</a></strong></font></p>
<p align="center"><strong><font color="#008000"> </font></strong></p><br><br></div><div><div></div><hr>Date: Wed, 20 Jun 2012 15:55:10 -0300<br>From: <a href="mailto:fernandohtoledo@gmail.com" target="_blank">fernandohtoledo@gmail.com</a><br>
To: <a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a><br>Subject: Re: [R-br] Ajuda na análise de seleção recorrente<div><div class="h5"><br><br><div>Maurício,</div><div><br></div><div>
Vendo o que você postou... e consultando meu oráculo, posso suspeitar que o "erro" na sua ANOVA seja que ele está interpretando suas variáveis como numérica e não categóricas!</div>
<div><br></div><div>Note que isso é um erro de semântica e não sintaxe!</div><div><br></div><div>Para ver se é esse o caso, reporte-nos a saída do seus </div><div><br></div><div>> str(couto)</div><div><br></div><div>A afirmativa que esses efeitos estejam sendo interpretados como variáveis numéricas, proceda:</div>
<div><br></div><div>> couto <- transform(couto,</div><div>> AMB = factor(AMB),</div><div>> ... outras variáveis a transformar ... )</div><div><br></div><div>
Então proceda de novo sua anova pelo aov(), fica a sugestão de usar o terms('modelo', keep = TRUE) no argumento formula da função... se for o caso de obter somas de quadrados ajustadas!</div><div><br></div><div>OBS: Em outra ocasião que postar, leia antes o guia de postagem e produza um CMR, isso facilita para quem se disponibiliza a ajudar!</div>
<div><br></div><div>sem mais...</div><div><br></div><div>à disposição (parafrazeando Walmes),</div><div>FH</div><div><br></div><div><br></div>Faltou você transformar suas variáveis em fator<br><br><div>
2012/6/20 Mauricio Couto <span dir="ltr"><<a href="mailto:couto_estatistica@hotmail.com" target="_blank">couto_estatistica@hotmail.com</a>></span><br><blockquote style="border-left:1px #ccc solid;padding-left:1ex">
<div><div dir="ltr">
Pessoal.<br>
<br>
gostaria da ajuda de vocês para analisar um modelo de seleção recorrente no R, o problema são os efeitos aninhados com interação (R/AS). Necessito estimar os componentes de variância para estimar por exemplo a herdabilidade entre outros.<br>
<br>
Modelo estatístico:<br>
<p><span style="font-family:Arial"><font size="3">Y<sub>ijkl</sub> = µ + A<sub>i</sub> + S<sub>j</sub> + AS<sub>ij</sub>+R/AS<sub>ijk</sub> +F/S<sub>jl</sub> + AS/S<sub>ijl</sub> + ε <sub>ijkl</sub><u></u><u></u><u></u></font></span></p>
<p><span style="line-height:115%;font-family:'Times New Roman';font-size:12pt">μ é a média<u></u><u></u></span></p>
<p><span style="line-height:115%;font-family:'Times New Roman';font-size:12pt">A<sub>i</sub> é o efeito do ambiente<u></u><u></u></span></p>
<p><span style="line-height:115%;font-family:'Times New Roman';font-size:12pt">S<sub>j</sub> é o efeito de set<u></u><u></u></span></p>
<p><span style="line-height:115%;font-family:'Times New Roman';font-size:12pt">AS<sub>ij</sub> é o efeito da interação ambiente e set<u></u><u></u></span></p>
<p><span style="line-height:115%;font-family:'Times New Roman';font-size:12pt">R/AS<sub>ijk</sub> é o efeito da repetição dentro de ambiente e set<u></u><u></u></span></p>
<p><span style="line-height:115%;font-family:'Times New Roman';font-size:12pt">F/S<sub>jl</sub> é o efeito da família dentro de set<u></u><u></u></span></p>
<p><span style="line-height:115%;font-family:'Times New Roman';font-size:12pt">AF/S<sub>ijl</sub> é o efeito da interação ambiente e família dentro de set<u></u><u></u></span></p>
<p><span style="line-height:115%;font-family:'MS Mincho';font-size:12pt">ɛ</span><sub><span style="line-height:115%;font-family:'Times New Roman';font-size:12pt">ijkl</span></sub><span style="line-height:115%;font-family:'Times New Roman';font-size:12pt"> é o erro experimental</span></p>
<p><span style="line-height:115%;font-family:'Times New Roman';font-size:12pt"></span> </p>
<p><span style="line-height:115%;font-family:'Times New Roman';font-size:12pt">Obs? Vejam como eu tentei fazer:</span></p><span style="line-height:115%;font-family:'Times New Roman';font-size:12pt">
<p><span style="font-family:'Times New Roman'" lang="EN-US">> couto<-read.table("C:/Users/Mauricio/Desktop/SR.txt", h=T)</span></p><span style="font-family:'Times New Roman'" lang="EN-US">
<p><span style="font-family:'Times New Roman'" lang="EN-US">>couto<-aov(CE~AMB*SET+REP/AMB*SET+Famílias/SET+AMB*Famílias/SET,data=couto)<u></u><u></u></span></p>
<u></u><u></u></span></span>
<p><span style="font-family:'Times New Roman'" lang="EN-US"><font size="3">> summary(couto)<u></u><u></u></font></span></p>
Obs: só que a ANOVA saiu errada.<br>
<br>
Desde de já agradeço.<br>
<br>
<br>
<br><br> <br>
<p align="center"><font color="#008000"><strong>Mauricio Farias Couto</strong></font></p>
<p align="center"><strong><font color="#008000">Doutorando em Genética e Melhoramento de Plantas (UENF-RJ)</font></strong></p>
<p align="center"><strong><font color="#008000">Mestre em Estatística Aplicada e Biometria (UFV-MG)</font></strong></p>
<p align="center"><font color="#008000"><strong>Engº Agronomo (UESC-BA)</strong></font></p>
<p align="center"><font color="#008000"><strong></strong></font> </p>
<p align="center"><font color="#008000"><strong>Cel: <a>(22) 8117-3725</a></strong></font></p>
<p align="center"><strong><font color="#008000"> </font></strong></p>
<br> </div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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></div></div>_______________________________________________
R-br mailing list
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<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>
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.</div> </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></div>