/* Analysis of a 3-way table Berkeley Admissions data using PROC FREQ */

options nocenter nodate nonumber linesize=72;
data berkeley;
   input D $ S $ A $ count;
   cards;
DeptA  Male    Reject  313
DeptA  Male    Accept  512
DeptA  Female  Reject   19
DeptA  Female  Accept   89
DeptB  Male    Reject  207
DeptB  Male    Accept  353
DeptB  Female  Reject    8
DeptB  Female  Accept   17
DeptC  Male    Reject  205
DeptC  Male    Accept  120
DeptC  Female  Reject  391
DeptC  Female  Accept  202
DeptD  Male    Reject  279
DeptD  Male    Accept  138
DeptD  Female  Reject  244
DeptD  Female  Accept  131
DeptE  Male    Reject  138
DeptE  Male    Accept   53
DeptE  Female  Reject  299
DeptE  Female  Accept   94
DeptF  Male    Reject  351
DeptF  Male    Accept   22
DeptF  Female  Reject  317
DeptF  Female  Accept   24
;

/*analysis of the three-way table including CMH test*/
proc freq data=berkeley order=data;
weight count;
tables D*S*A/ cmh chisq relrisk expected nocol norow;
tables S*A/chisq relrisk;
run;