macro nelson u v w x y zn zkm mconstant k1 k2 mcolumn u v w x y zkm zn Ltest u; Noparametric; KMEstimates; Survprobability zkm; Htime w; Hazard x; Censor v; Cvalue 0. let k1=count(x) let y(1)=x(1) do k2=2:k1 let y(k2)=y(k2-1)+x(k2) enddo name w 'Time' name x 'Haz Nelson' name y 'Cum Haz Nelson' name zn 'Survival Nelson' name zkm 'Survival KM' let zn = exp(-y) print w,y,zn,zkm plot zn*w; scale 1; min 0; scale 2; min 0; title "Nelson Survival Function". plot zkm*w; scale 1; min 0; scale 2; min 0; title "KM Survival Function". plot y*w; scale 1; min 0; scale 2; min 0; title "Nelson Plot". endmacro