options ls=72;
data a;
input li n r @@;
pihat=r/n;
cards;
8 2 0 10 2 0 12 3 0 14 3 0 16 3 0 18 1 1 20 3 2 22 2 1 24 1 0
26 1 1 28 1 1 32 1 0 34 1 1 38 3 2
proc plot;
plot pihat*li;;
proc logistic;
model r/n=li;
output out=esta p=phata reschi=resida;
proc print; var li phata;
proc plot; plot resida*li;
proc logistic;
model r/n=li/link=normit;
output out=estb p=phatb reschi=residb;
proc print; var li phatb;
proc logistic;
model r/n=li/link=cloglog;
output out=estc p=phatc reschi=residc;
proc print; var li phatc;
run;