Ejemplo de documentación 2 para PROC LIFETEST

El script utiliza el conjunto de datos 'Sashelp.BMT' (Trasplante de Médula Ósea) para realizar análisis de supervivencia. Demuestra cómo mostrar la tabla de sujetos en riesgo debajo del gráfico de supervivencia, realizar pruebas de Log-rank con ajuste de Sidak para comparaciones múltiples, y comparar grupos con un grupo de control específico. También muestra cómo usar PROC FORMAT y un paso DATA para reclasificar los estratos en un orden específico (no alfabético) para los gráficos, y cómo añadir bandas de confianza de Hall-Wellner.

Análisis de datos

Type : SASHELP

El script utiliza la tabla estándar 'Sashelp.BMT'. Se crea una tabla derivada temporal 'Bmt2' para ilustrar la manipulación del orden de los estratos mediante formatos.

1 Bloque de código

PROC PRINT

Explicación :
Muestra las 10 primeras observaciones de la tabla Sashelp.BMT para una vista previa de los datos.

¡Copiado!

1	PROC PRINT DATA=Sashelp.BMT(obs=10);
2	RUN;

2 Bloque de código

PROC LIFETEST

Explicación :
Activa ODS Graphics. Ejecuta un análisis de supervivencia de Kaplan-Meier. La opción 'plots=survival(atrisk=...)' muestra el número de sujetos en riesgo en intervalos específicos (cada 500 días) debajo del gráfico. La opción 'adjust=sidak' aplica un ajuste de Sidak a los valores p de la prueba Log-rank para comparaciones múltiples entre los grupos.

¡Copiado!

1	ods graphics on;
2
3	PROC LIFETEST DATA=sashelp.BMT plots=survival(atrisk=0 to 2500 BY 500);
4	time T * STATUS(0);
5	strata Group / test=logrank adjust=sidak;
6	RUN;

3 Bloque de código

PROC LIFETEST

Explicación :
Realiza un análisis comparativo sin producir gráficos ('plots=none'). La opción 'diff=control' compara todos los demás grupos con el grupo de control especificado ('AML-Low Risk') con ajuste de Sidak.

¡Copiado!

1	PROC LIFETEST DATA=sashelp.BMT notable plots=none;
2	time T * STATUS(0);
3	strata Group / test=logrank adjust=sidak diff=control('AML-Low Risk');
4	RUN;

4 Bloque de código

PROC FORMAT

¡Copiado!

1	PROC FORMAT;
2	invalue $bmtifmt 'ALL' = 1 'AML-Low Risk' = 2 'AML-High Risk' = 3;
3	value bmtfmt 1 = 'ALL' 2 = 'AML-Low Risk' 3 = 'AML-High Risk';
4	RUN;

5 Bloque de código

DATA STEP Data

¡Copiado!

1	DATA Bmt2;
2	SET sashelp.BMT(rename=(Group=G));
3	Group = INPUT(INPUT(G, $bmtifmt.), 1.);
4	label Group = 'Disease Group';
5	FORMAT Group bmtfmt.;
6	RUN;

6 Bloque de código

PROC LIFETEST

¡Copiado!

1	PROC LIFETEST DATA=Bmt2 plots=s(atrisk(outside maxlen=13)=0 to 2500 BY 500);
2	time T*STATUS(0);
3	strata Group / order=internal;
4	RUN;

7 Bloque de código

PROC LIFETEST

¡Copiado!

1	PROC LIFETEST DATA=Bmt2 plots=survival(cl cb=hw strata=panel);
2	time T * STATUS(0);
3	strata Group/order=internal;
4	RUN;
5
6	ods 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 : S A S S A M P L E L I B R A R Y

Retour à la liste