<div dir="ltr"><div><div>Olá Thiago<br><br></div>Porque você quer outra função para gera chutes iniciais. O manipulate não lhe atende não? <br></div>No blog do Walmes (Ridículas) na página inicial ele fornece um código para utilizar a função rpanel . De uma olhada lá e veja se é isso que você quer.<br>
<a href="http://ridiculas.wordpress.com/">http://ridiculas.wordpress.com/</a><br></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 31 de julho de 2013 06:31, Tiago Souza Marçal <span dir="ltr"><<a href="mailto:tiagosouzamarcal@hotmail.com" target="_blank">tiagosouzamarcal@hotmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div dir="ltr">Obrigado pela dica Augusto,<br> <br>eu conheço a nls(), entretanto, o problema são os chutes iniciais em start=list(), pois, um ajuste pode existir mas dependendo dos chutes iniciais não se consegue a convergência. Com outros modelos que tenho trabalhado a linearização tem dado certo para estimar os chutes iniciais, mas não consegui linearizar este modelo.<div class="im">
<br> <br>Att.<br> <br>Tiago.<br><br><div>#################################################################</div><div> </div><div>Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES</div><div> </div><div>Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas</div>
<div> </div><div>################################################################# </div><br> <br></div><div><hr>Date: Wed, 31 Jul 2013 00:23:26 -0400<br>From: <a href="mailto:ribas.aca@gmail.com" target="_blank">ribas.aca@gmail.com</a><br>
To: <a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a><br>Subject: Re: [R-br] Modelo não linear<div><div class="h5"><br><br><div dir="ltr">Opa, eu não entendo muito bem de regressão não linear, mas talvez a função nls possa fazer o que você tem interesse.<div>
<br></div><div>Tem um livro muito bom da serie use R explicando como usar, vale a pena a consulta caso a função nls se encaixe nas suas necessidades.</div>
<div><a href="http://www.amazon.com/Nonlinear-Regression-R-Use/dp/0387096159" target="_blank">http://www.amazon.com/Nonlinear-Regression-R-Use/dp/0387096159</a><br></div><div><br></div><div>Fora isso tem muitos manuais pela internet também.</div>
<div>Mas mais ou menos é isso que você quer, ajustar esse parâmetros aos dados?</div><div><br></div><div>########################################################################################</div><div>#exemplo</div><div>
tempo=1:12<br></div><div><div>y=c(8.81790512674818,8.83724295007473,8.84652971800299,8.85625726150244,8.86597721721444,8.86647798384983,</div><div> 8.87625601672093,8.89503521901844,8.91478268591035,8.93378059512321,8.97371264813241,9.00355384268932)</div>
<div><br></div><div>plot(y~tempo)</div><div><br></div><div>modelo<-nls(y~A*(1-B*exp(-C*tempo)) + exp(D + E*(tempo-tempo[1])),start=list(A=1,B=1,C=1,D=1,E=1))</div><div>summary(modelo)</div><div><br></div><div>curve(coef(modelo)[1]*(1-coef(modelo)[2]*exp(-coef(modelo)[3]*x)) + exp(coef(modelo)[4] + coef(modelo)[5]*(x-1)),1,12,</div>
<div> add=T,lty=2,col="red")</div></div><div><br></div><div>###########################################################################################</div><div>Mas é um bocado de parâmetros para 12 pontos.</div>
<div>Bem, espero ter ajudado.<br></div><div><br></div><div><br></div></div><div><br><br><div>Em 30 de julho de 2013 19:39, Tiago Souza Marçal <span dir="ltr"><<a href="mailto:tiagosouzamarcal@hotmail.com" target="_blank">tiagosouzamarcal@hotmail.com</a>></span> escreveu:<br>
<blockquote style="padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid">
<div><div dir="ltr">Boa noite pessoal,<br> <br>estou querendo modelar uma curva de embebição de sementes através do seguinte modelo não linear:<br> <br>y=A*(1-B*exp(-C*tempo)) + exp(D + E*(tempo-t0))<br> <br>t0 = 6<br> <br>
Gostaria de saber se existe um procedimento parecido com o do pacote manipulate do Rstudio no R para ajustar este modelo.<br> <br>CMR<br> <br>tempo=1:12<br>y=c(8.81790512674818,8.83724295007473,8.84652971800299,8.85625726150244,<br>
8.86597721721444,8.86647798384983,8.87625601672093,8.89503521901844,8.91478268591035,8.93378059512321,8.97371264813241,9.00355384268932)<br> <br>Att.<br> <br>Tiago.<br><table style="width:576pt;border-collapse:collapse" border="0" cellpadding="0" cellspacing="0" width="768">
<colgroup><col style="width:48pt" span="12" width="64">
</colgroup><tbody><tr style="min-height:15pt" height="20">
<td style="border:0px black;width:48pt;min-height:15pt;background-color:transparent" height="20" align="right" width="64"></td>
<td style="border:0px black;width:48pt;background-color:transparent" align="right" width="64"></td>
<td style="border:0px black;width:48pt;background-color:transparent" align="right" width="64"></td>
<td style="border:0px black;width:48pt;background-color:transparent" align="right" width="64"></td>
<td style="border:0px black;width:48pt;background-color:transparent" align="right" width="64"></td>
<td style="border:0px black;width:48pt;background-color:transparent" align="right" width="64"></td>
<td style="border:0px black;width:48pt;background-color:transparent" align="right" width="64"></td>
<td style="border:0px black;width:48pt;background-color:transparent" align="right" width="64"></td>
<td style="border:0px black;width:48pt;background-color:transparent" align="right" width="64"></td>
<td style="border:0px black;width:48pt;background-color:transparent" align="right" width="64"></td>
<td style="border:0px black;width:48pt;background-color:transparent" align="right" width="64"> <br> <br> <br> <br></td>
<td style="border:0px black;width:48pt;background-color:transparent" align="right" width="64"></td>
</tr>
</tbody></table> <br> <br><br><br> <br><div>#################################################################</div><div> </div><div>Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES</div><div> </div><div>Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas</div>
<div> </div><div>################################################################# </div> </div></div>
<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" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>
<div dir="ltr"><div>Grato<br>Augusto C. A. Ribas</div>
<div> </div>
<div>Site Pessoal: <a href="http://recologia.com.br/" target="_blank">http://recologia.com.br/</a><a href="http://augustoribas.heliohost.org" target="_blank"></a></div><div>Github: <a href="https://github.com/Squiercg" target="_blank">https://github.com/Squiercg</a></div>
<div>Lattes: <a href="http://lattes.cnpq.br/7355685961127056" target="_blank">http://lattes.cnpq.br/7355685961127056</a><br></div></div>
</div>
<br></div></div>_______________________________________________
R-br mailing list
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forne�a c�digo m�nimo reproduz�vel.</div> </div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>