options nodate nonumber ps=55 ls=80;

title 'Kuehl -- Example 2.1 Bacteria/Meat';

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 print;
run;

proc means mean std;
class package;
var logcount;
run;

proc glm;
class package;
model logcount=package / p;
run;

quit;