DEFINE:
          stb=CLUSTER_MEAN(st);
           CENTER st mav (GRANDMEAN);

  ANALYSIS: type = twolevel;

  MODEL:
          %WITHIN%
          int on mav
          st(a);
          %BETWEEN%
          int on stb(b);

          MODEL CONSTRAINT:
          NEW(effect);
          effect=b-a;