<html><body><div style="color:#000; background-color:#fff; font-family:verdana, helvetica, sans-serif;font-size:13px"><div style="" class="" id="yui_3_16_0_1_1422993873832_4128"><font style="" class="" id="yui_3_16_0_1_1422993873832_6738" size="2"># Boa noite! Novamente venho recorrer a lista para tirar duvidas sobre desdobramento na ANOVA.</font></div><div style="" class="" id="yui_3_16_0_1_1422993873832_4129" dir="ltr"><font style="" class="" id="yui_3_16_0_1_1422993873832_6737" size="2">#
Tenho resultados de materia seca de eucalipto de um experimento em
arranjo fatorial triplo (solo; nc; np) com # UE ditribuidas em DBC. <br style="" class=""></font></div><div style="" class="" id="yui_3_16_0_1_1422993873832_4134" dir="ltr"><font style="" class="" id="yui_3_16_0_1_1422993873832_6736" size="2"># Dados estão no: <a style="" class="" id="yui_3_16_0_1_1422993873832_4227" href="http://www.datafilehost.com/d/32efcd4b">http://www.datafilehost.com/d/32efcd4b</a></font></div><div style="" class="" id="yui_3_16_0_1_1422993873832_4228" dir="ltr"><font style="" class="" id="yui_3_16_0_1_1422993873832_6735" size="2"># Gostaria de corrigir os valores do Fcal no desdobramento das interações duplas solo:nc.</font></div><div style="" class="" id="yui_3_16_0_1_1422993873832_4253" dir="ltr"><font style="" class="" id="yui_3_16_0_1_1422993873832_6734" size="2"># Na anova geral (avi) o GL e SQ do residuo são obtidos normalmente:</font></div><div style="" class="" id="yui_3_16_0_1_1422993873832_4256" dir="ltr"><span id="yui_3_16_0_1_1422993873832_4280" class="" 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; orphans: 2; text-align: -webkit-left; text-indent: 0px; text-transform: none; white-space: pre-wrap; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; background-color: rgb(225, 226, 229); "></span></div><pre id="yui_3_16_0_1_1422993873832_4279" tabindex="0" class="" style="font-family: 'Lucida Console'; font-size: 10pt !important; outline-style: none; outline-width: initial; outline-color: initial; border-top-style: none; border-right-style: none; border-bottom-style: none; border-left-style: none; border-width: initial; border-color: initial; white-space: pre-wrap !important; word-break: break-all; margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; -webkit-user-select: text; line-height: 1.2; "><br style="" class=""><span class="" id="yui_3_16_0_1_1422993873832_6367" style="background-color: rgb(253, 248, 105);"><span id="yui_3_16_0_1_1422993873832_4314" class="" 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; orphans: 2; text-indent: 0px; text-transform: none; white-space: pre-wrap; widows: 2; word-spacing: 0px;"><pre id="yui_3_16_0_1_1422993873832_4313" tabindex="0" class="" style="font-family: 'Lucida Console'; font-size: 10pt !important; outline-style: none; outline-width: initial; outline-color: initial; border-top-style: none; border-right-style: none; border-bottom-style: none; border-left-style: none; border-width: initial; border-color: initial; white-space: pre-wrap !important; word-break: break-all; margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; -webkit-user-select: text; line-height: 1.2; "><span id="yui_3_16_0_1_1422993873832_4312" class="" style="color: blue; "><font style="" class="" id="yui_3_16_0_1_1422993873832_6819" color="#cd232c"><span style="" class="" id="yui_3_16_0_1_1422993873832_6242">avi<-aov(mspa~bloco + trat, data=a1)<br style="" class="">anova(avi)</span></font><br style="" class=""></span></pre></span></span><br style="" class="">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><pre id="yui_3_16_0_1_1422993873832_4279" tabindex="0" class="" style="font-family: 'Lucida Console'; font-size: 10pt !important; outline-style: none; outline-width: initial; outline-color: initial; border-top-style: none; border-right-style: none; border-bottom-style: none; border-left-style: none; border-width: initial; border-color: initial; white-space: pre-wrap !important; word-break: break-all; margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; -webkit-user-select: text; line-height: 1.2; "><br style="" class="">#refiz no excel e os resultados batem <br style="" class=""></pre><div style="" class="" id="yui_3_16_0_1_1422993873832_4254" dir="ltr"><br style="" class=""></div><div style="" class="" id="yui_3_16_0_1_1422993873832_4290" dir="ltr"><br style="" class=""></div><div style="" class="" id="yui_3_16_0_1_1422993873832_4338" dir="ltr"><font style="" class="" id="yui_3_16_0_1_1422993873832_6067" size="2" color="#cd232c">av1<-aov(mspa~bloco + solo*nc*np, data=a1) # com mais detalhes<br style="" class=""></font></div><div style="" class="" id="yui_3_16_0_1_1422993873832_4339" dir="ltr"><font style="" class="" id="yui_3_16_0_1_1422993873832_6191" size="2" color="#cd232c">anova(av1)</font></div><div style="" class="" id="yui_3_16_0_1_1422993873832_4341" dir="ltr"><span id="yui_3_16_0_1_1422993873832_4373" class="" 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; orphans: 2; text-align: -webkit-left; text-indent: 0px; text-transform: none; white-space: pre-wrap; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; background-color: rgb(225, 226, 229); "></span></div><pre id="yui_3_16_0_1_1422993873832_4372" tabindex="0" class="" style="font-family: 'Lucida Console'; font-size: 10pt !important; outline-style: none; outline-width: initial; outline-color: initial; border-top-style: none; border-right-style: none; border-bottom-style: none; border-left-style: none; border-width: initial; border-color: initial; white-space: pre-wrap !important; word-break: break-all; margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; -webkit-user-select: text; line-height: 1.2; "><br style="" class="">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 style="" class="" id="yui_3_16_0_1_1422993873832_4289" dir="ltr"><br style="" class=""></div><div style="" class="" id="yui_3_16_0_1_1422993873832_4135" dir="ltr"><span id="yui_3_16_0_1_1422993873832_4168" class="" 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; orphans: 2; text-align: -webkit-left; text-indent: 0px; text-transform: none; white-space: pre-wrap; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; background-color: rgb(225, 226, 229); "></span></div><pre id="yui_3_16_0_1_1422993873832_4167" tabindex="0" class="" style="font-family: 'Lucida Console'; font-size: 10pt !important; outline-style: none; outline-width: initial; outline-color: initial; border-top-style: none; border-right-style: none; border-bottom-style: none; border-left-style: none; border-width: initial; border-color: initial; white-space: pre-wrap !important; word-break: break-all; margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; -webkit-user-select: text; line-height: 1.2; "><span class="" style="white-space: pre; -webkit-user-select: text; color: blue; "></span># Após conferir no excel vi que: df, SS, MS, F e Pr estão corretos para av1<br style="" class=""># quando fui para o desdobramento conforme segue, obtive os resultados no R: <br style="" class=""><br style="" class=""><font id="yui_3_16_0_1_1422993873832_8842" color="#cd232c">av2 <- aov(mspa~bloco+solo/nc) #nc dentro solo<br style="" class="">anova(av2)<br style="" class=""> <br style="" class="">names(coef(av2)) <br style="" class="">length(names(coef(av2)))</font><br style="" class=""> <br style="" class=""> # buscando as posiçoes por expressões regulares<br style="" class=""> <br style="" class=""><font id="yui_3_16_0_1_1422993873832_8872" color="#cd232c">grep("ES", names(coef(av2))[9:14]) <br style="" class="">grep("EU_0-10", names(coef(av2))[9:14]) <br style="" class="">grep("EU_40-60", names(coef(av2))[9:14])<br style="" class="">grep("SE_0-10", names(coef(av2))[9:14])<br style="" class="">grep("SE_40-60", names(coef(av2))[9:14])<br style="" class="">grep("TM", names(coef(av2))[9:14])</font><br style="" class=""><br style="" class=""><br style="" class=""><font color="#cd232c">summary(av2, split=list("solo:nc"=list(<br style="" class=""> "ES"= c(1),<br style="" class=""> "EU_0-10"=c(2),<br style="" class=""> "EU_40-60"=c(3),<br style="" class=""> "SE_0-10"=c(4), <br style="" class=""> "SE_40-60"= c(5),<br style="" class=""> "TM"= c(6)<br style="" class=""> )))<br style="" class=""></font><span id="yui_3_16_0_1_1422993873832_4168" class="" 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; orphans: 2; text-align: -webkit-left; text-indent: 0px; text-transform: none; white-space: pre-wrap; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; background-color: rgb(225, 226, 229); "><pre id="yui_3_16_0_1_1422993873832_4167" tabindex="0" class="" style="font-family: 'Lucida Console'; font-size: 10pt !important; outline-style: none; outline-width: initial; outline-color: initial; border-top-style: none; border-right-style: none; border-bottom-style: none; border-left-style: none; border-width: initial; border-color: initial; white-space: pre-wrap !important; word-break: break-all; margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; -webkit-user-select: text; line-height: 1.2; "><span id="yui_3_16_0_1_1422993873832_4168" class="" 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; orphans: 2; text-align: -webkit-left; text-indent: 0px; text-transform: none; white-space: pre-wrap; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; background-color: rgb(225, 226, 229); "><pre id="yui_3_16_0_1_1422993873832_4167" tabindex="0" class="" style="font-family: 'Lucida Console'; font-size: 10pt !important; outline-style: none; outline-width: initial; outline-color: initial; border-top-style: none; border-right-style: none; border-bottom-style: none; border-left-style: none; border-width: initial; border-color: initial; white-space: pre-wrap !important; word-break: break-all; margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; -webkit-user-select: text; line-height: 1.2; "><br style="" class=""># saida da av2 abaixo: <br style="" class=""><span id="yui_3_16_0_1_1422993873832_4617" class="" 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; orphans: 2; text-align: -webkit-left; text-indent: 0px; text-transform: none; white-space: pre-wrap; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; background-color: rgb(225, 226, 229); "><pre id="yui_3_16_0_1_1422993873832_4616" tabindex="0" class="" style="font-family: 'Lucida Console'; font-size: 10pt !important; outline-style: none; outline-width: initial; outline-color: initial; border-top-style: none; border-right-style: none; border-bottom-style: none; border-left-style: none; border-width: initial; border-color: initial; white-space: pre-wrap !important; word-break: break-all; margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; -webkit-user-select: text; 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 class="" id="yui_3_16_0_1_1422993873832_4778" 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></span></pre></span></pre></span></pre><div style="" class="" id="yui_3_16_0_1_1422993873832_6401"> </div><div style="" class="" id="yui_3_16_0_1_1422993873832_6402" dir="ltr"><font style="" class="" id="yui_3_16_0_1_1422993873832_6638" size="2"># 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 style="" class=""></font></div><div style="" class="" id="yui_3_16_0_1_1422993873832_6495" dir="ltr"><font style="" class="" id="yui_3_16_0_1_1422993873832_6642" size="2"># saida do excel<br style="" class=""></font></div><div style="" class="" id="yui_3_16_0_1_1422993873832_6482" dir="ltr"><font style="" class="" size="2">
</font></div><div id="yui_3_16_0_1_1422993873832_6583" style="" class=""><font style="" class="" id="yui_3_16_0_1_1422993873832_6636" size="2"><span id="yui_3_16_0_1_1422993873832_6582" class="" style="mso-ansi-language:EN-US" lang="EN-US">FV<span id="yui_3_16_0_1_1422993873832_6641" class="" style="mso-tab-count:1"> </span><span id="yui_3_16_0_1_1422993873832_6581" class="" style="mso-spacerun:yes"> </span>Sum Sq<span id="yui_3_16_0_1_1422993873832_6612" class="" style="mso-tab-count:1"> </span>Mean Sq<span id="yui_3_16_0_1_1422993873832_6720" class="" style="mso-spacerun:yes">
</span>F value<span id="yui_3_16_0_1_1422993873832_6610" class="" style="mso-tab-count:1"> </span><span id="yui_3_16_0_1_1422993873832_6635" class="" style="mso-spacerun:yes"> </span>Pr(>F)</span></font></div><font style="" class="" size="2">
</font><div style="" id="yui_3_16_0_1_1422993873832_6529" class=""><font style="" class="" id="yui_3_16_0_1_1422993873832_6634" size="2">solo:nc: ES<span id="yui_3_16_0_1_1422993873832_6640" class="" style="mso-tab-count:1"> </span><span class="" id="yui_3_16_0_1_1422993873832_6528" style="mso-spacerun:yes"> </span>10.6<span id="yui_3_16_0_1_1422993873832_6598" class="" style="mso-tab-count:1"> </span>10.6<span id="yui_3_16_0_1_1422993873832_6611" class="" style="mso-tab-count:1"> </span><span id="yui_3_16_0_1_1422993873832_6679" class="" style="mso-spacerun:yes">
</span>1.740<span id="yui_3_16_0_1_1422993873832_6721" class="" style="mso-tab-count:1"> </span><span id="yui_3_16_0_1_1422993873832_6609" class="" style="mso-spacerun:yes"> </span>0.1897</font></div><font style="" class="" size="2">
</font><div id="yui_3_16_0_1_1422993873832_6585" style="" class=""><font style="" class="" id="yui_3_16_0_1_1422993873832_6625" size="2">solo:nc: EU_0-10<span id="yui_3_16_0_1_1422993873832_6599" class="" style="mso-tab-count:1"> </span>0.5<span id="yui_3_16_0_1_1422993873832_6597" class="" style="mso-tab-count:1"> </span>0.5<span id="yui_3_16_0_1_1422993873832_6584" class="" style="mso-tab-count:1"> </span><span id="yui_3_16_0_1_1422993873832_6595" class="" style="mso-spacerun:yes"> </span>0.079<span id="yui_3_16_0_1_1422993873832_6733" class="" style="mso-tab-count:1"> </span><span id="yui_3_16_0_1_1422993873832_6608" class="" style="mso-spacerun:yes"> </span>0.2213</font></div><font style="" class="" size="2">
</font><div id="yui_3_16_0_1_1422993873832_6594" style="" class=""><font style="" class="" id="yui_3_16_0_1_1422993873832_6624" size="2">solo:nc: EU_40-60<span id="yui_3_16_0_1_1422993873832_6600" class="" style="mso-tab-count:1"> </span>2.1<span id="yui_3_16_0_1_1422993873832_6648" class="" style="mso-tab-count:1"> </span>2.1<span id="yui_3_16_0_1_1422993873832_6596" class="" style="mso-tab-count:1"> </span><span id="yui_3_16_0_1_1422993873832_6593" class="" style="mso-spacerun:yes"> </span>0.343<span id="yui_3_16_0_1_1422993873832_6623" class="" style="mso-tab-count:1"> </span><span id="yui_3_16_0_1_1422993873832_6607" class="" style="mso-spacerun:yes"> </span>0.4406</font></div><font style="" class="" size="2">
</font><div id="yui_3_16_0_1_1422993873832_6592" style="" class=""><font style="" class="" id="yui_3_16_0_1_1422993873832_6622" size="2">solo:nc: SE_0-10<span id="yui_3_16_0_1_1422993873832_6601" class="" style="mso-tab-count:1"> </span>81.0<span id="yui_3_16_0_1_1422993873832_6647" class="" style="mso-tab-count:1"> </span>81.0<span id="yui_3_16_0_1_1422993873832_6646" class="" style="mso-tab-count:1"> </span><span id="yui_3_16_0_1_1422993873832_6591" class="" style="mso-spacerun:yes"> </span>13.347<span id="yui_3_16_0_1_1422993873832_6621" class="" style="mso-tab-count:1"> </span><span id="yui_3_16_0_1_1422993873832_6606" class="" style="mso-spacerun:yes"> </span>0.0004</font></div><font style="" class="" size="2">
</font><div id="yui_3_16_0_1_1422993873832_6589" style="" class=""><font style="" class="" id="yui_3_16_0_1_1422993873832_6629" size="2">solo:nc: SE_40-60<span id="yui_3_16_0_1_1422993873832_6602" class="" style="mso-tab-count:1"> </span>0.0<span class="" style="mso-tab-count:1"> </span>0.0<span id="yui_3_16_0_1_1422993873832_6588" class="" style="mso-tab-count:1"> </span><span id="yui_3_16_0_1_1422993873832_6590" class="" style="mso-spacerun:yes"> </span>0.001<span id="yui_3_16_0_1_1422993873832_6604" class="" style="mso-tab-count:1"> </span><span id="yui_3_16_0_1_1422993873832_6605" class="" style="mso-spacerun:yes"> </span>0.0236</font></div><font style="" class="" size="2">
</font><div id="yui_3_16_0_1_1422993873832_6586" style="" class=""><font style="" class="" id="yui_3_16_0_1_1422993873832_6631" size="2">solo:nc: TM<span id="yui_3_16_0_1_1422993873832_6644" class="" style="mso-tab-count:1"> </span>2.5<span id="yui_3_16_0_1_1422993873832_6603" class="" style="mso-tab-count:1"> </span>2.5<span id="yui_3_16_0_1_1422993873832_6632" class="" style="mso-tab-count:1"> </span>0.406<span id="yui_3_16_0_1_1422993873832_6587" class="" style="mso-tab-count:1"> </span><span id="yui_3_16_0_1_1422993873832_6633" class="" style="mso-spacerun:yes"> </span>0.4748</font></div>
<div style="" class="" id="yui_3_16_0_1_1422993873832_6482" dir="ltr"> </div><div id="yui_3_16_0_1_1422993873832_6691" style="" class="" dir="ltr"><font style="" class="" id="yui_3_16_0_1_1422993873832_6699" size="2"># Existe alguma maneira de informar o df e Sum Sq corretos (que estão na avi ou av1) no procedimento da av2 <br style="" class=""></font></div><div id="yui_3_16_0_1_1422993873832_6704" style="" class="" dir="ltr"><font style="" class="" id="yui_3_16_0_1_1422993873832_6699" size="2"># de modo que </font><font style="" class="" id="yui_3_16_0_1_1422993873832_6700" size="2">o F e Pr sejam corrigidos na saida do comando?</font></div><div id="yui_3_16_0_1_1422993873832_6873" style="" class="" dir="ltr"><font style="" class="" id="yui_3_16_0_1_1422993873832_6880" size="2">#
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 style="" class=""></div><div id="yui_3_16_0_1_1422993873832_6874" style="" class="" dir="ltr"><br style="" class=""></div><div id="yui_3_16_0_1_1422993873832_8868"><span></span></div><div id="yui_3_16_0_1_1422993873832_8869"> </div><div id="yui_3_16_0_1_1422993873832_8512" class="signature"><div id="yui_3_16_0_1_1422993873832_8681" style="text-align:center;">Odirley R. Campos<br>Engenheiro Agrônomo UFV/MG</div><div id="yui_3_16_0_1_1422993873832_8511" 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></div></body></html>