data serum;
 input y a b;
 cards;
143 1 1
169 1 1
95  1 1
111 1 1
132 1 1
150 1 1
141 1 1
152 2 1
83  2 1
91  2 1
86  2 1 
150 2 1
108 2 1
75  2 1
179 1 2
160 1 2
87  1 2
115 1 2
171 1 2
146 1 2
141 2 2
132 2 2
201 2 2
242 2 2
209 2 2
114 2 2
76  1 3
40  1 3
119 1 3
72  1 3
163 1 3
78  1 3
119 2 3
104 2 3
125 2 3
147 2 3
200 2 3
178 2 3
;
run;

proc glm data = serum;
  class a b;
  model y = a b a*b;
  lsmeans a b a*b;
  means a b a*b / duncan lines;  
title 'Example of Two-Way ANOVA with Interaction';
title2 'Data on Serum Fluorescence Ratings P. 440 of Fisher et al.';
run;