options nodate nonumber ls=76 ps=55;

title 'Kuehl -- Example 6.2 -- Concrete/Strength';

data casting;
input aggtype $ compmeth $ strength @@;
datalines;
B S 68 B S 63 B S 65 B R 126 B R 128 B R 133 B L 93 B L 101 B L 98 B V 56 B V 59 B V 57
S S 71 S S 66 S S 66 S R 107 S R 110 S R 116 S L 63 S L 60 S L 59 S V 40 S V 41 S V 44
;

proc means mean std;
class aggtype compmeth;
var strength;
run;

proc glm;
class aggtype compmeth;
model strength=aggtype compmeth aggtype*compmeth;
estimate 'B vs S in Low      ' aggtype 1 -1 aggtype*compmeth 1 0 0 0 -1 0 0 0;
estimate 'B vs S in Regular  ' aggtype 1 -1 aggtype*compmeth 0 1 0 0 0 -1 0 0;
estimate 'B vs S in Static   ' aggtype 1 -1 aggtype*compmeth 0 0 1 0 0 0 -1 0;
estimate 'B vs S in Very low ' aggtype 1 -1 aggtype*compmeth 0 0 0 1 0 0 0 -1;
run;



quit;