In medical and surgical research the hazard function is often used to estimate risk of an event across time. However the assumption of constant hazard is often not met and is not ideal when analyzing repeated events such are valve replacements, reoperations, reinterventions, multiple episodes of infection or repeated transplant rejections. I have read some information on modulated renewal theory and the Nelson estimator which involve calculating interfailure event times and segments/gaps between events within the same patient. Can anyone suggest a useful and relatively simple way of approaching this in SAS, Stata or R? Thank you!