Exemple PROC NLIN : Estimation de Paramètres pour un Modèle Non Linéaire

Le script commence par créer un dataset nommé 'Enzyme' contenant des données de concentration et de vélocité. Ensuite, il exécute PROC NLIN deux fois. La première exécution utilise la méthode de Marquardt pour estimer les paramètres d'un modèle de Michaelis-Menten. La deuxième exécution de PROC NLIN utilise la méthode de Newton pour estimer les paramètres d'un modèle exponentiel différent.

Analyse des données

Type : CREATION_INTERNE

Les données 'Enzyme' sont créées directement dans le script via une instruction DATA STEP et des DATALINES, simulant des observations de concentration et de vélocité pour une réaction enzymatique.

1 Bloc de code

DATA STEP Data

Explication :
Ce bloc DATA STEP crée un dataset nommé 'Enzyme'. Il définit deux variables, 'Concentration' et 'Velocity', et les initialise avec les valeurs fournies directement dans le script via l'instruction DATALINES. Ce dataset sera utilisé comme entrée pour les procédures NLIN.

Copié !

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 Bloc de code

PROC NLIN

Explication :
Cette procédure NLIN ajuste un modèle non linéaire aux données 'Enzyme' en utilisant la méthode d'optimisation de Marquardt. Elle estime les paramètres 'theta1' et 'theta2' pour le modèle de Michaelis-Menten: `Velocity = theta1*Concentration / (theta2 + Concentration)`. L'option `hougaard` est utilisée pour calculer les mesures d'asymétrie de Hougaard.

Copié !

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 Bloc de code

PROC NLIN

Explication :
Cette deuxième procédure NLIN ajuste un modèle non linéaire différent aux mêmes données 'Enzyme', cette fois en utilisant la méthode d'optimisation de Newton. Elle estime les paramètres 'x1' et 'x2' pour le modèle exponentiel: `Velocity = x1 * exp (x2 * Concentration)`. L'option `listcode` affiche le code source SAS pour le modèle et ses dérivées.

Copié !

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

Ce matériel est fourni "tel quel" par We Are Cas. Il n'y a aucune garantie, expresse ou implicite, quant à la qualité marchande ou à l'adéquation à un usage particulier concernant le matériel ou le code contenu dans les présentes. We Are Cas n'est pas responsable des erreurs dans ce matériel tel qu'il existe maintenant ou existera, et We Are Cas ne fournit pas de support technique pour celui-ci.

Informations de Copyright : S A S S A M P L E L I B R A R Y

Retour à la liste