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;

/* create new variables relating to fractional polynomials and normal range */

DATA pbc3; SET pbc3;
   logbili=log(bili);
   bili_05=1/sqrt(bili);
   bili_2=1/(bili*bili);
   lbilicent=logbili-log(17.1);
   bililarge=(bili>=17.1);
   logbilipos=bililarge*lbilicent;
RUN;

/* fit models from table */

PROC PHREG;
model followup*status(0)=tment;
RUN;

PROC PHREG;
model followup*status(0)=tment logbili;
RUN;

PROC PHREG;
model followup*status(0)=tment bili_05 bili_2;
RUN;

PROC PHREG;
model followup*status(0)=tment bililarge logbilipos;
RUN;