Rafael,
Talvez você esteja passando chutes iníciais ruins. Faça gráficos de dispersão da resposta conta as covariáveis para obter chutes melhores ou use outro meio de pré-estimação.
Envia sua resposta para lista de discussão do R, R-br, pois mais pessoas podem aproveitar-se da mesma dúvida, mais pessoas podem colaborar com você, você pode ser respondido em menos tempo, você pode contemplar diferentes soluções para o seu problema e sua mensagem fica no histórico da lista para consulta futura. Detalhes de inscrição na lista disponíveis em
http://www.leg.ufpr.br/doku.php/software:rbr
Não esqueça de ler o guia de postagem e submeter um CMR na sua mensagem disponíveis em
http://www.leg.ufpr.br/doku.php/software:rbr#guia_de_postagem_da_r-br_em_construcao
https://gist.github.com/1088208
À disposição.
Walmes.
==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: walmes@ufpr.br
skype: walmeszeviani
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
2012/10/1
<rafaelbrco@hotmail.com>
Olá, então to com o seguinte problema no R:
Tenho este modelo abaixo:
Gompertz=nls(Ht~SSgompertz(Dap,b0,b1,b2))
summary(Gompertz)
Com os seguintes resultados para os seus parâmetros
Formula: Ht ~ SSgompertz(Dap, b0, b1, b2)
Parameters:
Estimate Std. Error t value Pr(>|t|)
b0 16.35440 0.98953 16.53 <2e-16 ***
b1 1.67231 0.67153 2.49 0.0133 *
b2 0.81996 0.05814 14.10 <2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 4.704 on 310 degrees of freedom
Number of iterations to convergence: 7
Achieved convergence tolerance: 9.538e-07
Porem quando incluo a covariavel Idade (I) aparece o seguinte erro:
Gompertz.covariavel=gnls(Ht~SSgompertz(Dap,b0,b1,b2),params=list(b0~I,b1~I,b2~I),start=c(b0=16,0,b1=2,0,b2=1,0))
Erro em gnls(Ht ~ SSgompertz(Dap, b0, b1, b2), params = list(b0 ~ I, :
Step halving factor reduced below minimum in NLS step
Alguem poderia me dar algumas dicas