Kaplan Meier curve

[Nishan Lamichhane]

dear all,

Is it possible to make a small colored diamond for each event /failure in a Kaplan Meier curve on Stata 17?

stset weeks, fail(med_stop)
sts graph if diagnosis==1, risktable ytitle (Persistent rate) xtitle (Weeks since first med) title ("Fig 1A. Kaplan-Meier curve for XYZ", size(8pt) color(black))

[Andrew Musau]

If you have varying number of failures at different points in time, then it may be more informative to report the number of failures rather than just indicating failures. For some studies, e.g., dieting, you expect more failures closer to the beginning than the end as those who persist at the beginning are less likely to give up with the passage of time. Here is a way using the -addplot()- option.

webuse drugtr, clear
stset studytime, failure(died)
set scheme s1mono
*INDICATORS (MARKERS)
sts graph, addplot(scatter _d _t if _d==1, msy(d)) leg(off)

*NUMBER OF FAILURES
bys _t: egen howmany= total(_d)
sts graph, addplot(scatter _d _t if _d==1, msy(none) mlab(howmany) mlabpos(0)) leg(off)

[Nishan Lamichhane]

Thank you Andrew. It helped a lot.