option ls=76 nocenter;
data clinics;
   input clinic trt$ fav unfav;
   count=fav;response=1;output;
   count=unfav;response=0;output;
   /* if 'trt'='drug' then treat=1;else treat=0; */
   datalines;
1 drug 11 25  
1 cntl 10 27 
2 drug 16 4 
2 cntl 22 10 
3 drug 14 5 
3 cntl 7 12 
4 drug 2 14 
4 cntl 1 16 
5 drug 6 11 
5 cntl 0 12 
6 drug 1 10 
6 cntl 0 10 
7 drug 1 4 
7 cntl 1 8 
8 drug 4 2 
8 cntl 6 1 
;
PROC PRINT;
PROC FREQ DATA=clinics;
WEIGHT count;
TABLES clinic*trt*response/CMH EXPECTED CHISQ;
TITLE 'When cmh option in PROC FREQ is used for data analysis';
TITLE2 'Method suggested by Agresti, Categorical Data book, p. 640';
run;