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?
###cria a função gompertz e chama a função manipulate para encontrar valores do chute
plot(log(gestsimp$PESO.FETO)~gestsimp$Gest,data=gestsimp)
attach(gestsimp)
gompertz<-function(x,a,b,c){a-b*exp(-c*x)}
gompertz(c(90,110,130,140),1,5,0.0006)
library(manipulate)
start <- list()
manipulate({plot(log(PESO.FETO)~Gest,data=gestsimp)
curve(gompertz(x,a=a,b=b,c=c),add=TRUE)
start <- list(a=a,b=b,c=c)},
a=slider(0.00001, 10, initial=0.001),
b=slider(-0.0001, 50, initial=0.01),
c=slider(0.00007, 0.5, initial=0.0001))
##----------------------------------------------------------------------------
#calcula os parâmetros do modelo de gompertz
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)
summary(gompertz)#visualiza os parâmetros estimados
dput(gestsimp)
Provavelmente seus chutes estão errados! Acesse o site http://ridiculas.wordpress.com/ que provavelmente irás resolver o seu problema. Do contrário, nos forneça um CMR.(S,f,P)Allaman
\begin{signature}<<>>=Prof. Dr. Ivan Bezerra AllamanUniversidade Estadual de Santa CruzDepartamento de Ciências Exatas e TecnológicasIlhéus/BA - BrasilFone: +55 73 3680-5076@\end{signature}
_______________________________________________
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.