Tiempo de Recuperación Natural - Ejemplo PROC LIFEREG

Este script ilustra el uso del procedimiento LIFEREG para analizar datos de supervivencia con censura arbitraria (intervalos). Modela la relación entre el tiempo de recuperación, la edad y el sexo, incluyendo una interacción. Se realizan dos análisis para generar gráficos de probabilidad (probplots) basados en datos de referencia específicos (xrow1 y xrow2).

Análisis de datos

Type : CREATION_INTERNE

Todos los datos (tablas 'mice', 'xrow1', 'xrow2') se definen directamente en el script a través de bloques DATALINES.

1 Bloque de código

DATA STEP Data

Explicación :
Creación del conjunto de datos principal 'mice' que contiene las variables explicativas (sexo, edad) y los límites del intervalo de tiempo (time1, time2).

¡Copiado!

1	title 'Natural Recovery Time';
2	DATA mice;
3	INPUT sex age time1 time2;
4	DATALINES;
5	1 57 631 631
6	1 45 . 170
7	1 54 227 227
8	1 43 143 143
9	1 64 916 .
10	1 67 691 705
11	1 44 100 100
12	1 59 730 .
13	1 47 365 365
14	1 74 1916 1916
15	2 79 1326 .
16	2 75 837 837
17	2 84 1200 1235
18	2 54 . 365
19	2 74 1255 1255
20	2 71 1823 .
21	2 65 537 637
22	2 33 583 683
23	2 77 955 .
24	2 46 577 577
25	;

2 Bloque de código

DATA STEP Data

Explicación :
Creación de un conjunto de datos de referencia 'xrow1' para predicción o trazado (Sexo 1, Edad 50).

¡Copiado!

1	DATA xrow1;
2	INPUT sex age time1 time2;
3	DATALINES;
4	1 50 . .
5	;

3 Bloque de código

DATA STEP Data

Explicación :
Creación de un segundo conjunto de datos de referencia 'xrow2' (Sexo 2, Edad 60.6).

¡Copiado!

1	DATA xrow2;
2	INPUT sex age time1 time2;
3	DATALINES;
4	2 60.6 . .
5	;

4 Bloque de código

PROC LIFEREG

Explicación :
Ejecución del modelo de supervivencia de Weibull con los datos 'mice'. La opción 'xdata=xrow1' se utiliza para personalizar el gráfico de probabilidad. La instrucción 'model' define los límites de censura (time1, time2) y las covariables.

¡Copiado!

1	ods graphics on;
2	PROC LIFEREG DATA=mice xdata=xrow1;
3	class sex;
4	model (time1, time2) = age sex age*sex / dist=Weibull;
5
6	probplot / nodata
7	plower=.5
8	vref(intersect) = 75
9	vreflab = '75 Percent';
10	inset;
11	RUN;

5 Bloque de código

PROC LIFEREG

Explicación :
Segunda ejecución del modelo idéntica a la primera, pero utilizando 'xrow2' como datos de referencia para el gráfico.

¡Copiado!

1	PROC LIFEREG DATA=mice xdata=xrow2;
2	class sex;
3	model (time1, time2) = age sex age*sex / dist=Weibull;
4
5	probplot / nodata
6	plower=.5
7	vref(intersect) = 75
8	vreflab = '75 Percent';
9	inset;
10	RUN;
11	title;

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: LIFEREX4

Retour à la liste