data crab;
input color width satell;
color = color-1;
if satell>0 then y=1; if satell=0 then y=0;
datalines;
3  28.3  8  
4  22.5  0  
2  26.0  9  
4  24.8  0  
4  26.0  4  
3  23.8  0  
2  26.5  0  
4  24.7  0  
3  23.7  0  
4  25.6  0  
4  24.3  0  
3  25.8  0  
3  28.2  11 
5  21.0  0  
3  26.0  14 
2  27.1  8  
3  25.2  1  
3  29.0  1  
5  24.7  0  
3  27.4  5  
3  23.2  4  
2  25.0  3  
3  22.5  1  
4  26.7  2  
5  25.8  3  
5  26.2  0  
3  28.7  3  
3  26.8  5  
5  27.5  0  
3  24.9  0  
2  29.3  4  
2  25.8  0  
3  25.7  0  
3  25.7  8  
3  26.7  5  
5  23.7  0  
3  26.8  0  
3  27.5  6  
5  23.4  0  
3  27.9  6  
4  27.5  3  
2  26.1  5  
2  27.7  6  
3  30.0  5  
4  28.5  9  
4  28.9  4  
3  28.2  6  
3  25.0  4  
3  28.5  3  
3  30.3  3  
5  24.7  5  
3  27.7  5  
2  27.4  6  
3  22.9  4  
3  25.7  5  
3  28.3  15 
3  27.2  3  
4  26.2  3  
3  27.8  0  
5  25.5  0  
4  27.1  0  
4  24.5  5  
4  27.0  3  
3  26.0  5  
3  28.0  1  
3  30.0  8  
3  29.0  10 
3  26.2  0  
3  26.5  0  
3  26.2  3  
4  25.6  7  
4  23.0  1  
4  23.0  0  
3  25.4  6  
4  24.2  0  
3  22.9  0  
4  26.0  3  
3  25.4  4  
4  25.7  0  
3  25.1  5  
4  24.5  0  
5  27.5  0  
4  23.1  0  
4  25.9  4  
3  25.8  0  
5  27.0  3  
3  28.5  0  
5  25.5  0  
5  23.5  0  
3  24.0  0  
3  29.7  5  
3  26.8  0  
5  26.7  0  
3  28.7  0  
4  23.1  0  
3  29.0  1  
4  25.5  0  
4  26.5  1  
4  24.5  1  
4  28.5  1  
3  28.2  1  
3  24.5  1  
3  27.5  1  
3  24.7  4  
3  25.2  1  
4  27.3  1  
3  26.3  1  
3  29.0  1  
3  25.3  2  
3  26.5  4  
3  27.8  3  
3  27.0  6  
4  25.7  0  
3  25.0  2  
3  31.9  2  
5  23.7  0  
5  29.3  12 
4  22.0  0  
3  25.0  5  
4  27.0  6  
4  23.8  6  
2  30.2  2  
4  26.2  0  
3  24.2  2  
3  27.4  3  
3  25.4  0  
4  28.4  3  
5  22.5  4  
3  26.2  2  
3  24.9  6  
2  24.5  6  
3  25.1  0  
3  28.0  4  
5  25.8  10 
3  27.9  7  
3  24.9  0  
3  28.4  5  
4  27.2  5  
3  25.0  6  
3  27.5  6  
3  33.5  7  
3  30.5  3  
4  29.0  3  
3  24.3  0  
3  25.8  0  
5  25.0  8  
3  31.7  4  
3  29.5  4  
4  24.0  10 
3  30.0  9  
3  27.6  4  
3  26.2  0  
3  23.1  0  
3  22.9  0  
5  24.5  0  
3  24.7  4  
3  28.3  0  
3  23.9  2  
4  23.8  0  
4  29.8  4  
3  26.5  4  
3  26.0  3  
3  28.2  8  
5  25.7  0  
3  26.5  7  
3  25.8  0  
4  24.1  0  
4  26.2  2  
4  26.1  3  
4  29.0  4  
2  28.0  0  
5  27.0  0  
3  24.5  0  
;

proc genmod descending; class color;
model y = color width / dist=bin link=logit type3;
run;

proc genmod descending;
model y = color width / dist=bin link=logit type3;
run;



