%newsas(2by3);
title 'Data from a 2 x 3 Two-Factor Design';
data Two;
     Input A $  B $ @;
     Do Subject = 1 to 8;
        Input Y @;
        Output;
        End;
datalines;
A1  B1   7  33  26  27  21   6  14  19
A1  B2   6  11  11  18  14  18  19  14
A1  B3   9  12   6  24   7  10   1  10
A2  B1  42  25   8  28  30  22  17  32
A2  B2  28   6   1  15   9  15   2  37
A2  B3  13  18  23   1   3   4   6   2
proc print data=Two;

%meanplot(data=Two, var=Y, class=B A);

/*--------------------------------------------*
 | Carry out the anova, with CI for  B  means |
 *--------------------------------------------*/
proc glm data=Two;
     class A B;
     model Y = A B A*B;
     means B / bon lines;
     title2 'Proc GLM results';
run;