Análisis de supervivencia bayesiano con PROC LIFEREG (Ejemplo 2)

Este script ilustra un análisis de supervivencia sobre datos de fiabilidad de ventiladores industriales. Comienza creando el conjunto de datos 'Fan' que contiene los tiempos de funcionamiento y los indicadores de censura. La rutina LIFEREG se utiliza luego para ajustar un modelo paramétrico log-normal con un enfoque bayesiano (instrucción BAYES), generando muestras a posteriori en la tabla 'Post'. Un paso DATA posterior utiliza estas muestras para estimar la distribución de la probabilidad de fallo a 8000 horas. Finalmente, PROC MEANS proporciona estadísticas descriptivas (media y percentiles) sobre esta probabilidad estimada.

Análisis de datos

Type : CREACION_INTERNA

Los datos se definen directamente en el script a través de la instrucción 'datalines' en el paso DATA Fan.

1 Bloque de código

DATA STEP Data

Explicación :
Crea la tabla SAS 'Fan' leyendo datos internos (datalines). La variable 'Lifetime' representa la duración de la vida útil y 'Censor' indica si la observación está censurada (1) o no (0).

¡Copiado!

1	DATA Fan;
2	INPUT Lifetime Censor @code_sas_json/8_SAS_Intro_ReadFile_MultiCol_@@.json;
3	DATALINES;
4	450 0 460 1 1150 0 1150 0 1560 1
5	1600 0 1660 1 1850 1 1850 1 1850 1
6	1850 1 1850 1 2030 1 2030 1 2030 1
7	2070 0 2070 0 2080 0 2200 1 3000 1
8	3000 1 3000 1 3000 1 3100 0 3200 1
9	3450 0 3750 1 3750 1 4150 1 4150 1
10	4150 1 4150 1 4300 1 4300 1 4300 1
11	4300 1 4600 0 4850 1 4850 1 4850 1
12	4850 1 5000 1 5000 1 5000 1 6100 1
13	6100 0 6100 1 6100 1 6300 1 6450 1
14	6450 1 6700 1 7450 1 7800 1 7800 1
15	8100 1 8100 1 8200 1 8500 1 8500 1
16	8500 1 8750 1 8750 0 8750 1 9400 1
17	9900 1 10100 1 10100 1 10100 1 11500 1
18	;

2 Bloque de código

PROC LIFEREG

Explicación :
Ejecuta un análisis de regresión sobre los datos de supervivencia con una distribución log-normal. La instrucción 'bayes' activa el análisis bayesiano, fija una semilla aleatoria para la reproducibilidad (seed=1) y exporta las muestras de la distribución a posteriori a la tabla 'Post'.

¡Copiado!

1	ods graphics on;
2	PROC LIFEREG DATA=Fan;
3	model Lifetime*Censor( 1 )= / dist=lognormal;
4	bayes seed=1 outpost=Post;
5	RUN;

3 Bloque de código

DATA STEP Data

Explicación :
Utiliza los parámetros estimados (Intercept, Scale) almacenados en la tabla 'Post' para calcular la variable 'Frac'. Esta variable representa la probabilidad estimada de fallo a 8000 horas para cada iteración de la simulación bayesiana.

¡Copiado!

1	DATA Prob;
2	SET Post;
3	Frac = ProbNorm(( log(8000) - Intercept ) / Scale );
4	label Frac= 'Fraction Failing in 8000 Hours';
5	RUN;

4 Bloque de código

PROC MEANS

¡Copiado!

1
2	PROC MEANS
3	DATA = Prob(keep=Frac) n mean p10 p25 p50 p75 p90;
4	RUN;
5

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.

Retour à la liste