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
_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
