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_con... 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