data crab;
input color width satell;
color = color-1;
if color=4 then dark=0;
else dark=1;
if satell>0 then y=1; if satell=0 then y=0;
datalines;
3  28.3  8  
4  22.5  0  
2  26.0  9  
4  24.8  0  
4  26.0  4  
3  23.8  0  
2  26.5  0  
4  24.7  0  
3  23.7  0  
4  25.6  0  
4  24.3  0  
3  25.8  0  
3  28.2  11 
5  21.0  0  
3  26.0  14 
2  27.1  8  
3  25.2  1  
3  29.0  1  
5  24.7  0  
3  27.4  5  
3  23.2  4  
2  25.0  3  
3  22.5  1  
4  26.7  2  
5  25.8  3  
5  26.2  0  
3  28.7  3  
3  26.8  5  
5  27.5  0  
3  24.9  0  
2  29.3  4  
2  25.8  0  
3  25.7  0  
3  25.7  8  
3  26.7  5  
5  23.7  0  
3  26.8  0  
3  27.5  6  
5  23.4  0  
3  27.9  6  
4  27.5  3  
2  26.1  5  
2  27.7  6  
3  30.0  5  
4  28.5  9  
4  28.9  4  
3  28.2  6  
3  25.0  4  
3  28.5  3  
3  30.3  3  
5  24.7  5  
3  27.7  5  
2  27.4  6  
3  22.9  4  
3  25.7  5  
3  28.3  15 
3  27.2  3  
4  26.2  3  
3  27.8  0  
5  25.5  0  
4  27.1  0  
4  24.5  5  
4  27.0  3  
3  26.0  5  
3  28.0  1  
3  30.0  8  
3  29.0  10 
3  26.2  0  
3  26.5  0  
3  26.2  3  
4  25.6  7  
4  23.0  1  
4  23.0  0  
3  25.4  6  
4  24.2  0  
3  22.9  0  
4  26.0  3  
3  25.4  4  
4  25.7  0  
3  25.1  5  
4  24.5  0  
5  27.5  0  
4  23.1  0  
4  25.9  4  
3  25.8  0  
5  27.0  3  
3  28.5  0  
5  25.5  0  
5  23.5  0  
3  24.0  0  
3  29.7  5  
3  26.8  0  
5  26.7  0  
3  28.7  0  
4  23.1  0  
3  29.0  1  
4  25.5  0  
4  26.5  1  
4  24.5  1  
4  28.5  1  
3  28.2  1  
3  24.5  1  
3  27.5  1  
3  24.7  4  
3  25.2  1  
4  27.3  1  
3  26.3  1  
3  29.0  1  
3  25.3  2  
3  26.5  4  
3  27.8  3  
3  27.0  6  
4  25.7  0  
3  25.0  2  
3  31.9  2  
5  23.7  0  
5  29.3  12 
4  22.0  0  
3  25.0  5  
4  27.0  6  
4  23.8  6  
2  30.2  2  
4  26.2  0  
3  24.2  2  
3  27.4  3  
3  25.4  0  
4  28.4  3  
5  22.5  4  
3  26.2  2  
3  24.9  6  
2  24.5  6  
3  25.1  0  
3  28.0  4  
5  25.8  10 
3  27.9  7  
3  24.9  0  
3  28.4  5  
4  27.2  5  
3  25.0  6  
3  27.5  6  
3  33.5  7  
3  30.5  3  
4  29.0  3  
3  24.3  0  
3  25.8  0  
5  25.0  8  
3  31.7  4  
3  29.5  4  
4  24.0  10 
3  30.0  9  
3  27.6  4  
3  26.2  0  
3  23.1  0  
3  22.9  0  
5  24.5  0  
3  24.7  4  
3  28.3  0  
3  23.9  2  
4  23.8  0  
4  29.8  4  
3  26.5  4  
3  26.0  3  
3  28.2  8  
5  25.7  0  
3  26.5  7  
3  25.8  0  
4  24.1  0  
4  26.2  2  
4  26.1  3  
4  29.0  4  
2  28.0  0  
5  27.0  0  
3  24.5  0  
;

proc genmod descending order=data; class dark;
model y = dark width / dist=bin link=logit type3;
run;

proc genmod descending order=data; class dark;
model y = dark width dark*width / dist=bin link=logit type3;
run;