Veja se agora foi<br><br>###cria a função gompertz e chama a função manipulate para encontrar valores do chute<br>simple<-as.data.frame(cbind(Manejo=gestsimp$Manejo,Feto=gestsimp$Fetus,Gest=gestsimp$Gest,PESO.FETO=gestsimp$PESO.FETO))<br>
dput(simple)<br>plot(log(gestsimp$PESO.FETO)~gestsimp$Gest,data=simple)<br>attach(simple)<br>gompertz<-function(x,a,b,c){a-b*exp(-c*x)}<br>gompertz(c(90,110,130,140),1.5,20,0.026)<br>library(manipulate)<br>start <- list()<br>
manipulate({plot(log(PESO.FETO)~Gest,data=simple)<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=2.12,b=14.99,c=0.018),data=simple)<br>
summary(gompertz)#visualiza os parâmetros estimados<br><br><br><br><br><br><div class="gmail_quote">Em 10 de março de 2012 20:39, Benilton Carvalho <span dir="ltr"><<a href="mailto:beniltoncarvalho@gmail.com">beniltoncarvalho@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Fernando, faltou a saida do dput()...<br>
<br>
2012/3/10 Fernando Antonio de souza <<a href="mailto:nandodesouza@gmail.com">nandodesouza@gmail.com</a>>:<br>
<div class="HOEnZb"><div class="h5">> Segue o CMR, os datos encontra-se no dput(). Gostaria de acrescentar algo<br>
> mais. Meus  dados estão bastante dispersos. Pode ser por isto que não está<br>
> convergindo?Há alguma solução?<br>
><br>
> ###cria a função gompertz e chama a função manipulate para encontrar valores<br>
> 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>
> Em 10 de março de 2012 19:20, Ivan Bezerra Allaman <<a href="mailto:ivanalaman@yahoo.com.br">ivanalaman@yahoo.com.br</a>><br>
> escreveu:<br>
>><br>
>> Provavelmente seus chutes estão errados! Acesse o site<br>
>> <a href="http://ridiculas.wordpress.com/" target="_blank">http://ridiculas.wordpress.com/</a> que provavelmente irás resolver o seu<br>
>> problema. Do contrário, nos forneça um CMR.<br>
>><br>
>> (S,f,P)<br>
>> Allaman<br>
>><br>
>><br>
>> \begin{signature}<br>
>> <<>>=<br>
>> Prof. Dr. Ivan Bezerra Allaman<br>
>> Universidade Estadual de Santa Cruz<br>
>> Departamento de Ciências Exatas e Tecnológicas<br>
>> Ilhéus/BA - Brasil<br>
>> Fone: <a href="tel:%2B55%2073%203680-5076" value="+557336805076">+55 73 3680-5076</a><br>
>> E-mail: <a href="http://ivanalaman@yahoo.com.br/ivanalaman@gmail.com" target="_blank">ivanalaman@yahoo.com.br/ivanalaman@gmail.com</a><br>
>> @<br>
>> \end{signature}<br>
>><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<br>
>> código mínimo reproduzível.<br>
><br>
><br>
><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<br>
> mínimo reproduzível.<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>
</div></div></blockquote></div><br>