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

*Provides Wald based CIs in output by default;
proc genmod data=crab_logist descending; 
model Y = width / dist=binomial link=logit; 
run;

*Provides LRT based CIs using LRCI option;
proc genmod data=crab_logist descending; 
model Y = width / dist=binomial link=logit lrci; 
run;

*CONTRAST statement for LRT of beta=0;
proc genmod data=crab_logist descending; 
model Y = width / dist=binomial link=logit; 
contrast "LRT OUTPUT" width 1;
run;

*Exclude width as explanatory variable and
 use deviances for full and reduced models
 to generate the LRT of beta=0;
proc genmod data=crab_logist descending; 
model Y = / dist=binomial link=logit; 
run;