macro barlowpro t c nr tt ttt ttpl xx ww mconstant k1 k2 k3 k4 mcolumn t c nr tt ttt ttpl xx ww let k1=count(t) let nr(1)=k1 let tt(1)=k1*t(1) let ttt(1)=k1*t(1) let k3=c(1) if k3=1 let ttpl(1)=ttt(1) endif do k2=2:k1 let k3=k3+c(k2) let nr(k2)=k1+1-k2 let tt(k2)=nr(k2)*(t(k2)-t(k2-1)) let ttt(k2)=ttt(k2-1)+tt(k2) if c(k2)=1 let ttpl(k3)=ttt(k2) endif enddo do k4=1:k3 let xx(k4)=k4/k3 let ww(k4)=ttpl(k4)/ttpl(k3) enddo name t= 'Time' name c= 'Censor' name nr= 'No at risk' name tt = 'Total time' name ttt= 'Cum total time' name ttpl= 'Plot total time' name xx= 'i/n' name ww= 'TTT' plot ww*xx; title "TTT-plot censored data"; minimum 1 0; minimum 2 0. endmacro