SAS for Poisson and Negative Binomial GLMs for Crab Data of Table 4.3 
Agresti, p. 127 (number of male satellites by female's characteristics)
----------------------------------------------------------------------
data crab;
input color spine width satell weight; 
   weight=weight/1000; color=color-1;
datalines;

also try y=1 if satellites, 0 if not