El script comienza creando un conjunto de datos 'galaxies' que contiene las velocidades de varias galaxias. Luego, aplica el procedimiento HPFMM en tres pasos: 1) Búsqueda del número óptimo de componentes (de 3 a 7) con varianzas desiguales, basándose en el criterio AIC. 2) Misma búsqueda pero forzando varianzas iguales entre los componentes. 3) Ajuste de un modelo final de 5 componentes con una restricción sobre el valor de la varianza común.
Análisis de datos
Type : CREACIÓN_INTERNA
Los datos se crean directamente en el script a través de un paso DATA y una instrucción DATALINES. La variable 'velocity' se lee y se transforma en una nueva variable 'v' para el análisis.
1 Bloque de código
DATA STEP Data
Explicación : Este bloque DATA STEP lee los datos de velocidad de las galaxias proporcionados a través de 'datalines'. El operador '@@' (double trailing at) permite leer múltiples observaciones en la misma línea de datos. Se calcula una nueva variable 'v' dividiendo 'velocity' por 1000 para el escalado.
Explicación : Primer análisis con HPFMM para determinar el número óptimo de componentes (entre 3 y 7, opciones kmin y kmax) basado en el criterio de información de Akaike (AIC). Por defecto, las varianzas de los componentes normales se estiman por separado (desiguales). Los gráficos ODS están activados y algunas tablas de salida (historial de iteraciones, información de optimización) están ocultas.
¡Copiado!
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 Bloque de código
PROC HPFMM
Explicación : Segundo análisis con HPFMM, similar al primero, pero con la restricción de que las varianzas de los componentes sean iguales (opción EQUATE=SCALE). El criterio de convergencia sobre el gradiente está desactivado (gconv=0).
¡Copiado!
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 Bloque de código
PROC HPFMM
Explicación : Tercer y último análisis que ajusta un modelo específico de 5 componentes (K=5), con varianzas iguales (EQUATE=SCALE). La instrucción RESTRICT añade una restricción para fijar el valor de esta varianza común a 0.9025. Finalmente, los gráficos ODS están desactivados.
¡Copiado!
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;
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.
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.