version 9
capture log close
set more off
log using ehd7.log, replace

use rrdat1, clear

keep if noj <= 4   /*keep only episodes with serial number <= 4*/

gen des = tfin ~= ti               /*destination state*/
gen tf  = tfin - tstart + 1        /*ending time*/

gen coho2 = tb >= 468 & tb <= 504  /*cohort 2*/
gen coho3 = tb >= 588 & tb <= 624  /*cohort 3*/
gen lfx   = tstart - te            /*labor force experience*/
gen pnoj  = noj - 1                /*previous number of jobs*/

stset tf, failure(des)          /*define single episode data*/

stci, by (noj des) emean        /*descriptive statistics for the individuals' first four jobs*/

by noj, sort: streg edu coho2 coho3 lfx pres, dist(e) nohr   /*fit parametric survival models*/ 


log close