options nodate nonumber ls=76 ps=55;

title 'Kuehl -- Example 5.1 -- Casting/Strength';

data casting;
input cast strength @@;
datalines;
1 88.0 1 88.0 1 94.8 1 90.0 1 93.0 1 89.0 1 86.0 1 92.9 1 89.0 1 93.0
2 85.9 2 88.6 2 90.0 2 87.1 2 85.6 2 86.0 2 91.0 2 89.6 2 93.0 2 87.5
3 94.2 3 91.5 3 92.0 3 96.5 3 95.6 3 93.8 3 92.5 3 93.2 3 96.2 3 92.5
;

proc means mean std;
class cast;
var strength;
run;

proc glm;
class cast;
model strength=cast;
random cast;
run;

proc varcomp method=type1;
class cast;
model strength=cast;
run;

quit;