
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@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@gmail.com> escreveu:
Fernando, faltou a saida do dput()...
2012/3/10 Fernando Antonio de souza <nandodesouza@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@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@yahoo.com.br/ivanalaman@gmail.com @ \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.
_______________________________________________ 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.
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.
_______________________________________________ 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.