options ps=54 ls=76;

data amoeba;
infile 'entozamoeba.dat';
input trt yield;
run;


proc glm;
class trt;
model yield = trt;
means trt / bon tukey hovtest;
output out=amoebaout r=amoebares;
run;

proc chart;
var amoebares;
run;

quit;