Das Skript beginnt mit der Erstellung eines Datensatzes namens 'teratology' über einen DATA-Schritt und Inline-Daten (`cards`). Dieser Datensatz enthält Informationen über Würfe (litter), Behandlungsgruppen (group), die Gesamtzahl der Individuen (n) und die Anzahl der betroffenen Individuen (y). Indikatorvariablen (z2, z3, z4) werden für die Gruppen 2, 3 und 4 erstellt, um als Prädiktoren verwendet zu werden. Anschließend wird PROC GENMOD verwendet, um ein generalisiertes lineares Modell mit einer Binomialverteilung und einer Logit-Link-Funktion anzupassen. Die Option `repeated` mit `subject=litter` und `type=exch corrw` wird verwendet, um die Korrelation innerhalb des Subjekts (innerhalb des Wurfs) zu modellieren, was eine gängige Praxis in Teratologiestudien ist, um die Nicht-Unabhängigkeit der Beobachtungen innerhalb desselben Wurfs zu berücksichtigen.
Datenanalyse
Type : CREATION_INTERNE
Der Datensatz 'teratology' wird direkt im SAS-Skript unter Verwendung eines DATA-Schritts und der über die Anweisung `cards;` bereitgestellten Rohdaten erstellt.
1 Codeblock
DATA STEP Data
Erklärung : Dieser DATA-Schritt-Block erstellt den Datensatz 'teratology' aus Rohdaten, die direkt im Skript bereitgestellt werden. Die Variablen 'litter' (Wurf), 'group' (Behandlungsgruppe), 'n' (Gesamtzahl der Individuen) und 'y' (Anzahl der betroffenen Individuen) werden gelesen. Die binären Variablen 'z2', 'z3', 'z4' werden dann aus der Variablen 'group' abgeleitet, um die Behandlungsgruppen 2, 3 und 4 darzustellen. Diese Variablen dienen als Prädiktoren in der statistischen Analyse.
IF group=2THEN z2=1; IF group=3THEN z3=1; IF group=4THEN z4=1;
5
CARDS;
6
11101
7
21114
8
31129
9
4144
10
511010
11
61119
12
7199
13
811111
14
911010
15
101107
16
1111212
17
121109
18
13188
19
141119
20
15164
21
16197
22
1711414
23
181127
24
191119
25
201138
26
211145
27
2211010
28
2311210
29
241138
30
2511010
31
261143
32
2711313
33
28143
34
29188
35
301135
36
3111212
37
322101
38
33231
39
342131
40
35212 0
41
362144
42
37292
43
382132
44
392161
45
40211 0
46
4124 0
47
4221 0
48
43212 0
49
4438 0
50
453111
51
46314 0
52
473141
53
48311 0
54
4943 0
55
50413 0
56
51492
57
524172
58
53415 0
59
5442 0
60
554141
61
5648 0
62
5746 0
63
58417 0
64
;
2 Codeblock
PROC GENMOD
Erklärung : Dieser Block verwendet die PROC GENMOD, um ein generalisiertes lineares Modell an die Teratologiedaten anzupassen. Die Option `DATA=teratology` gibt den Eingabedatensatz an. `class litter;` deklariert 'litter' als Klassifikationsvariable, die für die Analyse von Wiederholungsdaten unerlässlich ist. Die Anweisung `model y/n= z2 z3 z4/dist=b link=logit;` spezifiziert, dass die Antwort 'y' unter 'n' Versuchen einer Binomialverteilung (`dist=b`) mit einer Logit-Link-Funktion (`link=logit`) folgt und dass die Prädiktoren 'z2', 'z3' und 'z4' sind. Die Option `repeated subject=litter/type=exch corrw;` ist entscheidend, um die Korrelation innerhalb des Subjekts oder innerhalb des Wurfs zu modellieren, wobei eine austauschbare Korrelationsstruktur (`type=exch`) angenommen und eine Arbeitskorrelationsmatrix (`corrw`) verwendet wird, was ein Standardansatz für diese Art von Studie ist.
Kopiert!
proc genmod DATA= teratology;
class litter;
model y/n= z2 z3 z4/dist=b link=logit;
repeated subject=litter/type=exch corrw;
run;
1
PROC GENMODDATA= teratology;
2
class litter;
3
model y/n= z2 z3 z4/dist=b link=logit;
4
repeated subject=litter/type=exch corrw;
5
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.
SAS und alle anderen Produkt- oder Dienstleistungsnamen von SAS Institute Inc. sind eingetragene Marken oder Marken von SAS Institute Inc. in den USA und anderen Ländern. ® zeigt die Registrierung in den USA an. WeAreCAS ist eine unabhängige Community-Site und nicht mit SAS Institute Inc. verbunden.
Diese Website verwendet technische und analytische Cookies, um Ihre Erfahrung zu verbessern.
Mehr erfahren.