*NKNW677.sas, one way anova using data in Table 16.1;
options nocenter;
data a1; 
   infile 'u:\.www\datasets512\CH16TA01.DAT';
   input cases design store;

proc print data=a1; 
run;

proc glm data=a1; 
   class design; 
   model cases=design;
   means design;
run;

symbol1 v=circle i=none;
proc gplot data=a1; 
   plot cases*design/frame;
run;

proc means data=a1; 
   var cases; by design;
   output out=a2 mean=avcases;
proc print data=a2;

symbol1 v=circle i=join;
proc gplot data=a2; 
   plot avcases*design/frame;
run;