options nocenter ps=64 ls=76;
data a;
   infile 'unbal.dat' firstobs=2;
   input T B y;
/* sequential fit of models for estimation of effects */
proc glm;
   class T B;
   model y = T / ss1 solution;
proc glm;
   class T B;
   model y = B / ss1 solution;
proc glm;
   class T B;
   model y = T B / ss1 solution;
proc glm;
   class T B;
   model y = T B T*B / ss1 solution;
/* examine types of sums of squares and LS vs marginal means */
proc glm;
   class T B;
   model y = T | B / e1 e2 e3;
   lsmeans T B T*B / stderr pdiff;
   means T B T*B / lsd lines;