     Hola Carlos,

El 24/01/12 15:54, Carlos Ortega escribió:
> Hola Argel,
>
> Como referencia adicional, por si puede dar algún resultado o ayudar a 
> enfocar el problema de una manera alternativa, hay un libro "Non 
> Linear Regression with R" que también tiene un paquete asociado en R 
> (nlrwr):
>
> http://www.amazon.es/Nonlinear-Regression-R-Use/dp/0387096159/ref=sr_1_sc_1?ie=UTF8&qid=1327441411&sr=8-1-spell 
> <http://www.amazon.es/Nonlinear-Regression-R-Use/dp/0387096159/ref=sr_1_sc_1?ie=UTF8&qid=1327441411&sr=8-1-spell> 
>
> http://cran.at.r-project.org/web/packages/nlrwr/index.html
>
> En este libro también se estudia el caso de esta ecuación y se 
> contempla lo que has comentado de P0 con diferentes valores al que se 
> le puede seguir aplicando la función nls() pero especificando que P0 
> indica realmente un factor. Extender el conjunto de datos con el caso 
> que comentas de P0=5 podría ayudar al modelo.

     Acabo de conseguir el libro, pero apenas lo voy comenzando a leer. 
Estudiaré con detalle la sección que comentas.

     Me queda la siguiente duda al tratar de ajustar el modelo mixto con 
nls, ¿sería válido quedarse con el resultado final que se obtiene con 
nls y "port", aunque el algoritmo señale que existe "Convergence 
failure: singular convergence (7)"?

     Lo anterior considerando que los valores de Km, Kic y k3 obtenidos 
con ambos sistemas de análisis son muy similares (excepto por Kiu que en 
SAS es muy grande y que en nls no converge pero aumenta con tendencia 
hacia valores muy grandes) la interpretación sería la misma que la de 
los autores, donde la constante Kiu que alcanza valores tan grandes hace 
que realmente carezca de sentido incluirla para modelar la reacción 
enzimática en estudio... quedando sólo las constantes Km, Kic y k3 en el 
modelo que ahora sería de inhibición competitiva.

> De todas formas llegado un punto de bloqueo también puedes escribir a 
> los autores del paper y preguntarles por detalles de cómo codificaron 
> su modelo, si han llegado o pueden utilizar R para replicar 
> resultados. El impacto de su paper ya está conseguido, vaya supongo 
> que te darían acceso a lo que han hecho...
>

     Buena idea, espero obtener respuesta y poder complementar todo lo 
que hasta ahora se ha discutido en este hilo.

> Saludos,
> Carlos Ortega
> www.qualityexcellence.es <http://www.qualityexcellence.es>

     De nuevo, muchas gracias por tus sugerencias.

     Saludos.

--
     Argel.

	[[alternative HTML version deleted]]

