DM 'LOG;CLEAR;OUT;CLEAR;';
options ls=79 pagesize=60 nodate nonumber;

DATA Raw;
  INFILE 'flea_beetle.dat';
  INPUT species y1-y4;

proc discrim data=raw method=normal manova crossvalidate out=norm_out;
class species;
priors prop;
var y1-y4;

proc print data = norm_out;


proc discrim data=raw method=normal pool=test slpool=.1 crossvalidate;
class species;
priors prop;
var y1-y4;



proc discrim data = raw distance can out=outcan;
class species;
var y1-y4;

proc print data = outcan;


goptions reset=global gunit=pct border ftext=swissb htitle=6 htext=3
         cback=white ;
title "Box Plots of Canonical Variable";
title2 h=4 'Flea Data: Table 3.8 in Rencher';
symbol interpol=boxt10 bwidth=6 bwidth=6 value=circle height=2;
axis1 label=none major=(number=2) value = ('Species HO' 'Species HC')
   offset=(5,5) length=50;
proc gplot data = outcan;
     plot can1*species /haxis=axis1;



proc discrim data=raw method=npar k=5 crossvalidate out=out1;
class species;
priors prop;
var y1-y4;

proc print data = out1;
run;