[R-es] Parámetros iniciales para ajustes no lineales

[Marcuzzi, Javier Rubén]

Estimado Roy Vera

No sabría exactamente pero:

str(nigra)

edad y ho son numérico, se sobreentiende pero por las dudas.

luego si no leo mal usted cambio el valor de a, es decir start=list(a=

¿sería muy problemático utilizar valores tipo 0.005 (el primero que usted 
usa), 0.01, 0.1, 0.5, 1, 2 ..., usted llega hasta 12.18?

Puede ser que sus valores iniciales están lejos de lo real, ¿en la 
bibliografía hay algún valor como para intentar a partir de este?

O realizar el "mismo modelo" pero lineal, tomar un valor de este y utilizar 
el resultado como valor inicial en el modelo no lineal.

Javier Marcuzzi


-----Mensaje original----- 
From: Roy Vera
Sent: Saturday, November 03, 2012 12:10 PM
To: r-help-es en r-project.org
Subject: [R-es] Parámetros iniciales para ajustes no lineales

Hola a todos

estoy aplicando la función polinómica de Hossfeld [1], y algunos otros 
modelos no lineales para tratar de ajustarlos a un grupo de datos 
forestales,


[1] Y= b*t*exp(c)/(t*exp(c)+a)

Al colocar la función en R con parámetros estimados, me devuelve los 
siguiente:

## model1 <- nls(ho ~ (b*edad*exp(c)/(edad*exp(c)+a)), data=nigra,
    start=list(a=0.005,b=0.08,c=-0.00006), trace=TRUE)##

## 108807.6 :   5e-03  8e-02 -6e-05
Error en nls(ho ~ (b * edad * exp(c)/(edad * exp(c) + a)), data = nigra,  :
  singular gradient ##


O si cambio de valores de los parámetros me suele dar otro tipo de mensaje:

## model1 <- nls(ho ~ (b*edad*exp(c)/(edad*exp(c)+a)), data=nigra,
     start=list(a=12.18,b=0.06,c=-0.00006), trace=TRUE) ##


##109135.3 :  12.18000  0.06000 -0.00006
Error en numericDeriv(form[[3L]], names(ind), env) :
  Missing value or an infinity produced when evaluating the model##

entiendo que esta pregunta va más en sentido de conocimiento de la 
estadística, pero si alguno tiene alguna idea para resolver este problema, 
se la agradecería..

saludos cordiales

Roy
[[alternative HTML version deleted]]







_______________________________________________
R-help-es mailing list
R-help-es en r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es