/* x6_5.sas */
/* Purpose:  */

title1 'Example 6.5.';

data a;
input plant$ dry1-dry6 /
      plant$ nit1-nit6 /
      plant$ gre1-gre6;
  block=1; dry=dry1; nit=nit1; gre=gre1; output;
  block=2; dry=dry2; nit=nit2; gre=gre2; output;
  block=3; dry=dry3; nit=nit3; gre=gre3; output;
  block=4; dry=dry4; nit=nit4; gre=gre4; output;
  block=5; dry=dry5; nit=nit5; gre=gre5; output;
  block=6; dry=dry6; nit=nit6; gre=gre6; output;
keep plant block dry nit gre;
cards;
marchigi   9.170  10.683  10.063   8.104  10.018   9.570
marchigi   0.286   0.335   0.315   0.259   0.319   0.304
marchigi  40.959  47.677  44.950  36.919  45.859  43.838
kayseri    9.403  10.914  11.018  11.385  13.387  12.848
kayseri    0.285   0.330   0.333   0.339   0.400   0.383
kayseri   42.475  49.546  50.152  51.718  60.758  58.334
atlantic  11.349  10.971   9.794   8.944  11.715  11.903
atlantic   0.369   0.357   0.319   0.291   0.379   0.386
atlantic  52.475  50.757  45.151  42.221  55.505  56.364
;

proc print data=a;
title2 'List of data';

proc anova data=a;
  class plant block;
  model dry nit gre = plant block;
  manova h=plant /printh printe summary;
  manova h=block /printh printe summary;
title2 'Multivariate ANOVA.';

run;