PROC NLIN Beispiel: Parameterschätzung für ein nichtlineares Modell

Das Skript beginnt mit der Erstellung eines Datasets namens 'Enzyme', das Konzentrations- und Geschwindigkeitsdaten enthält. Anschließend wird PROC NLIN zweimal ausgeführt. Die erste Ausführung verwendet die Marquardt-Methode, um die Parameter eines Michaelis-Menten-Modells zu schätzen. Die zweite Ausführung von PROC NLIN verwendet die Newton-Methode, um die Parameter eines anderen exponentiellen Modells zu schätzen.

Datenanalyse

Type : CREATION_INTERNE

Die Daten 'Enzyme' werden direkt im Skript über eine DATA STEP-Anweisung und DATALINES erstellt, die Konzentrations- und Geschwindigkeitsbeobachtungen für eine enzymatische Reaktion simulieren.

1 Codeblock

DATA STEP Data

Erklärung :
Dieser DATA STEP-Block erstellt ein Dataset namens 'Enzyme'. Er definiert zwei Variablen, 'Concentration' und 'Velocity', und initialisiert sie mit den direkt im Skript über die DATALINES-Anweisung bereitgestellten Werten. Dieses Dataset wird als Eingabe für die NLIN-Prozeduren verwendet.

Kopiert!

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 Codeblock

PROC NLIN

Erklärung :
Diese NLIN-Prozedur passt ein nichtlineares Modell an die 'Enzyme'-Daten an, indem sie die Marquardt-Optimierungsmethode verwendet. Sie schätzt die Parameter 'theta1' und 'theta2' für das Michaelis-Menten-Modell: `Velocity = theta1*Concentration / (theta2 + Concentration)`. Die Option `hougaard` wird zur Berechnung der Hougaard-Asymmetriemaße verwendet.

Kopiert!

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 Codeblock

PROC NLIN

Erklärung :
Diese zweite NLIN-Prozedur passt ein anderes nichtlineares Modell an dieselben 'Enzyme'-Daten an, diesmal unter Verwendung der Newton-Optimierungsmethode. Sie schätzt die Parameter 'x1' und 'x2' für das exponentielle Modell: `Velocity = x1 * exp (x2 * Concentration)`. Die Option `listcode` zeigt den SAS-Quellcode für das Modell und seine Ableitungen an.

Kopiert!

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

Dieses Material wird von We Are Cas "wie besehen" zur Verfügung gestellt. Es gibt keine ausdrücklichen oder stillschweigenden Garantien hinsichtlich der Marktgängigkeit oder Eignung für einen bestimmten Zweck in Bezug auf die hierin enthaltenen Materialien oder Codes. We Are Cas ist nicht verantwortlich für Fehler in diesem Material, wie es jetzt existiert oder existieren wird, noch bietet We Are Cas technischen Support dafür an.

Urheberrechtsinformationen : S A S S A M P L E L I B R A R Y

Retour à la liste