*NKNW878.sas, two-way anova with one observation per cell
   using data in Table 21.2;

options nocenter;
data a1; infile 'c:\Ageorge\Purdue\Stat512\Textdata\CH21TA02.DAT';
   input premium size region;
if size=1 then sizea='1_small ';
if size=2 then sizea='2_medium';
if size=3 then sizea='3_large ';
proc print data=a1; run;

proc glm data=a1;
   class sizea region;
   model premium=sizea region/solution;
   means sizea region / tukey;
   output out=a2 p=muhat;
run;
proc print data=a2; run;

symbol1 v='E' i=join c=black;
symbol2 v='W' i=join c=black;
title1 'Plot of the data';
proc gplot data=a2; 
   plot premium*sizea=region/frame;
run;

symbol1 v='E' i=join c=black;
symbol2 v='W' i=join c=black;
title1 'Plot of the model estimates';
proc gplot data=a2; 
   plot muhat*sizea=region/frame;
run;