*NKNW884.sas, two-way anova with one observation per cell,
   Tukey one df for additivity, using data in Table 21.2;
options nocenter;

data a1; infile 'c:\Ageorge\Purdue\Stat512\Textdata\CH21TA02.DAT';
   input premium size region;
proc print data=a1; run;

proc glm data=a1; 
   model premium=;
   output out=aall p=muhat;
proc print data=aall;
run;

proc glm data=a1; 
   class size;
   model premium=size;
   output out=aA p=muhatA;
proc print data=aA;
run;

proc glm data=a1; 
   class region;
   model premium=region;
   output out=aB p=muhatB;
proc print data=aB;
run;

data a2; merge aall aA aB;
   alpha=muhatA-muhat;
   beta=muhatB-muhat;
   atimesb=alpha*beta;
proc print data=a2;
   var size region atimesb;
run;

proc glm data=a2;
   class size region;
   model premium=size region atimesb/solution;
run;