Segue o CMR, os datos encontra-se no dput(). Gostaria de acrescentar algo mais. Meus dados estão bastante dispersos. Pode ser por isto que não está convergindo?Há alguma solução?<br><br>###cria a função gompertz e chama a função manipulate para encontrar valores do chute<br>
plot(log(gestsimp$PESO.FETO)~gestsimp$Gest,data=gestsimp)<br>attach(gestsimp)<br>gompertz<-function(x,a,b,c){a-b*exp(-c*x)}<br>gompertz(c(90,110,130,140),1,5,0.0006)<br>library(manipulate)<br>start <- list()<br>manipulate({plot(log(PESO.FETO)~Gest,data=gestsimp)<br>
curve(gompertz(x,a=a,b=b,c=c),add=TRUE)<br> start <- list(a=a,b=b,c=c)},<br> a=slider(0.00001, 10, initial=0.001),<br> b=slider(-0.0001, 50, initial=0.01),<br> c=slider(0.00007, 0.5, initial=0.0001))<br>
##----------------------------------------------------------------------------<br>#calcula os parâmetros do modelo de gompertz<br>gompertz<-nls(log(PESO.FETO)~a-b*exp(-c*Gest),start=list(a=1.68,b=24.599,c=0.026),subset=c(Gest>0,Fetos==1,Manejo==1),data=tabela)<br>
summary(gompertz)#visualiza os parâmetros estimados<br>dput(gestsimp)<br><br><div class="gmail_quote">Em 10 de março de 2012 19:20, Ivan Bezerra Allaman <span dir="ltr"><<a href="mailto:ivanalaman@yahoo.com.br">ivanalaman@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"><div><div style="font-size:12pt;font-family:arial,helvetica,sans-serif"><div><span>Provavelmente seus chutes estão errados! Acesse o site <a href="http://ridiculas.wordpress.com/" target="_blank">http://ridiculas.wordpress.com/</a> que provavelmente irás resolver o seu problema. Do contrário, nos forneça um CMR.</span></div>
<div><br><span></span></div><div><span>(S,f,P)</span></div><div><span>Allaman<br></span></div><div> </div><div style="background-color:transparent" align="center"><font style="background-color:transparent" face="comic sans ms" size="2"><font size="1"><b><br>
</b></font></font></div><div style="text-align:left;background-color:transparent" align="center"><font style="background-color:transparent"><font face="courier, monaco, monospace, sans-serif">\begin{signature}</font></font></div>
<div style="text-align:left;background-color:transparent" align="center"><font style="background-color:transparent"><font face="courier, monaco,
monospace, sans-serif"><<>>=</font></font></div><div style="text-align:left;background-color:transparent" align="center"><font face="courier, monaco, monospace, sans-serif">Prof. Dr. Ivan Bezerra Allaman</font></div>
<div style="text-align:left;background-color:transparent" align="center"><font face="courier, monaco, monospace, sans-serif">Universidade Estadual de Santa Cruz</font></div><div style="text-align:left;background-color:transparent" align="center">
<font face="courier, monaco, monospace, sans-serif">Departamento de Ciências Exatas e Tecnológicas</font></div><div style="text-align:left;background-color:transparent" align="center"><font face="courier, monaco, monospace, sans-serif">Ilhéus/BA - Brasil</font></div>
<div style="text-align:left;background-color:transparent" align="center"><font face="courier, monaco, monospace, sans-serif">Fone: <a href="tel:%2B55%2073%203680-5076" value="+557336805076" target="_blank">+55 73 3680-5076</a></font></div>
<div style="text-align:left;background-color:transparent" align="center"><font face="courier, monaco, monospace, sans-serif">E-mail: <a href="http://ivanalaman@yahoo.com.br/ivanalaman@gmail.com" target="_blank">ivanalaman@yahoo.com.br/ivanalaman@gmail.com</a></font></div>
<div style="text-align:left;background-color:transparent" align="center"><font style="background-color:transparent"><font face="courier, monaco, monospace, sans-serif">@</font></font></div><div style="text-align:left;background-color:transparent" align="center">
<font style="background-color:transparent"><font face="courier, monaco, monospace, sans-serif">\end{signature}</font></font></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>