FILENAME goptions URL "http://192.38.117.59/~linearpredictors/datafiles/goptions.sas";
%include goptions;

FILENAME url URL "http://192.38.117.59/~linearpredictors/datafiles/readPbc3.sas";
%include url;

DATA pbc3; SET pbc3;
     stage2=(stage=3 OR stage=4); if stage=. then stage2=.;
     logbili=log(bili);
     logalkph=log(alkph);
     logast=log(asptr);
     sex=1-sex;
     proc sort; BY tment;
RUN;    

/* Model 1 */

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=tment alb logbili logalkph logast age;
RUN;    

/* Model 2 - NB: SD for logast differs from that of Table 6.2.18 */

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=tment alb logbili logalkph logast;
RUN;    

/* Model 3 */

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=tment alb logbili age;
RUN;    

/* Model 4 */

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=tment alb logbili;
RUN;    

/* Model 5 */

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=tment alb logbili logalkph logast age stage2;
RUN;