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 / rl;
proc phreg;
        model survt*censor(0)= sex / rl;
proc phreg;
        model survt*censor(0)= bun / rl;
proc phreg;
        model survt*censor(0)= ca / rl;
proc phreg;
        model survt*censor(0)= hb / rl;
proc phreg;
        model survt*censor(0)= pc / rl;
proc phreg;
        model survt*censor(0)= bj / rl;
proc phreg;
        model survt*censor(0)= bun hb / rl;
proc phreg;
        model survt*censor(0)= hb bj / rl;
proc phreg;
        model survt*censor(0)= bun bj / rl;
proc phreg;
        model survt*censor(0)= bun hb bj / rl;
proc phreg;
        model survt*censor(0)= age bun hb / rl;
proc phreg;
        model survt*censor(0)= sex bun hb / rl;
proc phreg;
        model survt*censor(0)=  bun  hb ca / rl;
proc phreg;
        model survt*censor(0)= bun hb pc / rl;
run;