Hola,

Sí, prueba el paquete nls2 y dentro del paquete la función nls2().
Básicamente utiliza el mismo algoritmo que la función nls(), pero evaluado
en un intervalo que tú defines.
Mira los ejemplos de la función.

Si sigues con errores, creo que entonces es mejor que nos dejes el conjunto
de datos que quieres ajustar para hacer otras pruebas...

Saludos,
Carlos Ortega
www.qualityexcellence.es

El 3 de noviembre de 2012 20:45, Roy Vera <roy_2roger@yahoo.com.mx>escribió:

> Hola Javier
>
> gracias por tu respuesta, bueno la verdad es que he estado probando con
> varios valores de los parámetros, naturalmente que en la bibliografía
> existen parámetros al cual puedo referirme para poder iniciar, de hecho lo
> he intentado con algunos, pero igual me da el mismo error, por esa razón
> quisiera saber si existe otro método de encontrar los parámetros de
> partida...y también entender qué significa lo que me devuelve R como
> error...
>
>
> algunos datos adicionales:
>
> ho = variable numérica referida a altura dominante
> edad = variable edad numérica
> data = nigra (nombre del dataframe)
>
> saludos
>
> Roy
>
>
>
>
>
>
> ________________________________
>  De: "Marcuzzi, Javier Rubén" <javier.ruben.marcuzzi@gmail.com>
> Para: Roy Vera <roy_2roger@yahoo.com.mx>; r-help-es@r-project.org
> Enviado: Sábado, 3 de noviembre, 2012 18:13:44
> Asunto: Re: [R-es] Parámetros iniciales para ajustes no lineales
>
> 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@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@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>         [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>


-- 
Saludos,
Carlos Ortega
www.qualityexcellence.es

	[[alternative HTML version deleted]]

