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;

PROC UNIVARIATE DATA=pbc3;
    VAR alb bili alkph asptr age;
RUN;    

DATA pbc3; SET pbc3;
     alb30=(alb>30)*(alb-30); alb40=(alb>40)*(alb-40);
     logbili=log(bili);
     logbili10=(bili>10)*(logbili-log(10));
     logbili40=(bili>40)*(logbili-log(40));
     logalkph=log(alkph);
     logalkph500=(alkph>500)*(logalkph-log(500));
     logalkph1200=(alkph>1200)*(logalkph-log(1200));
     logast=log(asptr);
     logast60=(asptr>60)*(logast-log(60));
     logast120=(asptr>120)*(logast-log(120));
     age50=(age>50)*(age-50); age60=(age>60)*(age-60);
     sex=1-sex;
RUN;    

/* MODEL 1 */

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

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

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

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

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