/* x5_1.sas */
/* Purpose:  */

title1 'Example 5.1';

data a;
  input group$ rep strength;
cards;
A 1 3210
A 2 3000
A 3 3315
B 1 3225
B 2 3320
B 3 3165
B 4 3145
C 1 3220
C 2 3410
C 3 3320
C 4 3370
D 1 3545
D 2 3600
D 3 3580
D 4 3485
;
proc print;
title2 'List of data.';

proc varcomp method=type1;
  class group;
  model strength = group;

proc glm;
  class group;
  model strength = group /solution;
  means group;
  lsmeans group;

run;