version 9
set scheme sj
capture log close
set more off
log using ehc9.log, replace

use rrdat1, clear

gen     des = 2                          /*destination state*/ 
replace des = 1 if (presn/pres -1)>0.2   
replace des = 3 if (presn/pres -1)<0.0 
replace des = 0 if presn==-1

gen tf = tfin - tstart + 1               /*ending time*/

stset tf, f(des=1) /* upward moves */
sts list if pres <= 65, at(0 2 4 6/8 160 170 326 428) 
sts gen surv1 = s

stset tf, f(des=2) /* lateral moves */
sts list, at(0 3 4/7 184 194 209 350 428) 
sts gen surv2 = s

stset tf, f(des=3) /* downward moves */
sts list, at(0 2/6 275 293 312 332 428) 
sts gen surv3 = s

for any surv1 surv2 surv3 \ any "upward" "lateral" "downward": label var X "Y" 

graph twoway line surv1 surv2 surv3 _t, sort  /// 
ysc(r(0 1)) ylabel(0(0.2)1) xtitle("analysis time")legend(row(1)) 

log close