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

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

Generate a plot of the means;
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;

*Create 95% CI for each factor level.
  Uses each factor level standard std deviation;
proc means data=a1 
   mean std stderr clm
   maxdec=2; 
   class design;
   var cases;
run;

*Create 95% CI for each factor level.
  Uses pooled variance;
proc glm data=a1; 
   class design; 
   model cases=design;
   means design/t clm; 
run;

*Create 95% CIs using bonferonni adjustment;
proc glm data=a1; 
   class design; 
   model cases=design;
   means design/bon clm; 
run;

*Comparison of means using various adjustments;
proc glm data=a1; 
   class design; 
   model cases=design;
   means design/lsd tukey bon scheffe; 
   means design/lines tukey; 
run;