OPTIONS LS=72;
TITLE 'EXAMPLE 19. UNBALANCED TWO-WAY ANOVA: CATALYST DATA';
DATA REACTION;
DO CATALYST=1 TO 4;
DO BLOCK=1 TO 4;
INPUT TIME @@;
OUTPUT;
END;
END;
CARDS;
73 74  . 71
 . 75 67 72
73 75 68  .
75  . 72 75
;
PROC PRINT DATA=REACTION;
PROC GLM DATA=REACTION;
CLASS BLOCK CATALYST;
MODEL TIME=BLOCK CATALYST;

LSMEANS CATALYST/STDERR;
ESTIMATE 'C1VSC2' CATALYST 1 -1 0 0;
ESTIMATE 'C1VSC3' CATALYST 1 0 -1 0;
ESTIMATE 'C1VSC4' CATALYST 1 0 0 -1;
ESTIMATE 'C2VSC3' CATALYST 0 1 -1 0;
ESTIMATE 'C2VSC4' CATALYST 0 1 0 -1;
ESTIMATE 'C3VSC4' CATALYST 0 0 1 -1;