<div dir="ltr">Oi elias,<div><br></div><div>eu também já tentei deixar o a variável e o ajuste mas não tinha dado certo. Vou tentar novamente!</div><div><br></div><div>Muito obrigada!</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">Em 11 de fevereiro de 2016 18:23, Elias Teixeira Krainski <span dir="ltr"><<a href="mailto:eliaskrainski@yahoo.com.br" target="_blank">eliaskrainski@yahoo.com.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class=""><br>
On 11/02/16 18:43, Michelle Bau Graczyk wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Eu quero os valores de a para o ajuste y=x^a, entende? <br>
</blockquote>
<br></span>
A curva (para baixo ou para cima) não é muito clara devido à presenca de ruído. Note que o modelo estimado em R é<br>
   y_i = f(x_i) + e_i<br>
e no seu caso<br>
   f(x_i) = a - b(x_i-268)^c<br>
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:<br>
<br>
coef(ajuste2<-nls(yy2~a-b*(xx2-268)^c,start=list(a=9,b=1,c=1.4)))<br>
          a           b           c<br>
<a href="tel:8.069130938" value="+18069130938" target="_blank">8.069130938</a> 0.002174558 1.415845808<br>
> mean((yy2-predict(ajuste1))^2)<br>
[1] 0.2871361<br>
> mean((yy2-predict(ajuste2))^2)<br>
[1] 0.2340966<br>
<br>
lines(xx2,((predict(ajuste2))),col="blue",lwd=2)<br>
<br>
Elias<span class=""><br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br></span>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e fornea cdigo mnimo reproduzvel.</blockquote></div><br></div>