<div dir="ltr"><div class="gmail_default" style="font-family:'trebuchet ms',sans-serif">Você quer dizer que ajustou os polinômios ortogonais pelo R e fez isso na mão também usando os coeficientes dessas tabelas de final de livro? Pois bem, confira os coeficientes usados com os do R, via model.matrix(). Existem infinitas bases ortogonais que poder sem usadas. As vezes o R usa uma ortonormal e a sua feita na mão não é, ou vice versa. Independente do fator de escala (produto interno 1 ou não), espera-se que a significância não mude, pois só depende da ortogonalidade. Outra coisa que não irá mudar são os valores preditos.</div>
<div class="gmail_default" style="font-family:'trebuchet ms',sans-serif"><br></div><div class="gmail_default"><font face="courier new, monospace">dados <- structure(list(Bloco = structure(c(1L, 1L, 1L, 1L, 1L, 2L, 2L,</font></div>
<div class="gmail_default"><font face="courier new, monospace">2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L,</font></div><div class="gmail_default"><font face="courier new, monospace">5L, 6L, 6L, 6L, 6L, 6L), .Label = c("1", "2", "3", "4", "5", "6"),</font></div>
<div class="gmail_default"><font face="courier new, monospace">class = "factor"), Espacamento = c(1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2,</font></div><div class="gmail_default"><font face="courier new, monospace">3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5), Producao =</font></div>
<div class="gmail_default"><font face="courier new, monospace">c(3.02, 2.795, 2.966, 2.552, 2.822, 3.334, 3.1, 2.651, 2.687, 2.543,</font></div><div class="gmail_default"><font face="courier new, monospace">3.065, 2.741, 2.624, 2.84, 2.597, 3.109, 2.85, 2.759, 2.903, 2.57,</font></div>
<div class="gmail_default"><font face="courier new, monospace">3.181, 2.759, 2.759, 2.525, 2.66, 3.244, 2.723, 2.507, 2.417, 3.002)),</font></div><div class="gmail_default"><font face="courier new, monospace">.Names = c("Bloco", "Espacamento", "Producao"), row.names = c(NA,</font></div>
<div class="gmail_default"><font face="courier new, monospace">-30L), class = "data.frame")</font></div><div class="gmail_default"><font face="courier new, monospace"><br></font></div><div class="gmail_default">
<font face="courier new, monospace"># igualmente espaçado?</font></div><div class="gmail_default"><font face="courier new, monospace">diff(sort(unique(dados$Espacamento)))</font></div><div class="gmail_default"><font face="courier new, monospace"><br>
</font></div><div class="gmail_default"><font face="courier new, monospace"># 6 repetições</font></div><div class="gmail_default"><font face="courier new, monospace">xtabs(~Espacamento, dados)</font></div><div class="gmail_default">
<font face="courier new, monospace"><br></font></div><div class="gmail_default"><font face="courier new, monospace">dados$Bloco <- factor(dados$Bloco)</font></div><div class="gmail_default"><font face="courier new, monospace">dados$Espacamento <- factor(dados$Espacamento)</font></div>
<div class="gmail_default"><font face="courier new, monospace">modelo3 <- lm(Producao~Bloco+ordered(Espacamento), data=dados,</font></div><div class="gmail_default"><font face="courier new, monospace">              contrast=list(Bloco=contr.sum))</font></div>
<div class="gmail_default"><font face="courier new, monospace">summary(modelo3)</font></div><div class="gmail_default"><font face="courier new, monospace"><br></font></div><div class="gmail_default"><font face="courier new, monospace">nlevels(dados$Espacamento)</font></div>
<div class="gmail_default"><font face="courier new, monospace"><br></font></div><div class="gmail_default"><font face="courier new, monospace">X <- model.matrix(modelo3)</font></div><div class="gmail_default"><font face="courier new, monospace">i <- modelo3$assign%in%c(0,2)</font></div>
<div class="gmail_default"><font face="courier new, monospace"><br></font></div><div class="gmail_default"><font face="courier new, monospace"># os coeficientes usados</font></div><div class="gmail_default"><font face="courier new, monospace">round(unique(X[,i]), 3)</font></div>
<div class="gmail_default"><font face="courier new, monospace"><br></font></div><div class="gmail_default"><font face="courier new, monospace"># o gráficos do coeficientes, que inclusive é muito esclarecedor</font></div><div class="gmail_default">
<font face="courier new, monospace">matplot(1:5, unique(X[,i]), type="o")</font></div><div class="gmail_default"><font face="courier new, monospace"><br></font></div><div class="gmail_default"><font face="courier new, monospace"># 6 na diagonal, 6 repetições</font></div>
<div class="gmail_default"><font face="courier new, monospace">round(t(X[,i])%*%X[,i],2)</font></div><div class="gmail_default"><font face="courier new, monospace"><br></font></div><div class="gmail_default"><div class="gmail_default">
<font face="courier new, monospace">x <- ordered(gl(5,1))</font></div><div class="gmail_default"><font face="courier new, monospace">X <- model.matrix(~x)</font></div><div class="gmail_default"><font face="courier new, monospace">round(t(X)%*%X,2) # base ortonormal</font></div>
<div><br></div></div><div class="gmail_default"><font face="trebuchet ms, sans-serif">Recentemente saiu essa matéria no R bloggers.</font></div><div class="gmail_default"><font face="trebuchet ms, sans-serif"><a href="http://appliedpredictivemodeling.com/blog/2013/10/23/the-basics-of-encoding-categorical-data-for-predictive-models">http://appliedpredictivemodeling.com/blog/2013/10/23/the-basics-of-encoding-categorical-data-for-predictive-models</a></font><span style="font-family:'trebuchet ms',sans-serif">  </span></div>
<div class="gmail_default"><span style="font-family:'trebuchet ms',sans-serif"><br></span></div><div class="gmail_default"><span style="font-family:'trebuchet ms',sans-serif">À disposição.</span></div><div class="gmail_default">
<span style="font-family:'trebuchet ms',sans-serif">Walmes.</span></div><div class="gmail_extra"><br clear="all"><div><div dir="ltr"><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">skype: walmeszeviani<br style="font-family:'trebuchet ms',sans-serif"></span><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></div>
</div></div></div>