11/12 part2

data life;

input race$ gender $ belief$ count ;
datalines;
white f yes 371
white f und 49
white f no 74
white m yes 250
white m und 45
white m no 71
black f yes 64
black f und 9
black f no 15
black m yes 25
black m und 5
black m no 13
;
run;
proc print data=life;
run;


proc logistic data=life order=data;
class race gender;
freq count;
model belief= race gender /link=glogit aggregate scale=none;
output out=prob predprobs=I;
run;


proc logistic data=life order=data;
class race (ref="black")
gender (ref="m")/ param=ref;
freq count;
model belief= race gender /link=glogit aggregate scale=none;
output out=prob predprobs=I;
run;