* Myers and Montgomery, Example 7.6, Table 7.6, p. 332;

data mm_eg76;
 input x1 x2 block y;
 cblock=2*(block-1.5);
cards;
-1     -1     1 80.5
-1      1     1 81.5
 1     -1     1 82.0
 1      1     1 83.5
 0      0     1 83.9
 0      0     1 84.3
 0      0     1 84.0
 0      0     2 79.7
 0      0     2 79.8
 0      0     2 79.5
 1.414  0     2 78.4
-1.414  0     2 75.6
 0      1.414 2 78.5
 0     -1.414 2 77.0
;
run;

* Analysis without blocks;
proc rsreg data=mm_eg76;
 model y=x1 x2/nocode;
run;

* Analysis with blocks as in data;
proc rsreg data=mm_eg76;
 model y=block x1 x2/nocode covar=1;
run;

* Analysis with coded blocks;
proc sort data=mm_eg76 out=s_mm_eg76;
 by cblock x1 x2;
run;
proc rsreg data=s_mm_eg76;
 model y=cblock x1 x2/nocode covar=1 lackfit;
run;