yydy653.sas

高老师的经典多元统计分析教材源代码非常经典使用

/*  yydy653.sas  */

data d653;
  input group $ x1-x6;
  cards;
北京   190.33  43.77   9.73  60.54  49.01  9.04
天津   135.20  36.40  10.47  44.16  36.49  3.94
河北    95.21  22.83   9.30  22.44  22.81  2.80
山西   104.78  25.11   6.40   9.89  18.17  3.25
内蒙   128.41  27.63   8.94  12.58  23.99  3.27
辽宁   145.68  32.83  17.79  27.29  39.09  3.47
吉林   159.37  33.38  18.37  11.81  25.29  5.22
黑龙江 116.22  29.57  13.24  13.76  21.75  6.04
上海   221.11  38.64  12.53 115.65  50.82  5.89
江苏   144.98  29.12  11.67  42.60  27.30  5.74
浙江   169.92  32.75  12.72  47.12  34.35  5.00
安徽   153.11  23.09  15.62  23.54  18.18  6.39
福建   144.92  21.26  16.96  19.52  21.75  6.73
江西   140.54  21.50  17.64  19.19  15.97  4.94
山东   115.84  30.26  12.20  33.61  33.77  3.85
河南   101.18  23.26   8.46  20.20  20.50  4.30
;

proc standard mean=0 std=1 data=d653  out=sta653;
    var x1-x6;
run;

proc fastclus data=sta653 out=out653 maxc=5 list;
   var x1-x6;
   id group;
run;
proc sort data=out653;
   by cluster;
run;
proc print data=out653;
  var group cluster DISTANCE ;
run;

proc candisc  data=out653 out=can653;
    var x1-x6;
    class cluster;
run;
proc plot data=can653;
  plot can2*can1=cluster;
run;