OPTIONS NODATE NONUMBER PS=60 LS=80;

DATA in; 
    INPUT obs labno batch btu @@; OUTPUT;
CARDS;
 1  1  0  11655
 2  1  0  11690
 3  1  0  11654
 4  1  0  11653
 5  1  0  11665
 6  1  1  11754
 7  1  1  11676
 8  1  1  11712
 9  1  1  11711
10  1  1  11684
11  2  0  11665
12  2  0  11618
13  2  0  11572
14  2  0  11686
15  2  0  11725
16  2  1  11602
17  2  1  11582
18  2  1  11764
19  2  1  11618
20  2  1  11637
21  3  0  11586
22  3  0  11588
23  3  0  11599
24  3  0  11596
25  3  0  11622
26  3  1  11428
27  3  1  11429
28  3  1  11488
29  3  1  11632
30  3  1  11478
31  4  0  11651
32  4  0  11780
33  4  0  11663
34  4  0  11769
35  4  0  11800
36  4  1  11630
37  4  1  11587
38  4  1  11567
39  4  1  11531
40  4  1  11630
41  5  0  11591
42  5  0  11605
43  5  0  11549
44  5  0  11590
45  5  0  11609
46  5  1  11739
47  5  1  11747
48  5  1  11768
49  5  1  11768
50  5  1  11794
;

PROC GLM DATA=in;
     CLASS labno batch;
     MODEL btu = labno|batch / SS3;
     MEANS labno|batch;
     RANDOM batch labno*batch / TEST;
TITLE 'MIXED EFFECTS MODEL';

PROC VARCOMP DATA=in METHOD=REML;
     CLASS labno batch;
     MODEL btu = labno batch labno*batch / FIXED=1;
TITLE2 'VARIANCE COMPONENTS ANALYSIS';
RUN;