data bmt1715; filename bmtdata url 'http://publicifsv.sund.ku.dk/~pka/epidata/bmt1715.txt'; infile bmtdata firstobs=2; input disease timedxtx sex karnofsky stage time donor event; /* disease 10=ALL 20=AML 30=CML timedxtx time from diagnosis to transplant in months sex 0=males 1=females karnofsky =1 if Karn>90, =0 if Karn<=90 stage 1=early 2=intermediate 3=advanced time from transplant to event/cens. in months donor 1=HLA-id sib 2=HLA-matched 3=HLA-mismatched event 0=cens. 1= relapse 2=death in remission */ run; /* a */ proc phreg data=bmt1715 plot(overlay=row)=cumhaz; model time*event(0 2)=; strata donor; run; proc phreg data=bmt1715 plot(overlay=row)=cumhaz; model time*event(0 1)=; strata donor; run; /* b */ proc phreg data=bmt1715; class donor (ref="1"); model time*event(0 2)= donor/rl; run; proc phreg data=bmt1715; class donor (ref="1"); model time*event(0 1)= donor/rl; run; /* c */ proc phreg data=bmt1715; class donor (ref="1") stage (ref="1") disease (ref="10"); model time*event(0 2)= donor disease stage karnofsky; run; proc phreg data=bmt1715; class donor (ref="1") stage (ref="1") disease (ref="10"); model time*event(0 1)= donor disease stage karnofsky; run; /* d */ proc phreg data=bmt1715 plot(overlay=row)=cif; model time*event(0)=/eventcode=1; strata donor; run; proc phreg data=bmt1715 plot(overlay=row)=cif; model time*event(0)=/eventcode=2; strata donor; run; /* e */ proc phreg data=bmt1715; class donor (ref="1"); model time*event(0)= donor/eventcode=1; run; proc phreg data=bmt1715; class donor (ref="1"); model time*event(0)= donor/eventcode=2; run; /* f */ proc phreg data=bmt1715; class donor (ref="1") stage (ref="1") disease (ref="10"); model time*event(0)= donor disease stage karnofsky/eventcode=1; run; proc phreg data=bmt1715; class donor (ref="1") stage (ref="1") disease (ref="10"); model time*event(0)= donor disease stage karnofsky/eventcode=2; run;