options ls=80 nodate; libname fu '../../sdata'; data work; set fu.myeloma; /* in Table 1.3, 1 for male and 2 for female */ if sex=1 then sex=0; else sex=1; proc phreg; model survt*censor(0)= age sex bun ca hb pc bj /covb rl; run;