[R-br] Função Gnls

[Walmes Zeviani]

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 em ufpr.br
skype: walmeszeviani
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================


2012/10/1 <rafaelbrco em 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
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121001/f3453cb8/attachment.html>