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 < [email protected]> 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 > [email protected] > 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.
_______________________________________________ R-br mailing list [email protected] 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.
