options ls=76;
*** example 9.3 of Armitage Berry ***;
data AB93; 
input group $ group2 $ x y;
cards; 
A A1 39 4.62 
A A1 40 5.29 
A A1 41 5.52
A A1 41 3.71  
A A1 45 4.02 
A A1 49 5.09
A A1 52 2.70  
A A1 47 4.31     
A A1 61 2.70   
A A1 65 3.03       
A A1 58 2.73      
A A1 59 3.67     
A A2 29 5.21    
A A2 29 5.17   
A A2 33 4.88  
A A2 32 4.50 
A A2 31 4.47
A A2 29 5.12 
A A2 29 4.51  
A A2 30 4.85 
A A2 21 5.22
A A2 28 4.62
A A2 23 5.07
A A2 35 3.64
A A2 38 3.64    
A A2 38 5.09   
A A2 43 4.61  
A A2 39 4.73          
A A2 38 4.58         
A A2 42 5.12        
A A2 43 3.89       
A A2 43 4.62      
A A2 37 4.30     
A A2 50 2.70    
A A2 50 3.50
A A2 45 5.06
A A2 48 4.06
A A2 51 4.51
A A2 46 4.66
A A2 58 2.88
B B  27 5.29
B B  25 3.67
B B  24 5.82
B B  32 4.77
B B  23 5.71
B B  25 4.47
B B  32 4.55
B B  18 4.61
B B  19 5.86
B B  26 5.20
B B  33 4.44
B B  27 5.52
B B  33 4.97
B B  25 4.99
B B  42 4.89
B B  35 4.09
B B  35 4.24
B B  41 3.88
B B  38 4.85
B B  41 4.79
B B  36 4.36
B B  36 4.02
B B  41 3.77
B B  41 4.22
B B  37 4.94
B B  42 4.04
B B  39 4.51
B B  41 4.06
B B  43 4.02
B B  41 4.99
B B  48 3.86
B B  47 4.68
B B  53 4.74
B B  49 3.76
B B  54 3.98
B B  48 5.00
B B  49 3.31
B B  47 3.11
B B  52 4.76
B B  58 3.95
B B  62 4.60
B B  65 4.83
B B  62 3.18
B B  59 3.03
;
run;
proc glm; class group group2;
model y=x group;
run;
proc glm; class group;
model y=x group x*group;
run;
proc glm; class group2;
model y=x group2;
run;




data temp;
set  AB93; 
if (y > 4) then bin=1;
if (y <= 4) then bin=0;
run;
proc contents data = temp;
run;
proc logistic data=temp;
class group;
model bin=x group;
run;