options nocenter ps=64 ls=76;   /* split plot MJ 24.1 */
data split;
   input fert n1 n2 n3 n4;
   cards;
1 35.4 37.9 41.6 40.3
2 36.7 38.2 42.7 41.6
3 34.8 36.4 43.6 42.8
4 39.5 40   44.5 47.6
;
data b; set split;
   yield = n1; block=1; var=1; output;
   yield = n2; block=1; var=2; output;
   yield = n3; block=2; var=1; output;
   yield = n4; block=2; var=2; output;
   drop n1--n4;

proc glm;
   class block var fert;
   model yield = block fert block*fert var var*fert;
   random block block*fert;
   test h = block fert e = block*fert;
   lsmeans fert / pdiff stderr e = block*fert;
   lsmeans var var*fert / pdiff stderr;

proc mixed;
   class block var fert;
   model yield = fert var var*fert;
   random block block*fert;
   lsmeans fert | var;