Ejemplo PROC NLIN: Estimación de Parámetros para un Modelo No Lineal

El script comienza creando un conjunto de datos llamado 'Enzyme' que contiene datos de concentración y velocidad. Luego, ejecuta PROC NLIN dos veces. La primera ejecución utiliza el método de Marquardt para estimar los parámetros de un modelo de Michaelis-Menten. La segunda ejecución de PROC NLIN utiliza el método de Newton para estimar los parámetros de un modelo exponencial diferente.

Análisis de datos

Type : CREATION_INTERNE

Los datos 'Enzyme' se crean directamente en el script a través de una instrucción DATA STEP y DATALINES, simulando observaciones de concentración y velocidad para una reacción enzimática.

1 Bloque de código

DATA STEP Data

Explicación :
Este bloque DATA STEP crea un conjunto de datos llamado 'Enzyme'. Define dos variables, 'Concentration' y 'Velocity', y las inicializa con los valores proporcionados directamente en el script a través de la instrucción DATALINES. Este conjunto de datos se utilizará como entrada para los procedimientos NLIN.

¡Copiado!

1	DATA Enzyme;
2	INPUT Concentration Velocity;
3	DATALINES;
4	0.26 124.7 0.30 126.9
5	0.48 135.9 0.50 137.6
6	0.54 139.6 0.68 141.1
7	0.82 142.8 1.14 147.6
8	1.28 149.8 1.38 149.4
9	1.80 153.9 2.30 152.5
10	2.44 154.5 2.48 154.7
11	;

2 Bloque de código

PROC NLIN

Explicación :
Este procedimiento NLIN ajusta un modelo no lineal a los datos 'Enzyme' utilizando el método de optimización de Marquardt. Estima los parámetros 'theta1' y 'theta2' para el modelo de Michaelis-Menten: `Velocity = theta1*Concentration / (theta2 + Concentration)`. La opción `hougaard` se utiliza para calcular las medidas de asimetría de Hougaard.

¡Copiado!

1	PROC NLIN DATA=Enzyme method=marquardt hougaard;
2	parms theta1=155
3	theta2=0 to 0.07 BY 0.01;
4	model Velocity = theta1*Concentration / (theta2 + Concentration);
5	RUN;

3 Bloque de código

PROC NLIN

Explicación :
Este segundo procedimiento NLIN ajusta un modelo no lineal diferente a los mismos datos 'Enzyme', esta vez utilizando el método de optimización de Newton. Estima los parámetros 'x1' y 'x2' para el modelo exponencial: `Velocity = x1 * exp (x2 * Concentration)`. La opción `listcode` muestra el código fuente SAS para el modelo y sus derivadas.

¡Copiado!

1	PROC NLIN DATA=Enzyme method=newton listcode;
2	parms x1=4 x2=2;
3	model Velocity = x1 * exp (x2 * Concentration);
4	RUN;

Este material se proporciona "tal cual" por We Are Cas. No hay garantías, expresas o implícitas, en cuanto a la comerciabilidad o idoneidad para un propósito particular con respecto a los materiales o el código contenidos en este documento. We Are Cas no es responsable de los errores en este material tal como existe ahora o existirá, ni We Are Cas proporciona soporte técnico para el mismo.

Información de copyright : S A S S A M P L E L I B R A R Y

Retour à la liste