<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><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.<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><hr id="stopSpelling">Date: Wed, 31 Jul 2013 00:23:26 -0400<br>From: ribas.aca@gmail.com<br>To: r-br@listas.c3sl.ufpr.br<br>Subject: Re: [R-br] Modelo não linear<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 class="ecxgmail_extra"><br><br><div class="ecxgmail_quote">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 class="ecxgmail_quote" 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 width="768" style="width: 576pt; border-collapse: collapse;" border="0" cellspacing="0" cellpadding="0">
<colgroup><col width="64" style="width: 48pt;" span="12">
</colgroup><tbody><tr height="20" style="min-height: 15pt;">
<td width="64" height="20" align="right" style="border: 0px black; width: 48pt; min-height: 15pt; background-color: transparent;"></td>
<td width="64" align="right" style="border: 0px black; width: 48pt; background-color: transparent;"></td>
<td width="64" align="right" style="border: 0px black; width: 48pt; background-color: transparent;"></td>
<td width="64" align="right" style="border: 0px black; width: 48pt; background-color: transparent;"></td>
<td width="64" align="right" style="border: 0px black; width: 48pt; background-color: transparent;"></td>
<td width="64" align="right" style="border: 0px black; width: 48pt; background-color: transparent;"></td>
<td width="64" align="right" style="border: 0px black; width: 48pt; background-color: transparent;"></td>
<td width="64" align="right" style="border: 0px black; width: 48pt; background-color: transparent;"></td>
<td width="64" align="right" style="border: 0px black; width: 48pt; background-color: transparent;"></td>
<td width="64" align="right" style="border: 0px black; width: 48pt; background-color: transparent;"></td>
<td width="64" align="right" style="border: 0px black; width: 48pt; background-color: transparent;"> <br> <br> <br> <br></td>
<td width="64" align="right" style="border: 0px black; width: 48pt; background-color: transparent;"></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">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>_______________________________________________
R-br mailing list
R-br@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 forne�a c�digo m�nimo reproduz�vel.</div> </div></body>
</html>