options ls=80;
libname fu '../../sdata';
data work;
        set fu.dialysis;
proc phreg;
        model infectt*censor(0)= age sex;
        output out=outp wtressch = wschage wschsex;
data fu.phdial;
        set outp;
        agex = wschage + 0.03;
        sexx = wschsex - 2.711;
filename gsasfile 'phdial.gsf';
goptions reset=all gunit=pct border ftext=swissb htitle=6 htext=2.5
gaccess=gsasfile ROTATE=LANDSCAPE gsfmode=append device=ps;
proc gplot;
        plot agex*infectt;
        plot sexx*infectt;
proc reg;
        model agex = infectt;
proc reg;
        model sexx = infectt;
run;