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

use rrdat1, clear

gen     des = 2                          /*lateral moves*/
replace des = 1 if (presn/pres -1)>0.2   /*upward moves*/
replace des = 3 if (presn/pres -1)<0.0   /*downward moves*/
replace des = 0 if presn==-1

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*/

*define single episode data and fit parametric survival model 

stset tf, fail(des==1)   /*upward*/        
streg edu coho2 coho3 lfx pnoj pres if pres <=65, dist(e) nohr   

stset tf, fail(des==2)   /*lateral*/   
streg edu coho2 coho3 lfx pnoj pres, dist(e) nohr   

stset tf, fail(des==3)   /*downward*/
streg edu coho2 coho3 lfx pnoj pres, dist(e) nohr 

log close