options nodate nonumber ps=55 ls=76;

title 'Kuehl -- Example 7.3 -- Glucose Standards';

data one;
input day run glucose @@;
datalines;
1 1 42.5 1 1 43.3 1 1 42.9 1 2 42.2 1 2 41.4 1 2 41.8
2 3 48.0 2 3 44.6 2 3 43.7 2 4 42.0 2 4 42.8 2 4 42.8
3 5 41.7 3 5 43.4 3 5 42.5 3 6 40.6 3 6 41.8 3 6 41.8
;
run;

proc glm;
class day run;
model glucose=day run(day);
test h=day e=run(day);
random day run(day);
run;

proc varcomp method=type1;
class day run;
model glucose=day run(day);
run;

quit;