O problema que voce relata é pq o polinomio nao foi ajustado com x, x^2 ... x^5 e sim pelo geado pelo poly() POnrtanto seus coneficientes nao sao compatívies com as potencias de x na forma "bruta" Obtenbha os falores ajustados por fitted(r) e prediçoes em outros pontos usando predict() para fazer a predicão On Wed, 25 Jul 2012, Tito Conte wrote:
as dicas foram uteis e muito obrigado, muito embora quando eu dou o comando
r=lm(y~poly(x,5))
e depois dou a equação pelos coeficientes
b0=r$coefficient[1] b1=r$coefficient[2] b2=r$coefficient[3] b3=r$coefficient[4] b4=r$coefficient[5] b5=r$coefficient[6]
m_teste=b0+b1*x+b2*x^2+b3*x^3+b4*x^4+b5*x^5
o m_teste extrapola de novo
será que estou extraindo os coeficientes de forma correta? e aplicando a equação de forma correta?
vocês podem me dar outra luz?
obrigado