Ce script génère des données bivariées simulées et applique la procédure KDE pour estimer les densités. Il explore l'impact du multiplicateur de bande passante (bwm), demande des statistiques descriptives et des percentiles personnalisés, et produit divers graphiques univariés et bivariés via ODS Graphics.
Analyse des données
Type : CREATION_INTERNE
Les données 'bivnormal' sont générées par un Data Step utilisant la fonction rannor pour créer 1000 observations d'une distribution normale bivariée.
1 Bloc de code
DATA STEP Data
Explication : Création d'un jeu de données 'bivnormal' contenant 1000 observations simulées. Les variables x et y sont construites à partir de tirages aléatoires normaux pour introduire une corrélation.
Copié !
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 Bloc de code
PROC KDE
Explication : Exécution de l'estimation de densité bivariée avec modification de la bande passante (bandwidth). La première exécution utilise un multiplicateur global (bwm=2), la seconde spécifie des multiplicateurs différents pour x et y. ODS Graphics est activé pour visualiser les résultats.
Copié !
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;
1
ods graphics on;
2
PROC KDEDATA=bivnormal;
3
bivar x y / bwm=2;
4
RUN;
5
6
PROC KDEDATA=bivnormal;
7
bivar (x y) (x (bwm=0.5) y (bwm=2));
8
RUN;
9
ods graphics off;
3 Bloc de code
PROC KDE
Explication : Demande de tableaux de sortie supplémentaires : statistiques bivariées, niveaux de densité, percentiles et statistiques univariées.
Copié !
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 Bloc de code
PROC KDE
Explication : Estimation bivariée avec spécification explicite des niveaux de contours (levels) et des percentiles à calculer.
Copié !
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 Bloc de code
PROC KDE
Explication : Génération de graphiques univariés complets : densité seule, histogramme seul, et superposition histogramme/densité pour x. Une superposition des densités de x et y est également demandée.
Copié !
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;
Ce matériel est fourni "tel quel" par We Are Cas. Il n'y a aucune garantie, expresse ou implicite, quant à la qualité marchande ou à l'adéquation à un usage particulier concernant le matériel ou le code contenu dans les présentes. We Are Cas n'est pas responsable des erreurs dans ce matériel tel qu'il existe maintenant ou existera, et We Are Cas ne fournit pas de support technique pour celui-ci.
Informations de Copyright : SAS SAMPLE LIBRARY, NAME: KDEX2, PRODUCT: SAS
SAS et tous les autres noms de produits ou de services de SAS Institute Inc. sont des marques déposées ou des marques de commerce de SAS Institute Inc. aux États-Unis et dans d'autres pays. ® indique un enregistrement aux États-Unis. WeAreCAS est un site communautaire indépendant et n'est pas affilié à SAS Institute Inc.
Ce site utilise des cookies techniques et analytiques pour améliorer votre expérience.
En savoir plus.