options ls=80 nodate;
libname fu '../../sdata';
data work;
	set fu.prostat;
/* interactions */

tsize = treat*size;
tindex = treat * index;

proc phreg;
        model st*censor(0)= treat size index tsize tindex / rl;
run;