data one;
   do A=0,1;
      do B=0,1;
         do Y=1,2;
            input F @@;
            output;
         end;
      end;
   end;
   datalines;
23 63 31 70 67 100 70 104
;
run;



proc print data = one;
run;



*The following statements fit a full model to examine the main effects of A and B as well as the interaction effect of A and B;


proc logistic data=one;
   freq F;
   model Y=A B A*B;
run;


/* Now checking GOFs using the "AGGREGATE" option */

proc logistic data=one;
   freq F;
  model Y=A/ scale=none aggregate = (A B);
  *model Y=A/ scale=none aggregate;
  *model Y=A/ scale=none aggregate = (A);
  *model Y=A/ scale=none aggregate = (B);
run;