data prop_8_1;
do G='M', 'F';
 do I='S', 'O';
  do H='S', 'O';
   input count @;
   output;
  end;
 end;
end;
datalines;
76 160 6 25 114 181 11 48
;
proc genmod; class G H I;
 model count=G H I G*H G*I / dist=poisson link=log;
proc genmod; class G H I;
 model count=G H I G*H H*I / dist=poisson link=log;
proc genmod; class G H I;
 model count=G H I G*I H*I / dist=poisson link=log;
proc genmod; class G H I;
 model count=G H I G*H G*I H*I / dist=poisson link=log;
run;

data prop_8_5;
do use='S', 'N';
 do ejected='Y', 'N';
  do injury='N', 'F';
   input count @;
   output;
  end;
 end;
end;
datalines;
1105 14 411111 483 4624 497 157342 1008
;
proc genmod; class use ejected injury;
 model count=use ejected injury /dist=poisson link=log;
 run;
proc genmod; class use ejected injury;
 model count=use|ejected|injury /dist=poisson link=log;
 run;
proc genmod; class use ejected injury;
 model count=use|ejected|injury @2 /dist=poisson link=log;
 run;
data prop_8_5;
do use='S', 'N';
 do ejected='Y', 'N';
   input non_fatal_n fatal_n @;
   total_n = non_fatal_n + fatal_n;
   output;
  end;
 end;
datalines;
1105 14 411111 483 4624 497 157342 1008
;
proc print;
run;

 proc genmod; class use ejected;
 model fatal_n / total_n =use ejected /dist=binomial link=logit ;
 run;