OPTION LS=76 NOCENTER;
DATA hd;
INPUT bp disease n;
r=disease/n;y=log(r/(1-r));z=arsin(r);
DATALINES;
111.5  3  156
121.5  17 252
131.5  12 284
141.5  16 271
151.5  12 139
161.5  8  85
176.5  16 99
191.5  8  43
;
PROC GENMOD;
MODEL disease/n = bp / DIST = bin LINK = logit RESIDUALS;
title 'genmodm output for herat disease example, Agresti, p. 221';
run;

data pred;set hd;
p=1/(1+exp(6.082-0.0243*bp));np=n*p;
x=(disease-np)/(np*(1-p))**0.5;
PROC PRINT;VAR bp n disease np x;
TITLE 'To check the residual and X-square for chi-square test';
PROC MEANS N USS;VAR x;



PROC GLM DATA=hd;
MODEL y=bp/P;
title 'GLM output with logit tranformation';

PROC GLM DATA=hd;
MODEL z=bp/P;
title 'GLM output with arc sine tranformation';
run;