[R-br] Bons valores iniciais e não convergência - regressão não linear

Adriele Pereira adriele_tutoria em yahoo.com.br
Terça Março 8 10:28:31 BRT 2016


Olá, bom dia a todos.
Estou tentando ajustar o modelo não linear Richards ao seguinte conjunto de dados:  y<-c(32.2,45.1,63.3,87.7,107.6,120.4,131.8,148.9,159.3,166.8,173.4,187.4,195.5,198.8,211.0,221.4,226.6,231.2,237.7,243.1,246.6,250.3)
  x<-c(90,180,270,360,450,540,630,720,810,900,990,1080,1170,1260,1350,1440,1530,1620,1710,1800,1890,1980)

Ao utilizar o shiny, obtive os seguintes valores iniciais: alpha=270 , beta= 402, k=0.0016 , d=0.8. O problema é que o gráfico fica muito bom com esses valores, mas ao colocar esse valores iniciais na nls aparece a mensagem de erro:
 ajuste=nls(y~(alpha*(1 + (d - 1)*exp(-k*(x-beta)))^(1/(1-d))),
start=c(alpha= 270, beta= 402 , k= 0.0016, d=0.8))
Error in numericDeriv(form[[3L]], names(ind), env) : 
  Obtido valor faltante ou infinito quando avaliando o modelo


Já tentei mudar a parametrização do modelo e aumentar o número de iterações e mesmo assim não ajustou. Alguém já se deparou com um problema parecido? 

Desde já agradeço.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160308/7ab02251/attachment.html>


Mais detalhes sobre a lista de discussão R-br