options nocenter ls=79 ps=100;

data soil;
   infile 'treeave.dat' firstobs=2;
   input soil$ rep score height;

proc anova;
   class soil;
   model score = soil / ss1;
   means soil; /* means and SDs by group */
   means soil / lsd bon scheffe regwf; /* based on F */
   means soil / tukey snk regwq; /* based on range */
run;