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;

/* fit the basic Cox model with linear effect of bilirubin */
    
PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili;
RUN;

/* create variables for fractional polynomials of orders
   2, 3, 0.5, 0, -0.5, -1, -2, -3 */ 

DATA pbc3; SET pbc3;
    logbili=log(bili);
    bili2=bili*bili;
    bili3=bili2*bili;
    bili05=sqrt(bili);
    bili_05=1/bili05;
    bili_1=1/bili;
    bili_2=1/bili2;
    bili_3=1/bili3;
RUN;    

/* fit all the two-covariate Cox models from table */

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

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

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

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

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

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

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

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

/* fit all the three-covariate Cox models from table */

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili_3 bili_2;
RUN;

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili_3 bili_1;
RUN;

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili_3 bili_05;
RUN;

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

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili_3 bili05;
RUN;

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili_3 bili2;
RUN;

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili_3 bili3;
RUN;

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili_2 bili_1;
RUN;

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili_2 bili_05;
RUN;

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

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili_2 bili05;
RUN;

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili_2 bili2;
RUN;

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili_2 bili3;
RUN;

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili_1 bili_05;
RUN;

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

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili_1 bili05;
RUN;

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili_1 bili2;
RUN;

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili_1 bili3;
RUN;

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

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili_05 bili05;
RUN;

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili_05 bili2;
RUN;

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili_05 bili3;
RUN;

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

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

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

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili05 bili2;
RUN;

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili05 bili3;
RUN;

PROC PHREG DATA=pbc3;
    MODEL followup*status(0)=bili bili2 bili3;
RUN;

/* the LR tests from the table may now be obtained by subtraction */