options nodate nonumber ps=55 ls=76;

title 'Kuehl -- Example 8.1 -- Nitrate Nitrogen';

data one;
do trt=1 to 6;
do block=1 to 4;
input nitrate @@;
output;
end; end;
datalines;
34.98 41.22 36.94 39.97 40.89 46.69 46.65 41.90 42.07 49.42 52.68 42.91
37.18 45.85 40.23 39.20 37.99 41.99 37.61 40.45 34.89 50.15 44.57 43.29
;

proc glm;
class trt block;
model nitrate=trt block;
means trt / dunnett('4');
run;

quit;