Dieses Skript generiert simulierte bivariate Daten und wendet die KDE-Prozedur an, um Dichten zu schätzen. Es untersucht den Einfluss des Bandbreitenmultiplikators (bwm), fordert deskriptive Statistiken und benutzerdefinierte Perzentile an und erstellt verschiedene univariate und bivariate Grafiken über ODS Graphics.
Datenanalyse
Type : CREATION_INTERNE
Die Daten 'bivnormal' werden durch einen DATA Step generiert, der die Funktion rannor verwendet, um 1000 Beobachtungen einer bivariaten Normalverteilung zu erstellen.
1 Codeblock
DATA STEP Data
Erklärung : Erstellung eines 'bivnormal' Datensatzes mit 1000 simulierten Beobachtungen. Die Variablen x und y werden aus normalverteilten Zufallszahlen konstruiert, um eine Korrelation einzuführen.
Kopiert!
data bivnormal;
seed = 1283470;
do i = 1 to 1000;
z1 = rannor(seed);
z2 = rannor(seed);
z3 = rannor(seed);
x = 3*z1+z2;
y = 3*z1+z3;
output;
end;
drop seed;
run;
1
DATA bivnormal;
2
seed = 1283470;
3
DO i = 1 to 1000;
4
z1 = rannor(seed);
5
z2 = rannor(seed);
6
z3 = rannor(seed);
7
x = 3*z1+z2;
8
y = 3*z1+z3;
9
OUTPUT;
10
END;
11
drop seed;
12
RUN;
2 Codeblock
PROC KDE
Erklärung : Ausführung der bivariaten Dichteschätzung mit Bandbreitenmodifikation (bandwidth). Die erste Ausführung verwendet einen globalen Multiplikator (bwm=2), die zweite gibt unterschiedliche Multiplikatoren für x und y an. ODS Graphics ist zur Visualisierung der Ergebnisse aktiviert.
Kopiert!
ods graphics on;
proc kde data=bivnormal;
bivar x y / bwm=2;
run;
proc kde data=bivnormal;
bivar (x y) (x (bwm=0.5) y (bwm=2));
run;
ods graphics off;
proc kde data=bivnormal;
bivar x y / bivstats levels percentiles unistats;
run;
1
2
PROC KDE
3
DATA=bivnormal;
4
bivar x y / bivstats levels percentiles unistats;
5
RUN;
6
4 Codeblock
PROC KDE
Erklärung : Bivariate Schätzung mit expliziter Angabe der Konturniveaus (levels) und der zu berechnenden Perzentile.
Kopiert!
proc kde data=bivnormal;
bivar x y / levels=2.5, 50, 97.5
percentiles=2.5, 25, 50, 75, 97.5;
run;
1
PROC KDEDATA=bivnormal;
2
bivar x y / levels=2.5, 50, 97.5
3
percentiles=2.5, 25, 50, 75, 97.5;
4
RUN;
5 Codeblock
PROC KDE
Erklärung : Generierung vollständiger univariater Grafiken: Dichte allein, Histogramm allein und Histogramm/Dichte-Überlagerung für x. Eine Überlagerung der Dichten von x und y wird ebenfalls angefordert.
Kopiert!
ods graphics on;
proc kde data=bivnormal;
univar x / plots=(density histogram histdensity);
univar x y / plots=densityoverlay;
run;
ods graphics off;
1
ods graphics on;
2
PROC KDEDATA=bivnormal;
3
univar x / plots=(density histogram histdensity);
4
univar x y / plots=densityoverlay;
5
RUN;
6
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 : SAS SAMPLE LIBRARY, NAME: KDEX2, PRODUCT: SAS
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.