/* unbalanced split plot MJ 28.1 */
options nocenter ps=64 ls=76;
data split;
   infile 'usplit.dat';
   input herb plot fert height;
proc glm;
   class herb plot fert;
   model height = herb plot(herb) fert fert*herb;
   random plot(herb);
   test h = herb e = plot(herb);
   lsmeans herb / pdiff stderr e = plot(herb);
   lsmeans fert herb*fert / pdiff stderr;
   estimate 'nfert 1' intercept 180 herb 60 60 60
      fert 180 0 0 0 plot(herb) 20 20 20 12 12 12 12 12 15 15 15 15
      herb*fert 60 0 0 0 60 0 0 0 60 0 0 0 / divisor=180;