infile age alc tob y count using NewTuyns.dat

***
*** variable labels
***

label variable age "Age Group"
label variable alc "Alcohol"
label variable tob "Tobacco"
label variable y "Case/Control Status"

label define agegps 1 "25-34" 2 "35-44" 3 "45-54" 4 "55-64" 5 "65-74" 6 "75+"
label define alcgps 1 "<40g/day" 2 "40-79g/day" 3 "80-119g/day" 4 "120+g/day"
label define tobgps 1 "0-9g/day" 2 "10-19g/day" 3 "20-29g/day" 4 "30+g/day"
label define status 1 "Case" 0 "Control"

label values age agegps
label values alc alcgps
label values tob tobgps
label values y status

***
*** Univariate summaries
***
tabulate age [freq=count]
tabulate tob [freq=count]
tabulate alc [freq=count]

***
*** Bivariate summaries
***
tabulate y tob [freq=count], chi2 row
tabulate y alc [freq=count], chi2 row

***
*** test for trend
***
tabodds y tob [freq=count], or
tabodds y alc [freq=count], or

***
*** display row/column percents
***
tabulate age y [freq=count], chi2 col
tabulate age alc [freq=count], chi2 row
tabulate age tob [freq=count], chi2 row