Beispiel für PROC NLIN mit Parameterprofilierung und Bootstrap

Das Skript passt ein exponentielles Modell an einen Datensatz an. Es verwendet die PROC NLIN, um nicht nur die Modellparameter zu schätzen, sondern auch deren Stabilität und Empfindlichkeit zu bewerten. Die PROFILE-Anweisung wird verwendet, um die Profil-Likelihood von zwei Parametern zu untersuchen, während die BOOTSTRAP-Anweisung eine Resampling-Analyse durchführt, um robuste Konfidenzintervalle und grafische Diagnosen für die Schätzer zu liefern.

Datenanalyse

Type : CREATION_INTERNE

Der Datensatz 'clarke1987a' wird direkt im Skript mit einem DATA-Schritt und einer 'datalines'-Anweisung erstellt.

1 Codeblock

DATA STEP Data

Erklärung :
Dieser Codeblock erstellt die SAS-Tabelle 'clarke1987a' in der WORK-Bibliothek. Die Tabelle enthält zwei numerische Variablen, 'x' und 'y', deren Daten direkt aus dem Eingabestrom über die 'datalines'-Anweisung gelesen werden.

Kopiert!

1	DATA clarke1987a;
2	INPUT x y;
3	DATALINES;
4	1 3.183
5	2 3.059
6	3 2.871
7	4 2.622
8	5 2.541
9	6 2.184
10	7 2.110
11	8 2.075
12	9 2.018
13	10 1.903
14	11 1.770
15	12 1.762
16	13 1.550
17	;
18

2 Codeblock

PROC NLIN

Erklärung :
Dieser Block aktiviert die Generierung von ODS-Grafiken und führt dann die nichtlineare Regressionsprozedur (NLIN) aus. Er definiert die Anfangswerte der Parameter (parms), spezifiziert ein exponentielles Modell (model) und fordert zwei erweiterte Analysen an: eine Profil-Likelihood für 'theta1' und 'theta3' (profile) und eine Bootstrap-Analyse mit 2000 Resampling-Durchläufen, um Konfidenzintervalle und Diagnosegrafiken zu erhalten (bootstrap). Die ODS-Grafiken werden anschließend deaktiviert.

Kopiert!

1	ods graphics on;
2	PROC NLIN DATA=clarke1987a plots(stats=none)=diagnostics;
3	parms theta1=-0.15
4	theta2=2.0
5	theta3=0.80;
6	profile theta1 theta3 / range = -6 to 2 BY 0.2 all;
7	bootstrap / nsamples = 2000 seed=123 bootplots bootci bootcov;
8	model y = theta3 + theta2exp(theta1x);
9	RUN;
10	ods graphics off;

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