options nodate nonumber ps=200 ls=80 formdlim=' ';

/* Example 1: Horseshoe Crab Revisited */
data crabs;
input width cases satell; 
cards;
22.69 14  5
23.84 14  4
24.77 28 17
25.84 39 21
26.79 22 15
27.74 24 20 
28.67 18 15 
30.41 14 14 
;
proc genmod;  
  model satell/cases = width / dist=bin link=identity;
run;                 
proc genmod; /* lrci requests that two-sided confidence intervals for all model parameters be computed based on the profile likelihood function. */
   model satell/cases = width / dist=bin link=logit waldci lrci alpha=.01; 
run;                 
proc logistic;       
   model satell/cases = width / influence risklimits; 
output out=predict  p=pi_hat  lower=LCL  upper=UCL  pred=prob;
proc print  data=predict;
run;

   /* Example 1 (Cont'd) */

proc genmod data=crabs;    
   model satell/cases = width / dist=bin link=logit obstats residuals waldci lrci covb; 
run;                

proc logistic data=crabs;       
   model satell/cases = width / influence; 
output out=predict  p=pi_hat  lower=LCL  upper=UCL;
proc print  data=predict;
run;