options nocenter ps=50 ls=79;

data budding;
   infile 'budding.dat' firstobs=2;
   input locus days;

proc anova;
   class locus;
   model days = locus;
   means locus / lsd;

proc glm;
   class locus;
   model days = locus / ss1 solution;
   lsmeans locus / pdiff stderr;
   estimate 'grand mean' intercept 3 LOCUS 1 1 1 / divisor=3;