<font face="trebuchet ms,sans-serif">A ordem com que você declara não muda a precisão dos efeitos estimados. Ajuste com inversão de ordem e veja o erro padrão summary() do modelo. É o mesmo.<br><br><span style="font-family:courier new,monospace">> d1 <- read.table("<a href="http://www.leg.ufpr.br/~walmes/data/montgomery_14-1.txt">http://www.leg.ufpr.br/~walmes/data/montgomery_14-1.txt</a>",<br>
+ header=TRUE, sep="\t", colClasses=c("factor","factor","integer"))<br>> str(d1)<br>'data.frame': 18 obs. of 3 variables:<br> $ vidro : Factor w/ 2 levels "1","2": 1 1 1 2 2 2 1 1 1 2 ...<br>
$ fosforo : Factor w/ 3 levels "800","825","850": 1 1 1 1 1 1 2 2 2 2 ...<br> $ corrente: int 280 290 285 230 235 240 300 310 295 260 ...<br>> d1$corrente[1] <- NA<br>> m0 <- lm(corrente~vidro*fosforo, d1)<br>
> m1 <- lm(corrente~fosforo*vidro, d1)<br>> anova(m0)<br> Df Sum Sq Mean Sq F value Pr(>F) <br>vidro 1 14309.0 14309.0 264.1663 4.884e-09 ***<br>fosforo 2 <b>869.5</b> 434.7 8.0261 0.007088 ** <br>
vidro:fosforo 2 81.5 40.8 0.7527 0.493861 <br>Residuals 11 595.8 54.2 <br>> anova(m1)<br> Df Sum Sq Mean Sq F value Pr(>F) <br>fosforo 2 <b>1119.2</b> 559.6 10.3312 0.002983 ** <br>
vidro 1 14059.3 14059.3 259.5560 5.361e-09 ***<br>fosforo:vidro 2 81.5 40.8 0.7527 0.493861 <br>Residuals 11 595.8 54.2 <br><br>> summary(m0)<br>Coefficients:<br> Estimate Std. Error t value Pr(>|t|) <br>
(Intercept) 287.5000 5.2042 55.244 8.43e-15 ***<br>vidro2 -52.5000 <b>6.7185</b> -7.814 8.16e-06 ***<br>fosforo825 14.1667 6.7185 2.109 0.0587 . <br>fosforo850 0.8333 6.7185 0.124 0.9035 <br>
vidro2:fosforo825 -4.1667 9.0139 -0.462 0.6529 <br>vidro2:fosforo850 -10.8333 9.0139 -1.202 0.2547 <br><br>> summary(m1)<br>Coefficients:<br> Estimate Std. Error t value Pr(>|t|) <br>
(Intercept) 287.5000 5.2042 55.244 8.43e-15 ***<br>fosforo825 14.1667 6.7185 2.109 0.0587 . <br>fosforo850 0.8333 6.7185 0.124 0.9035 <br>vidro2 -52.5000 <b>6.7185</b> -7.814 8.16e-06 ***<br>
fosforo825:vidro2 -4.1667 9.0139 -0.462 0.6529 <br>fosforo850:vidro2 -10.8333 9.0139 -1.202 0.2547 <br></span><br>À disposição.<br>Walmes.<br><br></font><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Walmes Marques Zeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">fone: (+55) 41 3361 3573</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">VoIP: (3361 3600) 1053 1173</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a><br>
skype: walmeszeviani<br style="font-family:trebuchet ms,sans-serif"></span><span style="font-family:trebuchet ms,sans-serif">twitter: @walmeszeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">linux user number: 531218</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br>