[R-br] problemas com função nls

[Fernando Antonio de souza]

Veja se agora foi

###cria a função gompertz e chama a função manipulate para encontrar
valores do chute
simple<-as.data.frame(cbind(Manejo=gestsimp$Manejo,Feto=gestsimp$Fetus,Gest=gestsimp$Gest,PESO.FETO=gestsimp$PESO.FETO))
dput(simple)
plot(log(gestsimp$PESO.FETO)~gestsimp$Gest,data=simple)
attach(simple)
gompertz<-function(x,a,b,c){a-b*exp(-c*x)}
gompertz(c(90,110,130,140),1.5,20,0.026)
library(manipulate)
start <- list()
manipulate({plot(log(PESO.FETO)~Gest,data=simple)
            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=2.12,b=14.99,c=0.018),data=simple)
summary(gompertz)#visualiza os parâmetros estimados





Em 10 de março de 2012 20:39, Benilton Carvalho
<beniltoncarvalho em gmail.com>escreveu:

> Fernando, faltou a saida do dput()...
>
> 2012/3/10 Fernando Antonio de souza <nandodesouza em gmail.com>:
> > 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)
> >
> > Em 10 de março de 2012 19:20, Ivan Bezerra Allaman <
> ivanalaman em yahoo.com.br>
> > escreveu:
> >>
> >> 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 Allaman
> >> Universidade Estadual de Santa Cruz
> >> Departamento de Ciências Exatas e Tecnológicas
> >> Ilhéus/BA - Brasil
> >> Fone: +55 73 3680-5076
> >> E-mail: ivanalaman em yahoo.com.br/ivanalaman em gmail.com
> >> @
> >> \end{signature}
> >>
> >> _______________________________________________
> >> R-br mailing list
> >> R-br em 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.
> >
> >
> >
> > _______________________________________________
> > R-br mailing list
> > R-br em 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.
> _______________________________________________
> R-br mailing list
> R-br em 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120310/b50a1399/attachment.html>