options nodate nonumber ps=55 ls=80;

title 'Kuehl -- Example 2.2 Phlebitis/Rabbits';

data rabbits;
input trt tempdiff @@;
datalines;
1 2.2 1 1.6 1 0.8 1 1.8 1 1.4 1 0.4 1 0.6 1 1.5 1 0.5
2 0.3 2 0.0 2 0.6 2 0.0 2 -0.3 2 0.2
3 0.1 3 0.1 3 0.2 3 -0.4 3 0.3 3 0.1 3 0.1 3 -0.5
;

proc print;
run;

proc means mean std;
class trt;
var tempdiff;
run;

proc glm;
class trt;
model tempdiff=trt / p;
run;

quit;