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

Michelle Bau Graczyk mbgraczyk em gmail.com
Quinta Fevereiro 11 21:54:53 BRST 2016


Oi elias,

eu também já tentei deixar o a variável e o ajuste mas não tinha dado
certo. Vou tentar novamente!

Muito obrigada!



Em 11 de fevereiro de 2016 18:23, Elias Teixeira Krainski <
eliaskrainski em yahoo.com.br> escreveu:

>
> 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
> _______________________________________________
> R-br mailing list
> R-br em 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 fornea cdigo
> mnimo reproduzvel.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160211/9bbe659b/attachment.html>


Mais detalhes sobre a lista de discussão R-br