OPTIONS LS=72;
TITLE 'EXAMPLE 22. HIERARCHICAL MODEL: PURITY DATA';
DATA PURITY;
N=_N_;
INPUT SUPPLIER $ BATCH $ CPURITY @@;
CARDS;
S1 B1  1 S1 B1 -1 S1 B1  0 S1 B2 -2 S1 B2 -3 S1 B2 -4 
S1 B3 -1 S1 B3  0 S1 B3  1 S1 B4  1 S1 B4  4 S1 B4  0
S2 B1  1 S2 B1 -2 S2 B1 -3 S2 B2  0 S2 B2  4 S2 B2  2 
S2 B3 -1 S2 B3  0 S2 B3 -2 S2 B4  0 S2 B4  3 S2 B4  2
S3 B1  2 S3 B1  4 S3 B1  0 S3 B2 -2 S3 B2  0 S3 B2  2 
S3 B3  1 S3 B3 -1 S3 B3  2 S3 B4  3 S3 B4  2 S3 B4  1
;
PROC GLM;
CLASS SUPPLIER BATCH;
MODEL CPURITY=SUPPLIER BATCH(SUPPLIER);
RUN;
RANDOM BATCH(SUPPLIER)/TEST;
RUN;
TEST H=SUPPLIER E=BATCH(SUPPLIER);
RUN;