<html><body><div style="color:#000; background-color:#fff; font-family:verdana, helvetica, sans-serif;font-size:13px"><div dir="ltr" id="yui_3_16_0_1_1423085275275_10573"><span id="yui_3_16_0_1_1423085275275_10706">Cesar, obrigado pela resposta. <br></span></div><div id="yui_3_16_0_1_1423085275275_10707" dir="ltr"><br><span id="yui_3_16_0_1_1423085275275_10706"></span></div><div id="yui_3_16_0_1_1423085275275_10708" dir="ltr"><span style="" class="" id="yui_3_16_0_1_1423085275275_10706">Como eu disse, a SQ e o QM estão corretos ( na parte do desdobramento nc dentro de solo, em destaque, onde tenho interesse). <br></span></div><div id="yui_3_16_0_1_1423085275275_10813" dir="ltr"><span style="" class="" id="yui_3_16_0_1_1423085275275_10706">Porém o F e por consequência o p-valor foram calculados pelo R com o QMres=</span> 21.7, quando na verdade( pela primeira anova=av1) deveria ser 6.07. </div><div id="yui_3_16_0_1_1423085275275_10838" dir="ltr">No excel, ao fazer o desdobramento e calcular o F utilizei os mesmos valores de GL, SQ e QM da av1 (sejam: 118, 715.95 e 6.07). </div><div id="yui_3_16_0_1_1423085275275_10920" dir="ltr">A final, o QMR utilizado nos desdobramentos deverá ser o mesmo da "anova geral=av1", não é isto?</div><div id="yui_3_16_0_1_1423085275275_10700"> </div><div id="yui_3_16_0_1_1423085275275_10698" class="signature"><div id="yui_3_16_0_1_1423085275275_10699" style="text-align:center;">Odirley R. Campos<br>Engenheiro Agrônomo UFV/MG</div><div id="yui_3_16_0_1_1423085275275_10697" style="color:rgb(0, 0, 0);font-size:13px;font-family:arial, helvetica, clean, sans-serif;background-color:transparent;font-style:normal;text-align:center;">Doutorando em Solos e Nutrição de Plantas UFV/MG<br><br></div></div><br> <blockquote id="yui_3_16_0_1_1423085275275_10752" style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; margin-top: 5px; padding-left: 5px;"> <div id="yui_3_16_0_1_1423085275275_10758" style="font-family: verdana, helvetica, sans-serif; font-size: 13px;"> <div id="yui_3_16_0_1_1423085275275_10757" style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;"> <div id="yui_3_16_0_1_1423085275275_10762" dir="ltr"> <hr id="yui_3_16_0_1_1423085275275_10795" size="1"> <font id="yui_3_16_0_1_1423085275275_10763" size="2" face="Arial"> <b><span style="font-weight:bold;">De:</span></b> Cesar Rabak <cesar.rabak@gmail.com><br> <b><span style="font-weight: bold;">Para:</span></b> r-br@listas.c3sl.ufpr.br; Odirley Campos <camposagro@yahoo.com.br> <br> <b><span style="font-weight: bold;">Enviadas:</span></b> Quarta-feira, 4 de Fevereiro de 2015 14:06<br> <b><span style="font-weight: bold;">Assunto:</span></b> Re: [R-br] Anova, fatorial triplo, interacao dupla<br> </font> </div> <div id="yui_3_16_0_1_1423085275275_10756" class="y_msg_container"><br><div id="yiv2890314617"><div id="yui_3_16_0_1_1423085275275_10755"><div id="yui_3_16_0_1_1423085275275_10754" dir="ltr">Odirley,<div id="yui_3_16_0_1_1423085275275_10759"><br clear="none"></div><div id="yui_3_16_0_1_1423085275275_10753">Sua informação sobre as somas quadráticas "não baterem" em relação ao Excel não parecem corretas, posto que os valores são iguais em ambas as tabelas que você postou.<br clear="none"><br clear="none">Em relação aos valores F, não é possível discutir porque você não postou a linha de soma quadrática dos resíduos do Excel (e nem o nº de graus de liberdade) . </div><div id="yui_3_16_0_1_1423085275275_10856"><br clear="none"></div><div id="yui_3_16_0_1_1423085275275_10857"><br clear="none"></div></div><div id="yui_3_16_0_1_1423085275275_11041" class="yiv2890314617gmail_extra"><br clear="none"><div id="yui_3_16_0_1_1423085275275_11061" class="yiv2890314617gmail_quote">2015-02-03 19:07 GMT-02:00 Odirley Campos <span id="yui_3_16_0_1_1423085275275_11060" dir="ltr"><<a id="yui_3_16_0_1_1423085275275_11059" rel="nofollow" shape="rect" ymailto="mailto:camposagro@yahoo.com.br" target="_blank" href="mailto:camposagro@yahoo.com.br">camposagro@yahoo.com.br</a>></span>:<br clear="none"><blockquote id="yui_3_16_0_1_1423085275275_11065" class="yiv2890314617gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div class="qtdSeparateBR"><br><br></div><div class="yiv2890314617yqt0816671171" id="yiv2890314617yqt76197"><div id="yui_3_16_0_1_1423085275275_11064"><div id="yui_3_16_0_1_1423085275275_11063" style="color:#000;background-color:#fff;font-family:verdana, helvetica, sans-serif;font-size:13px;"><div><font># Boa noite! Novamente venho recorrer a lista para tirar duvidas sobre desdobramento na ANOVA.</font></div><div dir="ltr"><font>#
Tenho resultados de materia seca de eucalipto de um experimento em
arranjo fatorial triplo (solo; nc; np) com # UE ditribuidas em DBC. <br clear="none"></font></div><div dir="ltr"><font># Dados estão no: <a rel="nofollow" shape="rect" target="_blank" href="http://www.datafilehost.com/d/32efcd4b">http://www.datafilehost.com/d/32efcd4b</a></font></div><div dir="ltr"><font># Gostaria de corrigir os valores do Fcal no desdobramento das interações duplas solo:nc.</font></div><div id="yui_3_16_0_1_1423085275275_11062" dir="ltr"><font># Na anova geral (avi) o GL e SQ do residuo são obtidos normalmente:</font></div><div dir="ltr"><span style="border-collapse:separate;color:rgb(0,0,0);font-family:'Lucida Console';font-size:13px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:15px;text-indent:0px;text-transform:none;white-space:pre-wrap;word-spacing:0px;background-color:rgb(225,226,229);"></span></div><pre id="yui_3_16_0_1_1423085275275_11066" style="font-family:'Lucida Console';font-size:10pt!important;outline-style:none;outline-color:initial;border-top-style:none;border-right-style:none;border-bottom-style:none;border-left-style:none;border-color:initial;white-space:pre-wrap!important;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;line-height:1.2;"><br clear="none"><span style="background-color:rgb(253,248,105);"><span style="border-collapse:separate;color:rgb(0,0,0);"></span></span></pre><pre id="yui_3_16_0_1_1423085275275_11067" style="font-family:'Lucida Console';font-size:10pt!important;outline-style:none;outline-color:initial;border-top-style:none;border-right-style:none;border-bottom-style:none;border-left-style:none;border-color:initial;white-space:pre-wrap!important;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;line-height:1.2;"><span style="color:blue;"><font color="#cd232c"><span>avi<-aov(mspa~bloco + trat, data=a1)<br clear="none">anova(avi)</span></font><br clear="none"></span></pre><br clear="none">Response: mspa
Df Sum Sq Mean Sq F value Pr(>F)
bloco 2 2 0.993 0.1636 0.8493
trat 59 4963 84.119 13.8640 <2e-16 ***
Residuals 118 716 6.067<pre style="font-family:'Lucida Console';font-size:10pt!important;outline-style:none;outline-color:initial;border-top-style:none;border-right-style:none;border-bottom-style:none;border-left-style:none;border-color:initial;white-space:pre-wrap!important;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;line-height:1.2;"><br clear="none">#refiz no excel e os resultados batem <br clear="none"></pre><div dir="ltr"><br clear="none"></div><div dir="ltr"><br clear="none"></div><div dir="ltr"><font color="#cd232c">av1<-aov(mspa~bloco + solo*nc*np, data=a1) # com mais detalhes<br clear="none"></font></div><div dir="ltr"><font color="#cd232c">anova(av1)</font></div><div dir="ltr"><span style="border-collapse:separate;color:rgb(0,0,0);font-family:'Lucida Console';font-size:13px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:15px;text-indent:0px;text-transform:none;white-space:pre-wrap;word-spacing:0px;background-color:rgb(225,226,229);"></span></div><pre id="yui_3_16_0_1_1423085275275_11068" style="font-family:'Lucida Console';font-size:10pt!important;outline-style:none;outline-color:initial;border-top-style:none;border-right-style:none;border-bottom-style:none;border-left-style:none;border-color:initial;white-space:pre-wrap!important;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;line-height:1.2;"><br clear="none">Response: mspa
Df Sum Sq Mean Sq F value Pr(>F)
bloco 2 1.99 0.99 0.1636 0.84927
solo 5 1987.74 397.55 65.5217 < 2.2e-16 ***
nc 1 0.75 0.75 0.1232 0.72621
np 4 2350.35 587.59 96.8431 < 2.2e-16 ***
solo:nc 5 95.86 19.17 3.1599 0.01031 *
solo:np 20 395.63 19.78 3.2603 3.261e-05 ***
nc:np 4 10.54 2.64 0.4345 0.78348
solo:nc:np 20 122.12 6.11 1.0064 0.46028
Residuals 118 715.95 6.07 </pre><div dir="ltr"><br clear="none"></div><div dir="ltr"><span style="border-collapse:separate;color:rgb(0,0,0);font-family:'Lucida Console';font-size:13px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:15px;text-indent:0px;text-transform:none;white-space:pre-wrap;word-spacing:0px;background-color:rgb(225,226,229);"></span></div><pre id="yui_3_16_0_1_1423085275275_11069" style="font-family:'Lucida Console';font-size:10pt!important;outline-style:none;outline-color:initial;border-top-style:none;border-right-style:none;border-bottom-style:none;border-left-style:none;border-color:initial;white-space:pre-wrap!important;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;line-height:1.2;"><span style="white-space:pre-wrap;color:blue;"></span># Após conferir no excel vi que: df, SS, MS, F e Pr estão corretos para av1<br clear="none"># quando fui para o desdobramento conforme segue, obtive os resultados no R: <br clear="none"><br clear="none"><font color="#cd232c">av2 <- aov(mspa~bloco+solo/nc) #nc dentro solo<br clear="none">anova(av2)<br clear="none"> <br clear="none">names(coef(av2)) <br clear="none">length(names(coef(av2)))</font><br clear="none"> <br clear="none"> # buscando as posiçoes por expressões regulares<br clear="none"> <br clear="none"><font color="#cd232c">grep("ES", names(coef(av2))[9:14]) <br clear="none">grep("EU_0-10", names(coef(av2))[9:14]) <br clear="none">grep("EU_40-60", names(coef(av2))[9:14])<br clear="none">grep("SE_0-10", names(coef(av2))[9:14])<br clear="none">grep("SE_40-60", names(coef(av2))[9:14])<br clear="none">grep("TM", names(coef(av2))[9:14])</font><br clear="none"><br clear="none"><br clear="none"><font color="#cd232c">summary(av2, split=list("solo:nc"=list(<br clear="none"> "ES"= c(1),<br clear="none"> "EU_0-10"=c(2),<br clear="none"> "EU_40-60"=c(3),<br clear="none"> "SE_0-10"=c(4), <br clear="none"> "SE_40-60"= c(5),<br clear="none"> "TM"= c(6)<br clear="none"> )))<br clear="none"></font><span style="border-collapse:separate;color:rgb(0,0,0);font-family:'Lucida Console';font-size:13px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:15px;text-indent:0px;text-transform:none;white-space:pre-wrap;word-spacing:0px;background-color:rgb(225,226,229);"></span></pre><pre style="font-family:'Lucida Console';font-size:10pt!important;outline-style:none;outline-color:initial;border-top-style:none;border-right-style:none;border-bottom-style:none;border-left-style:none;border-color:initial;white-space:pre-wrap!important;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;line-height:1.2;"><span style="border-collapse:separate;color:rgb(0,0,0);font-family:'Lucida Console';font-size:13px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:15px;text-indent:0px;text-transform:none;white-space:pre-wrap;word-spacing:0px;background-color:rgb(225,226,229);"></span></pre><pre style="font-family:'Lucida Console';font-size:10pt!important;outline-style:none;outline-color:initial;border-top-style:none;border-right-style:none;border-bottom-style:none;border-left-style:none;border-color:initial;white-space:pre-wrap!important;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;line-height:1.2;"><br clear="none"># saida da av2 abaixo: <br clear="none"><span style="border-collapse:separate;color:rgb(0,0,0);font-family:'Lucida Console';font-size:13px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:15px;text-indent:0px;text-transform:none;white-space:pre-wrap;word-spacing:0px;background-color:rgb(225,226,229);"></span></pre><pre id="yui_3_16_0_1_1423085275275_11070" style="font-family:'Lucida Console';font-size:10pt!important;outline-style:none;outline-color:initial;border-top-style:none;border-right-style:none;border-bottom-style:none;border-left-style:none;border-color:initial;white-space:pre-wrap!important;margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;line-height:1.2;"> Df Sum Sq Mean Sq F value Pr(>F)
bloco 2 2 1.0 0.046 0.9552
solo 5 1988 397.5 18.359 1.73e-14 ***
solo:nc 6 97 16.1 0.744 0.6153
<span style="background-color:rgb(253,248,105);">solo:nc: ES 1 11 10.6 0.488 0.4859
solo:nc: EU_0-10 1 0 0.5 0.022 0.8817
solo:nc: EU_40-60 1 2 2.1 0.096 0.7570
solo:nc: SE_0-10 1 81 81.0 3.741 0.0548 .
solo:nc: SE_40-60 1 0 0.0 0.000 0.9875
solo:nc: TM 1 2 2.5 0.114 0.7362 </span>
Residuals 166 3595 21.7 </pre><div> </div><div dir="ltr"><font># vi que na parte que me interessa em destaque a SQ e QM estão corretas, mas o Fcal e P-valor estão errados pois o GL, SQ e QM do residuo não estão corretos. Conforme conferi pelo excel: <br clear="none"></font></div><div dir="ltr"><font># saida do excel<br clear="none"></font></div><div dir="ltr"><font>
</font></div><div id="yui_3_16_0_1_1423085275275_11071"><font><span lang="EN-US">FV<span> </span><span> </span>Sum Sq<span> </span>Mean Sq<span>
</span>F value<span> </span><span> </span>Pr(>F)</span></font></div><font>
</font><div><font>solo:nc: ES<span> </span><span> </span>10.6<span> </span>10.6<span> </span><span>
</span>1.740<span> </span><span> </span>0.1897</font></div><font>
</font><div><font>solo:nc: EU_0-10<span> </span>0.5<span> </span>0.5<span> </span><span> </span>0.079<span> </span><span> </span>0.2213</font></div><font>
</font><div><font>solo:nc: EU_40-60<span> </span>2.1<span> </span>2.1<span> </span><span> </span>0.343<span> </span><span> </span>0.4406</font></div><font>
</font><div><font>solo:nc: SE_0-10<span> </span>81.0<span> </span>81.0<span> </span><span> </span>13.347<span> </span><span> </span>0.0004</font></div><font>
</font><div><font>solo:nc: SE_40-60<span> </span>0.0<span> </span>0.0<span> </span><span> </span>0.001<span> </span><span> </span>0.0236</font></div><font>
</font><div><font>solo:nc: TM<span> </span>2.5<span> </span>2.5<span> </span>0.406<span> </span><span> </span>0.4748</font></div>
<div dir="ltr"> </div><div dir="ltr"><font># Existe alguma maneira de informar o df e Sum Sq corretos (que estão na avi ou av1) no procedimento da av2 <br clear="none"></font></div><div id="yui_3_16_0_1_1423085275275_11072" dir="ltr"><font># de modo que </font><font>o F e Pr sejam corrigidos na saida do comando?</font></div><div dir="ltr"><font>#
Ou como eu poderia extrair a coluna do Sum Sq ou Mean Sq na forma de
vector ou data.frame para fazer as correções manualmente no R?</font><br clear="none"></div><div dir="ltr"><br clear="none"></div><div><span></span></div><div> </div><div id="yui_3_16_0_1_1423085275275_11074"><div id="yui_3_16_0_1_1423085275275_11073" style="text-align:center;">Odirley R. Campos<br clear="none">Engenheiro Agrônomo UFV/MG</div><div style="color:rgb(0,0,0);font-size:13px;font-family:arial, helvetica, clean, sans-serif;background-color:transparent;font-style:normal;text-align:center;">Doutorando em Solos e Nutrição de Plantas UFV/MG<br clear="none"><br clear="none"></div></div></div></div></div><br clear="none">_______________________________________________<br clear="none">
R-br mailing list<br clear="none">
<a rel="nofollow" shape="rect" ymailto="mailto:R-br@listas.c3sl.ufpr.br" target="_blank" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br clear="none">
<a rel="nofollow" shape="rect" target="_blank" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br clear="none">
Leia o guia de postagem (<a rel="nofollow" shape="rect" target="_blank" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br clear="none"></blockquote></div><br clear="none"></div></div></div><br><br></div> </div> </div> </blockquote> </div></body></html>