** Armitage and Berry example 12.1*; 
data AB121;
input score yes num;
cards;
0  59  156
1  10  41
2  12  48
3  5   33
;
proc logistic data=AB121;
class score; 
model yes/num=score;
run;


data AB121N;
input score visit $ num;
cards;
0   yes 59
0   no  97
1   yes 10
1   no  31
2   yes 12 
2   no  36
3   yes  5 
3   no  28
;
proc freq; 
weight num;
tables score * visit / exact;
run;

proc logistic 
model 

proc catmod; direct score;
model yes/num=score;
run;



* Example of table 4.2 of Agresti <<Categorical Data Analysis>> *
* remiss = 0 not achieve Remission in Cancer patients           *
* remiss = 1 Remission in Cancer patients                       *
* li = percentage of cell that are "labeled" after              *
* injection of tritiated thymidine                              *;

options ls=70;
data logreg;
input remiss li;
cards;
0 8 
0 8 
0 10
0 10
0 12
0 12
0 12
0 14
0 14
0 14
0 16
0 16
0 16
1 18
1 20
1 20
0 20
1 22
0 22
0 24
1 26
1 28
0 32
1 34
1 38
1 38
0 38
;
run;
proc logistic;
model remiss=li;
run;



proc logistic descending;
model remiss=li;
run; 


* adapted data from Table 4.4 of Agresti *
* of Framingham Men Blood pressure Data  *;

data two;
input bp num hd;
cards;
112  156 3
122  252 17
132  284 12
142  271 16
152  139 12
162   85  8
177   99 16
197   43  8 
180    .  . 
;
run;
proc logistic data=two descending;
model hd/num=bp;
output out=rates predicted=phat lower=phatlow upper=phatup;
run;

proc print data=rates;
run;