regression

logisticLackfit

Descripción

La acción logisticLackfit calcula la prueba de bondad de ajuste de Hosmer y Lemeshow para un modelo de regresión logística binaria. Esta prueba es útil para evaluar si las probabilidades predichas por el modelo coinciden con las probabilidades observadas en los datos. La acción requiere una tabla de entrada que contenga las variables del modelo y una tabla de restauración que contenga el modelo ajustado previamente por la acción 'logistic'.

regression.logisticLackfit / restore={...} table={...} <binEps=double> <cutpt=double-list> <df=double> <dfReduce=integer> <display={...}> <nGroups=integer> <noncentrality=double> <outputTables={...}> <powerAdj=boolean>;
Parámetros
ParámetroDescripción
restoreEspecifica la tabla de almacenamiento del modelo (ítem store) que contiene el modelo ajustado previamente. Este parámetro es obligatorio.
tableEspecifica la tabla de datos de entrada para la cual se evaluará la bondad de ajuste. Este parámetro es obligatorio.
binEpsEspecifica la precisión de las probabilidades predichas que se utilizan para la clasificación. El valor debe estar entre 0 y 1. El valor predeterminado es 1E-05.
cutptEspecifica los puntos de corte para las particiones de Hosmer y Lemeshow. Puede ser un único valor o una lista de valores.
dfEspecifica los grados de libertad a utilizar para la prueba de Hosmer y Lemeshow. El valor mínimo es 0.
dfReduceEspecifica la reducción en los grados de libertad para la prueba de Hosmer y Lemeshow. El valor predeterminado es 2.
displayEspecifica una lista de tablas de resultados para enviar al cliente para su visualización.
nGroupsEspecifica el número máximo de grupos a crear para la prueba de Hosmer y Lemeshow. El valor predeterminado es 10 y el mínimo es 5.
noncentralityEspecifica el parámetro de no centralidad para la prueba de Hosmer y Lemeshow. El valor predeterminado es 0.
outputTablesEnumera los nombres de las tablas de resultados para guardar como tablas CAS en el servidor.
powerAdjCuando se establece en True, ajusta el número de grupos para que la prueba de Hosmer y Lemeshow pueda mantener la potencia. El valor predeterminado es False.
Nota sobre la creación de datos

La documentación proporcionada no incluye un ejemplo explícito de creación de datos. La acción 'logisticLackfit' normalmente opera sobre una tabla de datos existente y un ítem store creado por la acción 'logistic'.

¡Copiado!
1/* No se proporcionó código de creación de datos en la documentación de origen. */

Ejemplos

Este ejemplo muestra cómo realizar una prueba de bondad de ajuste de Hosmer y Lemeshow básica utilizando un modelo previamente ajustado y almacenado en 'myModelStore' y datos de la tabla 'myTable'.

Código SAS® / CAS Código en espera de validación por la comunidad
¡Copiado!
1PROC CAS;
2regression.logisticLackfit /
3 restore={name='myModelStore'},
4 TABLE={name='myTable'};
5RUN;
Resultado :
Se espera como salida una tabla 'LackfitChiSq' que contiene el estadístico chi-cuadrado de Hosmer y Lemeshow, los grados de libertad y el p-valor asociado.

Este ejemplo realiza la prueba de Hosmer y Lemeshow especificando un número diferente de grupos (nGroups=8) y ajustando los grados de libertad (dfReduce=1). También genera la tabla de particiones como una tabla de salida CAS llamada 'MyPartitions'.

Código SAS® / CAS Código en espera de validación por la comunidad
¡Copiado!
1PROC CAS;
2regression.logisticLackfit /
3 restore={name='myModelStore'},
4 TABLE={name='myTable'},
5 nGroups=8,
6 dfReduce=1,
7 outputTables={names={LackfitPartitions='MyPartitions'}}$
8RUN;
Resultado :
Además de la tabla 'LackfitChiSq', se creará una nueva tabla CAS llamada 'MyPartitions' en la caslib activa, que contiene los detalles de las particiones utilizadas para la prueba (observados, esperados, etc.).

En este ejemplo, en lugar de dejar que la acción determine los grupos, se especifican explícitamente los puntos de corte de probabilidad para formar las particiones de la prueba.

Código SAS® / CAS Código en espera de validación por la comunidad
¡Copiado!
1PROC CAS;
2regression.logisticLackfit /
3 restore={name='myModelStore'},
4 TABLE={name='myTable'},
5 cutpt={0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9};
6RUN;
Resultado :
La salida será la tabla 'LackfitChiSq', pero los cálculos se basarán en los 10 grupos definidos por los 9 puntos de corte proporcionados, en lugar de los grupos basados en cuantiles que se usan por defecto.

FAQ

¿Cuál es el propósito de la acción `logisticLackfit`?
¿Qué es la prueba de Hosmer y Lemeshow?
¿Cuáles son los parámetros clave para usar la acción `logisticLackfit`?
¿Cómo se especifica el modelo que se va a evaluar?
¿Qué determina el parámetro `nGroups`?