options nocenter;
data sim;
   mu=2;
   do rep=1 to 5;
   r=rannor(528911);
   do block=1 to 3;
   b=.5* rannor(63281004);
   do unit=1 to 2;
   input trt @@;
   if trt=4 then t=-.6;
   if trt=3 then t=-.5;
   if trt=5 then t=-.4;
   if trt=1 then t=0;
   if trt=2 then t=.5;
   if trt=6 then t=1.0;
   epsilon=.7071*rannor(20981323);
   y=mu+r+b+t+epsilon;
   output;
   end;
   end;
   end;
   cards;
1  2
3  4
5  6
1  3
2  5
4  6
1  4
2  6
3  5
1  5
2  4
3  6
1  6
2  3
4  5
proc sort data=sim;
   by rep block trt ;
data sim;
   set sim;
   y= round(y,.1);
proc print data=sim;
   var rep block trt y;

Simulated Data. Compare 6 treatments
Balanced incomplete block design, resolvable.
Reps and blocks(rep) are random effects

Obs    rep    block    trt     y

  1     1       1       1     1.0
  2     1       1       2     1.5
  3     1       2       3     1.5
  4     1       2       4     0.5
  5     1       3       5     0.6
  6     1       3       6     2.4
  7     2       1       1     2.9
  8     2       1       3     1.0
  9     2       2       2     2.2
 10     2       2       5     1.5
 11     2       3       4     1.5
 12     2       3       6     2.6
 13     3       1       1     2.9
 14     3       1       4     3.0
 15     3       2       2     4.1
 16     3       2       6     4.3
 17     3       3       3     3.8
 18     3       3       5     3.5
 19     4       1       1     1.4
 20     4       1       5     2.1
 21     4       2       2     0.4
 22     4       2       4     -0.2
 23     4       3       3      1.7
 24     4       3       6      2.5
 25     5       1       1      0.6
 26     5       1       6      2.1
 27     5       2       2      2.5
 28     5       2       3      2.3
 29     5       3       4      2.3
 30     5       3       5      2.6