************************************************************************** * * * athens.do * * * * PURPOSE: analysis of Tzonou et al. (1998) data investigating * * CHD and iron consumption * * * * DATE: 11 April 2002 * * * ************************************************************************** *** *** read the data and format *** infile female age iron case count using athens.dat label define fnames 0 "male" 1 "female" label values female fnames label define anames 1 "<= 49" 2 "50-59" 3 "60-69" 4 ">= 70" label values age anames label define inames 1 "<= 250mg/month" 2 "251-300 mg/month" 3 "301-350 mg/month" 4 "351-400 mg/month" 5 "> 400 mg/month" label values iron inames label define ynames 0 "control" 1 "case" label values case ynames *** *** create a single record per subject *** drop if count==0 expand count *** *** recode Iron *** generate newiron=iron recode newiron 1/3=0 4/5=1 label define nnames 0 "<=350 mg/month" 1 "> 350 mg/month" label values newiron nnames *** *** univariate and bivariate distributions *** table case table newiron table age table female ***** Note: we use column summaries since it's a case-control study tabulate newiron case, col tabulate age case, col tabulate female case, col tabulate newiron female, row tabulate newiron age, row *** *** crude analysis *** cc case newiron *** *** adjusted analysis *** cc case newiron, by(age) bd cc case newiron, by(female) bd mhodds case newiron age female mhodds case newiron female, by(age)