Publicado el :
Estadística CREATION_INTERNE

Ejemplo introductorio de PROC LOESS: Incidencia de melanomas

Este código también está disponible en: Deutsch English Français
En espera de validación
Este script crea un conjunto de datos 'Melanoma' que contiene registros anuales. Comienza visualizando los datos brutos con PROC SGPLOT. Luego, utiliza PROC LOESS de forma progresiva: primero con los parámetros predeterminados, luego solicitando detalles estadísticos, comparando diferentes parámetros de suavizado (opción smooth), exportando los resultados a través de ODS OUTPUT y, finalmente, mostrando intervalos de confianza (CLM).
Análisis de datos

Type : CREATION_INTERNE


Los datos se incluyen directamente en el script a través de la instrucción DATALINES (conjunto de datos 'Melanoma').

1 Bloque de código
DATA STEP Data
Explicación :
Creación del conjunto de datos 'Melanoma'. La instrucción INPUT está configurada para leer múltiples pares 'Año Incidencia' por línea (nota: el código fuente contiene un artefacto de referencia de archivo donde se espera el símbolo ' @@' para la lectura continua).
¡Copiado!
1DATA Melanoma;
2 INPUT Year Incidences @code_sas_json/8_SAS_Intro_ReadFile_MultiCol_@@.json;
3 FORMAT Year d4.0;
4 DATALINES;
51936 0.9 1937 0.8 1938 0.8 1939 1.3
61940 1.4 1941 1.2 1942 1.7 1943 1.8
71944 1.6 1945 1.5 1946 1.5 1947 2.0
81948 2.5 1949 2.7 1950 2.9 1951 2.5
91952 3.1 1953 2.4 1954 2.2 1955 2.9
101956 2.5 1957 2.6 1958 3.2 1959 3.8
111960 4.2 1961 3.9 1962 3.7 1963 3.3
121964 3.7 1965 3.9 1966 4.1 1967 3.8
131968 4.7 1969 4.4 1970 4.8 1971 4.8
141972 4.8
15;
2 Bloque de código
PROC SGPLOT
Explicación :
Generación de un diagrama de dispersión (scatter plot) para visualizar la evolución de la incidencia en función del año antes de la modelización.
¡Copiado!
1 
2PROC SGPLOT
3DATA=Melanoma;
4scatter y=Incidences x=Year;
5RUN;
6 
3 Bloque de código
PROC LOESS
Explicación :
Activación de gráficos ODS. Primera ejecución de PROC LOESS con las opciones predeterminadas para ajustar una curva de regresión local.
¡Copiado!
1ods graphics on;
2 
3PROC LOESS DATA=Melanoma;
4 model Incidences=Year;
5RUN;
4 Bloque de código
PROC LOESS
Explicación :
Ejecución de PROC LOESS solicitando detalles adicionales en la salida: resumen del modelo y estadísticas de salida.
¡Copiado!
1 
2PROC LOESS
3DATA=Melanoma;
4model Incidences=Year / details(ModelSummary OutputStatistics);
5RUN;
6 
5 Bloque de código
PROC LOESS Data
Explicación :
Modelización con especificación explícita de varios parámetros de suavizado (smooth) para comparación. La opción 'residual' solicita el cálculo de los residuos. Las estadísticas de salida se guardan en la tabla SAS 'Results' a través de ODS OUTPUT.
¡Copiado!
1PROC LOESS DATA=Melanoma;
2 model Incidences=Year/smooth=0.1 0.25 0.4 0.6 residual;
3 ods OUTPUT OutputStatistics=Results;
4RUN;
6 Bloque de código
PROC PRINT
Explicación :
Visualización de las 5 primeras observaciones de la tabla 'Results' creada anteriormente.
¡Copiado!
1PROC PRINT DATA=Results(obs=5);
2 id obs;
3RUN;
7 Bloque de código
PROC LOESS
Explicación :
Ejecución destinada a producir un gráfico específico de los residuos por parámetro de suavizado (ResidualsBySmooth).
¡Copiado!
1 
2PROC LOESS
3DATA=Melanoma plots=ResidualsBySmooth(smooth);
4model Incidences=Year/smooth=0.1 0.25 0.4 0.6;
5RUN;
6 
8 Bloque de código
PROC LOESS
Explicación :
Último modelo solicitando los límites de confianza (CLM) con un umbral alfa de 0.1 (intervalo de confianza del 90%). Desactivación de los gráficos ODS al final.
¡Copiado!
1PROC LOESS DATA=Melanoma;
2 model Incidences=Year/clm alpha=0.1;
3RUN;
4 
5ods 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