options pagesize=50;
data;
input quality season epith @@;
cards;
1 1 24.2 2 1 23.2 3 1 19.8
1 2 16.5 2 2 13.2 3 2 10.4
1 3 17.9 2 3 12.5 3 3 10.1
1 4 17.8 2 4 16.3 3 4 15.1
proc anova;
class quality season;
model epith=quality season;
means quality/bon tukey scheffe alpha=.05;
run;