/********************************************************
 THIS COMMAND FILE READS IN THE ALLIGATOR DATA SET
 AND DOES SOME SIMPLE DESCRIPTIVE STATISTICS.
 EXAMPLE OF LOG-LINEAR MODEL.

 FILENAME: ALLIGATOR.SAS
**************************************************************/
data alligat;
  input lake size food count;
  cards;
1 1 1 23
1 1 2 4
1 1 3 12
1 2 1 7
1 2 2 0
1 2 3 9
2 1 1 5 
2 1 2 11
2 1 3 4
2 2 1 13
2 2 2 8
2 2 3 7
3 1 1 5
3 1 2 11
3 1 3 8
3 2 1 8
3 2 2 7
3 2 3 14
4 1 1 16 
4 1 2 19
4 1 3 8
4 2 1 17
4 2 2 1
4 2 3 4
;
proc freq;
  tables size*lake*food / nocol nopercent norow;
  weight count;
run;

/*Lab Ex 5: Saturated Model*/

proc genmod data=alligat;
  class lake size food;
  model count = lake | size | food / dist=poisson type3;
  title "Saturated Model";
run;

/*Lab Ex 5: Model with 2-way interactions*/

proc genmod data=alligat;
  class lake size food;
  model count = lake  size lake*size  food 
                lake*food size*food / dist=poisson type3;
  output out=preddat2 p=predict;
  title "2-Way Interactions Model";
run;

proc print data=preddat2;
run;