[R-br] Erro no ajuste usando nls()

[Elias Teixeira Krainski]

On 11/02/16 18:43, Michelle Bau Graczyk wrote:
> Eu quero os valores de a para o ajuste y=x^a, entende? 

A curva (para baixo ou para cima) não é muito clara devido à presenca de 
ruído. Note que o modelo estimado em R é
    y_i = f(x_i) + e_i
e no seu caso
    f(x_i) = a - b(x_i-268)^c
sendo que você fixou 'a' é um parâmetro de limiar e você fixou em 9.548. 
Porém, note que limiar nesse modelo é uma média, não um valor máximo. 
Assim, observando os dados me parece que esse valor é muito alto... 
Considerando a possibilidade de estimar o limiar 'a', o resultado é uma 
curva com concavidade invertica (c>1) e cujo ajuste tem erro quadrático 
médio menor. Veja:

coef(ajuste2<-nls(yy2~a-b*(xx2-268)^c,start=list(a=9,b=1,c=1.4)))
           a           b           c
8.069130938 0.002174558 1.415845808
 > mean((yy2-predict(ajuste1))^2)
[1] 0.2871361
 > mean((yy2-predict(ajuste2))^2)
[1] 0.2340966

lines(xx2,((predict(ajuste2))),col="blue",lwd=2)

Elias