capture log close oxy1 *F:\OxyHeroin\cox1.do clear set more off set mem 300m use "F:\OxyHeroin\oxyherfin.dta" compress gen oxstatus=3 replace oxstatus=1 if afult18 ==1 replace oxstatus=2 if afult18 ==2 gen a1217=0 replace a1217=1 if yeligher3 ==1 & newage==1 gen a1825=0 replace a1825=1 if yeligher3 ==1 & newage==2 gen a2649=0 replace a2649=1 if yeligher3 ==1 & newage==3 gen coxpop=0 replace coxpop=1 if age >12 & age <50 & ynewoxyev3 ==1 char nevmar [omit] 2 char ynewoxyev3 [omit] 2 char afult18 [omit] 2 char rurorurb [omit] 2 char newage [omit] 3 char irsex [omit] 2 char white [omit] 2 char afult18 [omit] 2 char ilt20k [omit] 2 char xmstatus [omit] 4 char oxstatus [omit] 3 char xmrjflag [omit] 2 char xcocflag [omit] 2 char xinhflag [omit] 2 char xhalflag [omit] 2 char xanlflag [omit] 2 tab1 ynewoxyev3 newage nevmar xanlflag t_event if yeligher3 ==1 log using "F:\OxyHeroin\cox_models.log", name (oxy1) replace svyset verep[pweight=awt0408], strata (vestr) singleunit(scaled) stset t_event [pweight=awt0408], fail (yrecher3=1) foreach indvar of varlist newage irsex white nevmar rurorurb ilt20k afult18 oxstatus ynewoxyev3 xmrjflag xcocflag /// xinhflag xhalflag xanlflag { capture noisily xi: svy, subpop(yeligher3): stcox i.`indvar' } svyset verep[pweight=awt0408], strata (vestr) singleunit(scaled) stset t_event [pweight=awt0408], fail (yrecher3=1) capture noisily xi: svy, subpop(yeligher3): stcox i.newage i.nevmar svyset verep[pweight=awt0408], strata (vestr) singleunit(scaled) stset t_event [pweight=awt0408], fail (yrecher3=1) capture noisily xi: svy, subpop(yeligher3): stcox i.newage i.nevmar i.xcocflag i.xanlflag svyset verep[pweight=awt0408], strata (vestr) singleunit(scaled) stset t_event [pweight=awt0408], fail (yrecher3=1) capture noisily xi: svy, subpop(yeligher3): stcox i.ynewoxyev3 i.newage i.nevmar i.xanlflag log close oxy1