options nodate pagesize=54 linesize=80; /* Read in the input file */ data ex0507; * infile '../../../data/biostat/ex0507.dat'; infile 'C:\biostat\data\ex0507.dat'; input age smoke death numcases smoke_r death_r; run; /* Obtain crosstabs and mantel-haenszel test based on reversed smoking and death (SAS default is to model the probability in the "lowest" category) */ /* We weight the cases by the variable numcases which represents number of people in each table cell */ proc freq; table age*smoke_r*death_r / cmh; weight numcases; run; quit;