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;
     logbili=log(bili);
     logalkph=log(alkph);
     logast=log(asptr);
     sex=1-sex;
     proc sort; BY tment;
RUN;    

PROC FREQ DATA=pbc3;
    TABLES sex*tment/NOROW NOPERCENT;
RUN;    

PROC UNIVARIATE DATA=pbc3;
    BY tment;
    VAR alb logbili logalkph logast age;
RUN;    

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

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

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

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

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=logalkph;
RUN;    

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

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=logast;
RUN;    

/* NB effect of tment differs from Table 6.2.17 */

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

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

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

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=sex;
RUN;    

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