On 03/10/2012 06:16 PM, Fernando Antonio de souza wrote:
caros amigos estou tentando rodar um modelo não linear de Gormpertz da
forma Log(Y) = a-b*exp(-c*X) utiliando a função nls. A função é como segue:
caros amigos estou tentando rodar um modelo não linear de Gormpertz da
forma Log(Y) = a-b*exp(-c*X) utiliando a função nls. A função é como segue:
gompertz-nls(log(PESO.FETO)~a-b*exp(-c*Gest),start=list(a=3.4,b=22.99,c=0.018),data=gestsimp)
dput(gestsimp)
Acontece que toda vez que rodo, recebo o
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
Fernando, faltou a saida do dput()...
2012/3/10 Fernando Antonio de souza nandodeso...@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
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,
note q nao entro na discussao da validade do modelo mas a questao
eh a selecao dos valores iniciais... b
blah = simple[complete.cases(simple),]
nls(log(PESO.FETO)~a-b*exp(-c*Gest),start=list(a=1.68,b=-24.599,c=-0.026),
data=blah)
curve(-10.48+6.79*exp(.003915*x), 90, 140)
with(blah,
Muito obrigado Benilton! Estou iniciando a estudar modelos não lineares e
pequenos problemas de início parece se tornar grande. Mas eu fiquei com uma
dúvida agora. Quando eu logaritmizo o modelo de Gompertz eu o torno linear?
Pergunto isto porque no gráfico o modelo assume a forma de uma reta e
O modelo q vc empregou a-b*exp(c*X) nao e' linearizavel por meio de
log() para valores gerais (seria, por exemplo, com a=0). O fato de o
grafico parecer uma reta, nao quer dizer que seja uma reta... Note que
se vc fizer um grafico de uma parabola e der um zoom bastante grande,
eventualmente vc