options nodate nonumber ps=55 ls=80;

title 'Kuehl -- Example 2.1 (Extended) Bacteria/Meat p. 77';

data bacteria;
input package logcount @@;
datalines;
1 7.66 1 6.98 1 7.80
2 5.26 2 5.44 2 5.80
3 7.41 3 7.33 3 7.04
4 3.51 4 2.91 4 3.66
;
proc glm;
class package;
model logcount=package;
estimate '1 vs others   ' package 3 -1 -1 -1 / divisor=3;
estimate '2 vs (3 & 4)  ' package 0 2 -1 -1/ divisor=2;
estimate '3 vs 4        ' package 0 0 1 -1;
run;

proc glm;
class package;
model logcount=package;
contrast '1 vs others   ' package 3 -1 -1 -1;
contrast '2 vs (3 & 4)  ' package 0 2 -1 -1;
contrast '3 vs 4        ' package 0 0 1 -1;
run;

quit;