Sem testar código segue um comentário genérico:
Estes polinomios de alto grau podem gerar variáveis de valores muito
elevados e com isto instabiliades numéricas nos calculos matriciais
tente usar polinomios ortonormalizados fazendo:
y ~ poly(x, 5)
Os coeficientes serão diferentes mas
Leonardo, fiz a troca e não consegui rodar. Muito obrigado.
Edson Lira
Estatístico
Manaus-Amazonas
De: Leonard de Assis assis.leon...@gmail.com
Para: r-br@listas.c3sl.ufpr.br
Enviadas: Terça-feira, 24 de Julho de 2012 13:56
Assunto: Re: [R-br] Tabelas em sql
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
Tito
seu exemplo ilustra um o problema numérico
que pode ocorrer em ajustes.
O perigo é que o resultado é prouzido e se não avaliao o erro pode passar
desapercebido e a curva de ajuste erraa pode acabar sendo usada.
Veja o efeito no gráfico.
Teoricamente os ajustes deveriam ser iguais,
Paulo agora outra dúvida, tenho dez mil pontos, como comparar os que estão
a direita e a esquerda dessa curva (como se fosse um dentro ou fora do
sistema), já que não consigo extrair o coeficiente de da função poly?
Necessito dos dados numericamente pois estes pontos serão submetidos a
novos
Seja mais claro na sua colocação e forneça um CMR ilustrando com
comentários os pontos em dúvida.
À disposição.
Walmes.
==
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Walmes vou criar um exemplo para facilitar o entendimento
o que eu preciso é comparar pontos amostrados com os resultados de um
modelo de interpolação polinomial
por exemplo
#minhas amostras
y=seq(1:9)
x=seq(3:11)
#pontos do modelo
xm=seq(2,11)
ym=seq(0:8)
#tirando os valores do modelo da