data a; input a b y; datalines; 1 1 3.5 1 1 3.4 1 2 3.2 1 2 3.3 1 2 4.1 1 3 2.6 2 1 4.5 2 1 4.2 2 2 3.9 2 2 4.6 2 3 5.0 2 3 6.1 proc glm; class a b; model y=a b/e1 e2 e3; run; proc glm; class a b; model y=a|b/e1 e2 e3; run;