version 9
set scheme sj
capture log close
set more off
log using ehg1.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(gomp) nohr         /*fit parametric survival model*/

stcurve, hazard ytick(0(0.005)0.02) ylabel(0(0.01)0.02)   /*gompertz transition rate*/

log close