Startbeispiel für PROC NLMIXED (nlmgs1)

Dieses SAS^©-Programm verwendet die NLMIXED-Prozedur, um eine logistische Wachstumskurve zu modellieren. Die Daten, die aus Draper und Smith (1981) über das Wachstum von Orangenbäumen stammen, werden online über einen DATA STEP erstellt. Das Modell enthält feste Parameter (b1, b2, b3) und einen Zufallseffekt (u1) auf Baumebene (Subjekt), um die Variabilität zwischen den Bäumen zu erfassen. Das Modell ist eine Normalverteilung, bei der der Mittelwert durch eine logistische Funktion und die Varianz durch s2e definiert wird.

Datenanalyse

Type : INTERNE_ERSTELLUNG

Die Daten mit dem Namen 'tree' werden direkt im SAS-Skript unter Verwendung einer DATALINES-Anweisung erstellt, die die Variablen 'tree', 'day' und 'y' enthält.

1 Codeblock

DATA STEP Data

Erklärung :
Dieser DATA STEP-Block erstellt ein SAS-Dataset namens 'tree'. Er verwendet die INPUT-Anweisung, um die Variablen 'tree' (Baum-ID), 'day' (Tag) und 'y' (gemessener Wert) zu definieren. Die Daten sind direkt im Skript nach der DATALINES-Anweisung enthalten, was bedeutet, dass das Dataset aus literalen Daten erstellt wird.

Kopiert!

1	DATA tree;
2	INPUT tree day y;
3	DATALINES;
4	1 118 30
5	1 484 58
6	1 664 87
7	1 1004 115
8	1 1231 120
9	1 1372 142
10	1 1582 145
11	2 118 33
12	2 484 69
13	2 664 111
14	2 1004 156
15	2 1231 172
16	2 1372 203
17	2 1582 203
18	3 118 30
19	3 484 51
20	3 664 75
21	3 1004 108
22	3 1231 115
23	3 1372 139
24	3 1582 140
25	4 118 32
26	4 484 62
27	4 664 112
28	4 1004 167
29	4 1231 179
30	4 1372 209
31	4 1582 214
32	5 118 30
33	5 484 49
34	5 664 81
35	5 1004 125
36	5 1231 142
37	5 1372 174
38	5 1582 177
39	;
40

2 Codeblock

PROC NLMIXED

Erklärung :
Dieser Block führt die NLMIXED-Prozedur aus, um ein nichtlineares Modell mit gemischten Effekten an das 'tree'-Dataset anzupassen.
- Die PARMS-Anweisung initialisiert die Modellparameter: b1, b2, b3 (feste Parameter der logistischen Kurve), s2u (Varianz der Zufallseffekte) und s2e (Varianz des Restfehlers).
- Die Gleichungen 'num', 'ex', 'den' definieren die logistische Form des Modellmittelwerts.
- Die MODEL-Anweisung legt fest, dass die Antwortvariable 'y' einer Normalverteilung mit einem Mittelwert 'num/den' und einer Varianz 's2e' folgt.
- Die RANDOM-Anweisung deklariert 'u1' als einen Zufallseffekt, der einer Normalverteilung mit Mittelwert 0 und Varianz 's2u' folgt, wobei 'tree' die Subjektvariable ist, was darauf hinweist, dass die Zufallseffekte spezifisch für jeden Baum sind.

Kopiert!

1	PROC NLMIXED DATA=tree;
2	parms b1=190 b2=700 b3=350 s2u=1000 s2e=60;
3	num = b1+u1;
4	ex = exp(-(day-b2)/b3);
5	den = 1 + ex;
6	model y ~ normal(num/den,s2e);
7	random u1 ~ normal(0,s2u) subject=tree;
8	RUN;
9

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