proc format;
         value specname
            1='Setosa    '
            2='Versicolor'
            3='Virginica ';
      run;
   
      data iris;
         title 'Fisher (1936) Iris Data';
         input SepalLength SepalWidth PetalLength PetalWidth Species @@;
         format Species specname.;
         label SepalLength='Sepal Length in mm.'
               SepalWidth ='Sepal Width in mm.'
               PetalLength='Petal Length in mm.'
               PetalWidth ='Petal Width in mm.';
         symbol = put(species, specname10.);
         datalines;
      50 33 14 02 1 64 28 56 22 3 65 28 46 15 2 67 31 56 24 3
      63 28 51 15 3 46 34 14 03 1 69 31 51 23 3 62 22 45 15 2
      59 32 48 18 2 46 36 10 02 1 61 30 46 14 2 60 27 51 16 2
      65 30 52 20 3 56 25 39 11 2 65 30 55 18 3 58 27 51 19 3
      68 32 59 23 3 51 33 17 05 1 57 28 45 13 2 62 34 54 23 3
   
   ... more lines ...   

      55 23 40 13 2 66 30 44 14 2 68 28 48 14 2 54 34 17 02 1
      51 37 15 04 1 52 35 15 02 1 58 28 51 24 3 67 30 50 17 2
      63 33 60 25 3 53 37 15 02 1
      ;

   proc fastclus data=iris maxc=2 maxiter=10 out=clus;
      var SepalLength SepalWidth PetalLength PetalWidth;
   run;
   
   proc freq;
      tables cluster*species;
   run;
   
   proc fastclus data=iris maxc=3 maxiter=10 out=clus;
      var SepalLength SepalWidth PetalLength PetalWidth;
   run;
   
   proc freq;
      tables cluster*Species;
   run;
   
   proc candisc anova out=can;
      class cluster;
      var SepalLength SepalWidth PetalLength PetalWidth;
      title2 'Canonical Discriminant Analysis of Iris Clusters';
   run;
   
   proc sgplot data=Can;
      scatter y=Can2 x=Can1 /group=Cluster ;
      title2 'Plot of Canonical Variables Identified by Cluster';
   run;