Este script genera datos bivariados simulados y aplica el procedimiento KDE para estimar las densidades. Explora el impacto del multiplicador de ancho de banda (bwm), solicita estadísticas descriptivas y percentiles personalizados, y produce varios gráficos univariados y bivariados a través de ODS Graphics.
Análisis de datos
Type : CREATION_INTERNE
Los datos 'bivnormal' son generados por un Data Step utilizando la función rannor para crear 1000 observaciones de una distribución normal bivariada.
1 Bloque de código
DATA STEP Data
Explicación : Creación de un conjunto de datos 'bivnormal' que contiene 1000 observaciones simuladas. Las variables x e y se construyen a partir de extracciones aleatorias normales para introducir una correlación.
¡Copiado!
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 Bloque de código
PROC KDE
Explicación : Ejecución de la estimación de densidad bivariada con modificación del ancho de banda. La primera ejecución utiliza un multiplicador global (bwm=2), la segunda especifica multiplicadores diferentes para x e y. ODS Graphics está activado para visualizar los resultados.
¡Copiado!
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 Bloque de código
PROC KDE
Explicación : Solicitud de tablas de salida adicionales: estadísticas bivariadas, niveles de densidad, percentiles y estadísticas univariadas.
¡Copiado!
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 Bloque de código
PROC KDE
Explicación : Estimación bivariada con especificación explícita de los niveles de contorno (levels) y los percentiles a calcular.
¡Copiado!
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 Bloque de código
PROC KDE
Explicación : Generación de gráficos univariados completos: densidad sola, histograma solo, y superposición histograma/densidad para x. También se solicita una superposición de las densidades de x e y.
¡Copiado!
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;
Este material se proporciona "tal cual" por We Are Cas. No hay garantías, expresas o implícitas, en cuanto a la comerciabilidad o idoneidad para un propósito particular con respecto a los materiales o el código contenidos en este documento. We Are Cas no es responsable de los errores en este material tal como existe ahora o existirá, ni We Are Cas proporciona soporte técnico para el mismo.
Información de copyright : SAS SAMPLE LIBRARY, NAME: KDEX2, PRODUCT: SAS
SAS y todos los demás nombres de productos o servicios de SAS Institute Inc. son marcas registradas o marcas comerciales de SAS Institute Inc. en los EE. UU. y otros países. ® indica registro en los EE. UU. WeAreCAS es un sitio comunitario independiente y no está afiliado a SAS Institute Inc.
Este sitio utiliza cookies técnicas y analíticas para mejorar su experiencia.
Saber más.