options ps=64 ls=80 nocenter;
data a;
   infile 'alfalfa.dat' firstobs=2;
   input pop $ sample rep seedwt harv1 harv2;

proc plot;
   plot harv2*harv1=sample;
proc glm;
   class pop sample rep;
   model seedwt harv1 harv2 = rep pop sample(pop);
   lsmeans pop / stderr pdiff e=sample(pop);
   lsmeans pop sample(pop) / stderr pdiff;
proc glm;
   class pop sample rep;
   model harv2 = harv1 rep pop sample(pop);

data b; set a;
   yield = harv1; harvest = 1; output;
   yield = harv2; harvest = 2; output;
   drop harv1 harv2;

proc plot; 
   plot yield*seedwt=sample;
proc glm;
   class pop sample rep harvest;
   model yield = rep pop sample(pop) rep*sample(pop) harvest harvest*pop;
   lsmeans pop / stderr pdiff e=sample(pop);
proc glm;
   class pop sample rep harvest;
   model yield = seedwt rep pop sample(pop) seedwt*rep*sample(pop)
      harvest harvest*pop;
   test h = pop e = sample(pop);
   test h = pop sample(pop) e = seedwt*rep*sample(pop);
proc sort; by pop sample rep;
proc means noprint; by pop sample rep;
   var yield seedwt;
   output out=mean mean=yield seedwt;
proc glm;
   class pop sample rep;
   model yield = seedwt rep pop sample(pop);
   lsmeans pop / stderr pdiff e=sample(pop);
   lsmeans pop sample(pop) / stderr pdiff;