Das Skript beginnt mit der Erstellung eines Datensatzes 'galaxies', der die Geschwindigkeiten mehrerer Galaxien enthält. Anschließend wird die HPFMM-Prozedur in drei Schritten angewendet: 1) Suche nach der optimalen Anzahl von Komponenten (3 bis 7) mit ungleichen Varianzen, basierend auf dem AIC-Kriterium. 2) Gleiche Suche, jedoch mit erzwungenen gleichen Varianzen zwischen den Komponenten. 3) Anpassung eines finalen Modells mit 5 Komponenten und einer Einschränkung des Wertes der gemeinsamen Varianz.
Datenanalyse
Type : CREATION_INTERNE
Die Daten werden direkt im Skript über einen DATA-Schritt und eine DATALINES-Anweisung erstellt. Die Variable 'velocity' wird gelesen und in eine neue Variable 'v' für die Analyse umgewandelt.
1 Codeblock
DATA STEP Data
Erklärung : Dieser DATA STEP-Block liest die über 'datalines' bereitgestellten Galaxien-Geschwindigkeitsdaten. Der Operator ' @@' (double trailing at) ermöglicht das Lesen mehrerer Beobachtungen in einer einzigen Datenzeile. Eine neue Variable 'v' wird durch Teilen von 'velocity' durch 1000 zur Skalierung berechnet.
Erklärung : Erste Analyse mit HPFMM zur Bestimmung der optimalen Anzahl von Komponenten (zwischen 3 und 7, Optionen kmin und kmax) basierend auf dem Akaike-Informationskriterium (AIC). Standardmäßig werden die Varianzen der normalen Komponenten separat (ungleich) geschätzt. ODS-Grafiken werden aktiviert und einige Ausgabetabellen (Iterationsverlauf, Optimierungsinformationen) werden ausgeblendet.
Kopiert!
title2 "Three to Seven Components, Unequal Variances";
ods graphics on;
proc hpfmm data=galaxies criterion=AIC;
model v = / kmin=3 kmax=7;
ods exclude IterHistory OptInfo ComponentInfo;
run;
1
title2 "Three to Seven Components, Unequal Variances";
2
ods graphics on;
3
PROC HPFMMDATA=galaxies criterion=AIC;
4
model v = / kmin=3 kmax=7;
5
ods exclude IterHistory OptInfo ComponentInfo;
6
RUN;
3 Codeblock
PROC HPFMM
Erklärung : Zweite Analyse mit HPFMM, ähnlich der ersten, jedoch mit der Einschränkung, dass die Varianzen der Komponenten gleich sind (Option EQUATE=SCALE). Das Konvergenzkriterium für den Gradienten wird deaktiviert (gconv=0).
Kopiert!
title2 "Three to Seven Components, Equal Variances";
proc hpfmm data=galaxies criterion=AIC gconv=0;
model v = / kmin=3 kmax=7 equate=scale;
run;
1
title2 "Three to Seven Components, Equal Variances";
2
PROC HPFMMDATA=galaxies criterion=AIC gconv=0;
3
model v = / kmin=3 kmax=7 equate=scale;
4
RUN;
4 Codeblock
PROC HPFMM
Erklärung : Dritte und letzte Analyse, die ein spezifisches Modell mit 5 Komponenten (K=5) mit gleichen Varianzen (EQUATE=SCALE) anpasst. Die RESTRICT-Anweisung fügt eine Einschränkung hinzu, um den Wert dieser gemeinsamen Varianz auf 0.9025 festzulegen. Schließlich werden die ODS-Grafiken deaktiviert.
Kopiert!
title2 "Five Components, Equal Variances = 0.9025";
proc hpfmm data=galaxies;
model v = / K=5 equate=scale;
restrict int 0 (scale 1) = 0.9025;
run;
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.
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.