options ls=78;
title "Box's Test - Pottery Data";

data pottery;
  infile "D:\Statistics\STAT 505\data\pottery.csv" delimiter=',' firstobs=2;
  input site $ al fe mg ca na;
  run;

 /* The pool option specifies that covariance matrices should be pooled.
  * The class statement specifies the categorical variable site.
  * The var statement specifies the five variables to be used to 
  * calculate the covariance matrix.
  */

proc discrim data=pottery pool=test;
  class site;
  var al fe mg ca na;
  run;