<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'>
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 style="" face="Courier New"> Df Sum Sq Mean Sq F value Pr(>F) </font><font style="" face="Courier New"><br></font><font style="" face="Courier New">AMB 1 10 10.46 1.287 0.25726 </font><font style="" face="Courier New"><br></font><font style="" face="Courier New">SET 7 403 57.52 7.077 5.85e-08 ***</font><font style="" face="Courier New"><br></font><font style="" face="Courier New">REP 1 57 56.56 6.959 0.00868 ** </font><font style="" face="Courier New"><br></font><font style="" face="Courier New">Famílias 192 13551 70.58 8.684 < 2e-16 ***</font><font style="" face="Courier New"><br></font><font style="" face="Courier New">AMB:SET 7 98 13.97 1.719 0.10298 </font><font style="" face="Courier New"><br></font><font style="" face="Courier New">AMB:REP 1 8 7.84 0.965 0.32661 </font><font style="" face="Courier New"><br></font><font style="" face="Courier New">SET:REP 7 53 7.58 0.933 0.48060 </font><font style="" face="Courier New"><br></font><font style="" face="Courier New">AMB:Famílias 192 2012 10.48 1.289 0.01911 * </font><font style="" face="Courier New"><br></font><font style="" face="Courier New">AMB:SET:REP 7 43 6.11 0.751 0.62845 </font><font style="" face="Courier New"><br></font><font style="" face="Courier New">Residuals 384 3121 8.13 </font><font style="" face="Courier New"><br></font><font style="" face="Courier New">---</font><font style="" face="Courier New"><br></font><font style="" face="Courier New">Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 </font><br><br><br>Ps: desculpe o incomodo.<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: (22) 8117-3725</strong></font></p>
<p align="center"><strong><font color="#008000"> </font></strong></p><br><br><div><div id="SkyDrivePlaceholder"></div><hr id="stopSpelling">Date: Wed, 20 Jun 2012 15:55:10 -0300<br>From: fernandohtoledo@gmail.com<br>To: r-br@listas.c3sl.ufpr.br<br>Subject: Re: [R-br] Ajuda na análise de seleção recorrente<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 class="ecxgmail_quote">
2012/6/20 Mauricio Couto <span dir="ltr"><<a href="mailto:couto_estatistica@hotmail.com">couto_estatistica@hotmail.com</a>></span><br><blockquote class="ecxgmail_quote" 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 style="" class="ecxMsoNormal"><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 style="" class="ecxMsoNormal"><span style="line-height:115%;font-family:'Times New Roman';font-size:12pt">μ é a média<u></u><u></u></span></p>
<p style="" class="ecxMsoNormal"><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 style="" class="ecxMsoNormal"><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 style="" class="ecxMsoNormal"><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 style="" class="ecxMsoNormal"><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 style="" class="ecxMsoNormal"><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 style="" class="ecxMsoNormal"><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 style="" class="ecxMsoNormal"><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 style="" class="ecxMsoNormal"><span style="line-height:115%;font-family:'Times New Roman';font-size:12pt"></span> </p>
<p style="" class="ecxMsoNormal"><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 style="" class="ecxMsoNormal"><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 style="" class="ecxMsoNormal"><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 style="" class="ecxMsoNormal"><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 target="_blank">(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">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></body>
</html>