*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;

*95% CI with Bonferonni adjustment;
proc glm data=a1; 
   class design; 
   model cases=design;
   means design/t bon clm; 
run;

*Comparison of means using tukey adjustment;
proc glm data=a1; 
   class design; 
   model cases=design;
   means design/tukey lines; 
run;

*Constrast and estimate statements;
proc glm data=a1; 
   class design; 
   model cases=design;
   contrast '1&2 v 3&4' design .5 .5 -.5 -.5; 
   estimate '1&2 v 3&4' design .5 .5 -.5 -.5; 
run;

*Joint test of several contrasts;
proc glm data=a1; 
   class design; 
   model cases=design;
   contrast '1 v 2&3&4' design 1 -.3333 -.3333 -.3333; 
   estimate '1 v 2&3&4' design 3 -1 -1 -1 /divisor=3;
   contrast '2 v 3 v 4' design 0 1 -1  0,
                        design 0 0  1 -1; 
run;