options linesize=80 pagesize=65 nodate;

title1 'Completely Randomized Design - Balanced Data';
title2 'Nitrogen Content of Red Clover Plants';
title3 'Steel & Torrie Text 2nd ed. - pg. 140';

data clover;
  length strain $ 9;
  do strain = '3DOk1', '3DOk5', '3DOk4', '3DOk7', '3DOk13', 'Composite';
    do rep = 1 to 5;
      input mg_n @@;
      output;
    end;
  end;
cards;
19.4 32.6 27.0 32.1 33.0
17.7 24.8 27.9 25.2 24.3
17.0 19.4  9.1 11.9 15.8
20.7 21.0 20.5 18.8 18.6
14.3 14.4 11.8 11.6 14.2
17.3 19.4 19.1 16.9 20.8
;;;;
proc print;
proc glm;
  classes strain;
  model mg_n = strain;
  means strain / lsd lines;

proc sort;
  by strain;
proc means;
  by strain;
  output out=aves mean=ave_n;
  var mg_n;
proc plot;
  plot ave_n*strain;