Regresión de Poisson exacta con PROC GENMOD

El script comienza creando un conjunto de datos llamado 'ingots' que contiene el número de lingotes defectuosos ('Notready') para diferentes combinaciones de tiempo de calentamiento ('Heat') y tiempo de remojo ('Soak'), así como el total probado ('Total'). Luego calcula el logaritmo del total ('lnTotal') para usarlo como variable de compensación en el modelo. Posteriormente se llama al procedimiento PROC GENMOD para ajustar un modelo de regresión de Poisson, donde 'Notready' es la variable dependiente. Se solicita un análisis condicional exacto para los parámetros 'Heat' y 'Soak' para obtener estimaciones y pruebas más precisas, lo que es particularmente útil con muestras pequeñas.

Análisis de datos

Type : CREATION_INTERNE

Los datos se crean directamente en el script utilizando un paso DATA y la instrucción 'datalines'.

1 Bloque de código

DATA STEP Data

Explicación :
Este bloque de datos crea la tabla 'ingots'. Lee las variables 'Heat', 'Soak', 'Notready' y 'Total' de los datos en línea (datalines), luego calcula una nueva variable 'lnTotal' que es el logaritmo de 'Total'. Esta tabla es la fuente para el análisis estadístico.

¡Copiado!

1	DATA ingots;
2	INPUT Heat Soak Notready Total @code_sas_json/8_SAS_Intro_ReadFile_MultiCol_@@.json;
3	lnTotal= log(Total);
4	DATALINES;
5	7 1.0 0 10 14 1.0 0 31 27 1.0 1 56 51 1.0 3 13
6	7 1.7 0 17 14 1.7 0 43 27 1.7 4 44 51 1.7 0 1
7	7 2.2 0 7 14 2.2 2 33 27 2.2 0 21 51 2.2 0 1
8	7 2.8 0 12 14 2.8 0 31 27 2.8 1 22 51 4.0 0 1
9	7 4.0 0 9 14 4.0 0 19 27 4.0 1 16
10	;
11

2 Bloque de código

PROC GENMOD

Explicación :
Este bloque utiliza PROC GENMOD para realizar una regresión de Poisson. 'Heat' y 'Soak' se definen como variables de clase. El modelo predice 'Notready' en función de 'Heat' y 'Soak', con 'lnTotal' como variable de compensación para normalizar la exposición. La instrucción 'exact' solicita un análisis condicional exacto para las variables 'Heat' y 'Soak', y 'exactoptions' define un tiempo límite para el estado de este análisis.

¡Copiado!

1	PROC GENMOD DATA=ingots;
2	class Heat Soak / param=ref;
3	model Notready=Heat Soak / offset=lnTotal dist=Poisson link=log;
4	exact Heat Soak / joint estimate;
5	exactoptions statustime=10;
6	RUN;

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: GENMEX11, TITLE: Example 11 for PROC GENMOD, PRODUCT: STAT

Retour à la liste