<div><font face="trebuchet ms,sans-serif">Para ter esses efeitos você tem que usar o contraste tipo polinomial para o seu fator. Tem diversas formas de fazer isso. Como você não deu um CMR eu tô mandando um com lm() mesmo.</font></div>
<div><font face="trebuchet ms,sans-serif"><br></font></div><div><div><font face="'courier new', monospace">x <- factor(rep(1:4,3))</font></div><div><font face="'courier new', monospace">y <- rnorm(x)</font></div>
<div><font face="'courier new', monospace"><br></font></div><div><font face="'courier new', monospace"># convertendo dentro</font></div><div><font face="'courier new', monospace">m0 <- lm(y~x, contrast=list(x=contr.poly))</font></div>
<div><font face="'courier new', monospace">coef(m0)</font></div><div><font face="'courier new', monospace"><br></font></div><div><font face="'courier new', monospace">m0 <- lm(y~x)</font></div><div>
<font face="'courier new', monospace">coef(m0)</font></div><div><font face="'courier new', monospace"><br></font></div><div><font face="'courier new', monospace"># convertendo fora</font></div><div>
<font face="'courier new', monospace">contrasts(x) <- contr.poly(nlevels(x))</font></div><div><font face="'courier new', monospace">m0 <- lm(y~x)</font></div><div><font face="'courier new', monospace">coef(m0)</font></div>
<div style="font-family:'trebuchet ms',sans-serif"><br></div><div style="font-family:'trebuchet ms',sans-serif">Lembrando que usar polinômios ortogonais só tem sentido para **níveis igualmente espaçados** do fator. Porque eles são/eram úteis? Porque a matriz X'X fica ortogonal e dessa forma você faz inferência para um efeito marginal aos demais (embora não faça o mínimo de sentido julgar o termo linear depois do quadrático etc). Na minha opinião os testes de hipótese sequênciais são mais apropriados e não dependem dessa condição restritiva de níveis igualmente espaçados.</div>
<div style="font-family:'trebuchet ms',sans-serif"><br></div><div style="font-family:'trebuchet ms',sans-serif">À disposição.</div><div style="font-family:'trebuchet ms',sans-serif">Walmes.</div><div style="font-family:'trebuchet ms',sans-serif">
<br></div></div><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></span><br style="font-family:trebuchet ms,sans-serif"><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>