option ls=76 nocenter;
data clinics;
   input clinic trt$ fav unfav;
   n=fav+unfav;y=log((fav+0.5)/(unfav+0.5));
   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 genmod;class clinic trt;
model fav/n=clinic trt/dist=bin link=logit RESIDUALS;
TITLE 'When PROC GENMOD is used with trt as discrete';

DATA detail;SET clinics;
array cli(8) cli1-cli8;
DO i=1 to 8;
IF clinic=i THEN cli(i)=1; ELSE cli(i)=0;
END;
treat=0;IF trt='drug' THEN treat=1;
PROC PRINT;VAR clinic cli1-cli8 trt treat fav unfav;
TITLE 'Using PROC GENMOD with hospital and trt treated as continuous';
proc genmod;
model fav/n=cli1-cli8 treat/NOINT dist=bin link=logit RESIDUALS;

PROC MIXED DATA=clinics;CLASS clinic trt;
MODEL y=trt/ddfm=satterth;
RANDOM clinic;
TITLE 'PROC MIXED with logit transformation';


run;