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/clodds=wald;
output out=m2 p=prob xbeta=logit;
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';
title2 'from Table 9.14 P. 338 of Woolson and Clarke';
run;

proc sort data=m2;
by medage;
run;

proc gplot data=m2;
 plot logit*medage;
 plot prob*medage;
 title 'maoidata with logistic Regression';
run;

quit;