     Hola Carlos, muchas gracias por seguir revisando el tema. Pongo 
algunas notas abajo.

El 27/01/12 04:24, Carlos Ortega escribió:
> Hola Argel,
>
> He continuado mirando este tema, pero buscando la función que hay que 
> ajustar podía expresarse analíticamente de otra manera. También he 
> mirado de qué forma queda expresada la función si se considera que P0 
> es cero y que Kiu tiende a infinito.
>

     Aquí no creo conveniente considerar que siempre P0 = 0 para 
simplificar el modelo, porque en realidad puede tomar valores iniciales. 
No lo puse en los datos que envié inicialmente porque el ajuste que no 
me daba convergencia era con este set de datos donde P0 = 0 (en el que 
los autores reportan que sí converge), pero el set de datos completo 
incluye una serie de 7 experimentos adicionales donde P0 = 5, y el 
análisis completo involucra el ajuste con los dos grupos de datos.

> Pero en el camino me he encontrado con un detalle que no sé si es un 
> error o es algo que no lo estoy yo haciendo bien.
> En la expresión:
>
> ((Km+Km*P0/Kic+Km*S0/Kic)*log(
> (S0-(Pt-P0))/S0)+(-Km/Kic+1+P0/Kiu+S0/Kiu)*((S0-(Pt-P0))-S0)+(-1/(2*Kiu))*((S0-(Pt-P0))^2-S0^2))*(-1/(E*k3))
>
> Mira el término, que aunque tiene muchos paréntesis:
>
> ((S0-(Pt-P0))-S0)
>
> Esto desarrollado es:
> ( (S0 - Pt + P0) - S0) = ( S0 - Pt + P0 - S0 ) = (P0 - Pt)    .... ¿Se 
> cancelan los S0?
>

     Tienes razón, se cancelan los S0 y queda sólo (P0-Pt). Así reportan 
los autores que ingresaron el modelo en SAS. Ya lo había notado y lo 
cambié para ingresarlo en nls dejando sólo (P0-Pt) y el resultado es el 
mismo, no converge, aunque los valores finales de "port" aunque son 
similares a los reportados por los autores con SAS, al parecer ajustan 
mejor, de acuerdo con las sumas de cuadrados que comentaba en el mensaje 
anterior.

> No he probado a ajustar la función si esto se confirma, por si hay 
> cambios en el ajuste nls...
>
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es <http://www.qualityexcellence.es>

     Saludos.

--
     Argel.

	[[alternative HTML version deleted]]

