data accidents;
input gender$ location$ seatbelt$ y1-y5;
resp=1; count=y1; output;
resp=2; count=y2; output;
resp=3; count=y3; output;
resp=4; count=y4; output;
resp=5; count=y5; output;
drop y1-y5;
datalines;
female urban no 7287 175 720 91 10
female urban yes 11587 126 577 48 8
female rural no 3246 73 710 159 31
female rural yes 6134 94 564 82 17
male urban no 10381 136 566 96 14
male urban yes 10969 83 259 37 1
male rural no 6123 141 710 188 45
male rural yes 6693 74 353 74 12
;
proc print;run;
proc genmod data=accidents;
class gender location seatbelt;
model resp = gender location seatbelt location*seatbelt/ dist=multinomial link=clogit lrci type3;
freq count;
run;