data maoi;
input medage improve total;
label medage='median age'
      improve='number improved'
      total='total number of patients';
lines;
30 3 27
34 3 21
38 6 32
42 6 19
46 11 26
50 12 25
54 10 20
58 12 23
62 12 17
run;       



proc logistic;
 model improve/total=medage/influence;
output out=m2 p=prob lower=LCL upper=UCL;
run;

Proc print data=m2;

run;

title 'example for logistic regression with residual diagnosgtics';
title2 'from Table 9.14 P. 338 of Woolson and Clarke';
run;

proc genmod;
model improve/total=medage /dist=bin link =logit obstats
residuals walci lrci;

run;

quit;