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

[Benilton Carvalho]

Ainda nao foi dessa vez... mas esta' facil.. Execute o seu codigo,
dai' veja no R o que vai aparecer... vai aparecer algo esquisito como
abaixo:

> dput(data.frame(x=rnorm(10)))
structure(list(x = c(-0.543332683249016, -1.72683892829404, -1.46539343776498,
-1.39513469551024, 0.809122589671721, -1.58459611384001, -0.818826156554395,
-0.391288255253397, -0.579453920817577, 0.473369713399808)), .Names =
"x", row.names = c(NA,
-10L), class = "data.frame")


dai' vc copia pra gente esse resultado estranho (q no meu caso comeca
com "structure")... :)

b

2012/3/11 Fernando Antonio de souza <nandodesouza em gmail.com>:
> 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.
>
>
>
> _______________________________________________
> 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.