version 9 capture log close set more off log using ehd1.log, replace use rrdat1, clear gen des = tfin ~= ti /*destination state*/ gen tf = tfin - tstart + 1 /*ending time*/ stset tf, failure(des) /*define single episode data*/ streg, dist(e) nohr /*fit parametric survival model*/ log close