ANALYSIS
SURVIVAL HANDOUT #23
DATA HMOHIV;
INFILE 'A:\DATA4.TXT';
INPUT ID @13 ENTDATE DATE7. @23 ENDDATE DATE7. TIME AGE DRUG CENSOR;
MONTH=(ENDDATE-ENTDATE)/30;
DATA AA; INPUT DRUG @@; CARDS;
0 1
;
PROC LIFETEST DATA=HMOHIV OUTS=KM NOTABLE; STRATA DRUG;
TIME MONTH*CENSOR(0);
DATA KM; SET KM;
IF DRUG=0 THEN SURV0=SURVIVAL; IF DRUG=1 THEN SURV1=SURVIVAL;
PROC PHREG DATA=HMOHIV;
MODEL MONTH*CENSOR(0)=DRUG;
BASELINE OUT=BB COVARIATES=AA SURVIVAL=SURVIVAL;
DATA CC; SET BB; IF DRUG=0 OR DRUG=1;
IF DRUG=0 THEN SUR0=SURVIVAL; IF DRUG=1 THEN SUR1=SURVIVAL;
DATA FINAL; SET KM CC;
DATA EST;
INPUT DRUG;
CONTROL=1;
MONTH=0; CENSOR=0;
CARDS;
0
1
;
DATA SET LAST; HMOHIV EST;
PROC LIFEREG DATA=LAST;
MODEL MONTH*CENSOR(0)=DRUG/D=WEIBULL;
OUTPUT OUT=WEIB QUANTILES=.01 .05 .10 .15 .20 .25 .30 .35 .40 .45 .50
.55 .60 .65 .70 .75 .80 .85 .90 .95 .99
P=PRED CONTROL=CONTROL ;
DATA WEIB; SET WEIB; SURV=1-_PROB_; MONTH=PRED;
IF DRUG=0 THEN SURVV0=SURV;
IF DRUG=1 THEN SURVV1=SURV;
DATA FINAL; SET FINAL WEIB;
PROC PRINT;
RUN;
PROC GPLOT...