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);
RUN;

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili/RISKLIMITS;
RUN;

/* NB: slightly different Wald statistic due to rounding in book */

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

PROC MEANS DATA=pbc3 MEAN NWAY;
     CLASS biligroup;
     VAR bili logbili;
     OUTPUT OUT=means MEAN=bilimean logbilimean;
RUN;

PROC SORT DATA=means; BY biligroup; RUN;

PROC SORT DATA=pbc3; BY biligroup; RUN;    

DATA pbc3; MERGE pbc3 means;
     BY biligroup;
RUN;     

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bilimean/RISKLIMITS;
RUN;

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=logbilimean/RISKLIMITS;
RUN;