Generalized Linear Model for Binary Data - Reocities

data obstats1;
set obstat1 (rename=(yvar1=y yvar2=n));
pred1=pred*n;
run;
/* Using Probit model */
proc genmod data=binary;
make 'OBSTATS' out=obstat2;
model killed/n = dose / dist=bin
link=probit
lrci
waldci
obstat;
run;
data obstats2;
set obstat2 (rename=(yvar1=y yvar2=n));
pred2=pred*n;
run;
/* Using extrme value model */
proc genmod data=binary;
make 'OBSTATS' out=obstat3;
model killed/n = dose / dist=bin
link=cll
lrci
waldci
obstat;
run;
data obstats3;
set obstat3 (rename=(yvar1=y yvar2=n));
pred3=pred*n;
run;
data comp;
merge obstats1 obstats2 obstats3;
keep y n pred1 pred2 pred3;
run;
proc print;
title ‘Comparison of various dose-response models for the beetle
mortality data’;
run;