* Data file:  addicts_sub.dat
*
* survival data matrix for westmead study - missing doses
* follow-up until 1/1/89 but corrected during 1/7/90 follow-up
* Column 1: Addict ID
*        2: Clinic (1 or 2)
*        3: Censor status (0=censored, 1=complete)
*        4: Survival time (days)
*        5: Sex (0=male, 1=female)
*        6: Age
*        7: Employment status
*        8: History of benzodiazepine/barbiturate abuse (0=no, 1=yes)
*        9: History of amphetamine/cocaine abuse (0=no, 1=yes)
*       10: History of alcohol abuse (0=no, 1=yes)
*       11: Maximum dose
*

infile id clinic status time female age employ benzo amph alc dose using addicts_sub.dat

label variable time "time (days)"
label variable status "status"
label variable clinic "study clinic"
label variable dose "methadone dose"

label define ylab 0 "no" 1 "yes"
label values employ ylab
label values benzo ylab
label values amph ylab
label values alc ylab

*** recode clinic ***
recode clinic 1=0 2=1
label define clab 0 "clinic 1" 1 "clinic 2"
label values clinic clab

***
*** predictors:  female, age, employ, benzo, amph, alc
***

generate dosecat = dose
recode dosecat min/59=1 60/79=2 80/max=3
label define dlab 1 "dose <= 59" 2 "dose 60-79" 3 "80 <= dose"
label values dosecat dlab

generate agecat = age
recode agecat min/24.999=1 25/29.999=2 30/max=3
label define alab 1 "age < 25" 2 "age 25-29" 3 "30 <= age"
label values agecat alab