/*Ott 'n' Longnecker problem 9-17*/
options ls = 80;
data plant;
input nutrient$ nct;
lines;
A	19.4
A	32.6
A	27
A	32.1
A	33
B	18.2
B	24.6
B	25.5
B	19.4
B	21.7
B	20.8
C	20.7
C	21
C	20.5
C	18.8
C	18.6
C	20.1
C	21.3
;
proc univariate data = plant normal plot;
var nct;
by nutrient;
run;

proc glm data = plant;
class nutrient;
model nct = nutrient;
means nutrient/ duncan scheffe snk lsd tukey alpha = .05;
contrast '3DOK5 takes on all comers' nutrient -.5 1 -.5;
run;

quit;