OPTION LS=76 NOCENTER;
DATA rat;INPUT drug $ x2 y1-y10;
ARRAY in (10) y1-y10;
do i=1 to 10;
  y=in(i);
x1=0;if drug='B' then x1=1;
combined=x1+x2;output;
end;
CARDS;
A 5  15 16 16 15 18 16 13 17 19 15  
A 10 18 16 17 15 18 19 19 18 20 16 
A 20 20 17 19 18 21 21 18 20 19 17 
B 5  16 15 17 15 18 18 17 17 15 16 
B 10 19 18 21 20 22 21 23 22 20 19 
B 20 24 23 25 24 23 22 25 26 25 24 
;
PROC PRINT;
TITLE 'Print all the data';
PROC GLM;MODEL y=x1 x2 x1*x2/P CLI;
TITLE 'For full model y=x1 x2 x1*x2'; 
 
PROC GLM;MODEL y=x1 x2;
TITLE 'For partial model y=x1 x2';
 
PROC GLM;CLASS combined;
MODEL y=combined;
TITLE 'For pure error estimation';