options nodate nonumber ps=55 ls=76;

title 'Kuehl -- Example 14.2 -- Turfgrass Data';

data one;
do nitrogen=1 to 4;
do block=1 to 2;
do years=2 to 8 by 3;
input chlrfl @@;
output;
end; end; end;
datalines;
3.8 5.3 5.9 3.9 5.4 4.3 5.2 5.6 5.4 6.0 6.1 6.2
6.0 5.6 7.8 7.0 6.4 7.8 6.8 8.6 8.5 7.9 8.6 8.4
;

proc glm;
class block nitrogen years;
model chlrfl = block nitrogen block*nitrogen years years*nitrogen;
test h=nitrogen e=block*nitrogen;
run;


quit;