*NKNW817b.sas, two-way anova using data in Table 19.7;
options nocenter;
data a1; infile 'c:\Ageorge\Purdue\Stat512\Textdata\CH19TA07.DAT';
   input sales height width;
proc print data=a1; run;

proc glm data=a1;
   class height width;
   model sales=height width height*width/solution;
   means height*width;
run;

proc glm data=a1;
   class height width;
   model sales=height width height*width;
   output out=a2 r=resid;
run;

proc rank data=a2 
   out=a3 normal=blom;
  var resid;
  ranks zresid;
proc sort data=a3; 
   by zresid;
symbol1 v=circle i=sm70;
proc gplot data=a3;
   plot resid*zresid/frame;
run;

proc glm data=a1;
   class height width;
   model sales=height width;
   means height / tukey lines;
run;