El script comienza con la creación de un conjunto de datos interno llamado 'Males'. Este dataset se alimenta de datos brutos (datalines) que describen el número de fallos y retiros por período de tiempo. Los datos se transforman para crear observaciones individuales con un indicador de censura (0 para eventos, 1 para casos censurados) y una variable de frecuencia. Luego, se utiliza el procedimiento PROC LIFETEST para analizar estos datos de supervivencia. Aplica el método de la tabla de vida ('method=lt'), define intervalos de tiempo específicos para el análisis y genera una serie de gráficos (funciones de supervivencia, log-supervivencia, log-log supervivencia, riesgo acumulado, densidad de probabilidad) a través de la opción 'plots'.
Análisis de datos
Type : CREATION_INTERNE
El conjunto de datos 'Males' se crea y se rellena directamente en el script utilizando un bloque DATA STEP y datos brutos especificados mediante la instrucción DATALINES. Las variables 'Freq', 'Years' y 'Censored' se construyen a partir de estos datos brutos.
1 Bloque de código
DATA STEP Data
Explicación : Este bloque DATA STEP crea y rellena el conjunto de datos 'Males'. Lee el número de fallos ('fail') y retiros ('withdraw') de las líneas de datos proporcionadas (datalines). Para cada par 'fail'/'withdraw', se generan dos observaciones: una para los eventos (Censored=0) con 'Freq' igual a 'fail', y otra para los casos censurados (Censored=1) con 'Freq' igual a 'withdraw'. La variable 'Years' se incrementa para definir los períodos de seguimiento. La variable 'Years' se inicializa en -0.5 para que la siguiente incrementación comience en 0.5.
¡Copiado!
title 'Survival of Males with Angina Pectoris';
data Males;
keep Freq Years Censored;
retain Years -.5;
input fail withdraw;
Years + 1;
Censored=0;
Freq=fail;
output;
Censored=1;
Freq=withdraw;
output;
datalines;
456 0 226 39 152 22 171 23 135 24 125 107
83 133 74 102 51 68 42 64 43 45 34 53
18 33 9 27 6 23 0 30
;
1
title 'Survival of Males with Angina Pectoris';
2
DATA Males;
3
keep Freq Years Censored;
4
retain Years -.5;
5
INPUT fail withdraw;
6
Years + 1;
7
Censored=0;
8
Freq=fail;
9
OUTPUT;
10
Censored=1;
11
Freq=withdraw;
12
OUTPUT;
13
DATALINES;
14
456 0 22639152221712313524125107
15
83133741025168426443453453
16
1833927623 0 30
17
;
2 Bloque de código
PROC LIFETEST
¡Copiado!
ods graphics on;
proc lifetest data=Males method=lt intervals=(0 to 15 by 1)
plots=(s,ls,lls,h,p);
time Years*Censored(1);
freq Freq;
run;
ods graphics off;
1
ods graphics on;
2
PROC LIFETESTDATA=Males method=lt intervals=(0 to 15BY1)
3
plots=(s,ls,lls,h,p);
4
time Years*Censored(1);
5
freq Freq;
6
RUN;
7
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 (NAME: LIFTEX3, PRODUCT: STAT)
SAS y todos los demás nombres de productos o servicios de SAS Institute Inc. son marcas registradas o marcas comerciales de SAS Institute Inc. en los EE. UU. y otros países. ® indica registro en los EE. UU. WeAreCAS es un sitio comunitario independiente y no está afiliado a SAS Institute Inc.
Este sitio utiliza cookies técnicas y analíticas para mejorar su experiencia.
Saber más.