options ls=78;
title "Two-Way MANOVA: Rice Data";

data rice;
  infile "D:\Statistics\STAT 505\data\rice.csv" firstobs=2 delimiter=',';
  input block variety $ height tillers;
  run;

 /*
  * The class statement specifies block and variety as categorical variables.
  * The model statement specifies responses height and tillers with predictors
  * block and variety. lsmeans displays the least-squares means for variety.
  * The manova statement requests the test for response mean vectors across
  * levels of variety, and the printe and printh options display sums of 
  * squares and cross products for error and the hypothesis, respectively.
  */

proc glm data=rice;
  class block variety;
  model height tillers=block variety;
  lsmeans variety;
  manova h=variety / printe printh;
  run;