yydy633.sas

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

/*   yydy633.sas  */

title;
 goptions  ftext="宋体";
data d633;
   input num group $ x1-x3;
   cards;
  1    北京       9.30     30.55      8.70
  2    天津       4.67     29.38      8.92
  3    河北       0.96     24.69     15.21
  4    山西       1.38     29.24     11.30
  5    内蒙       1.48     25.47     15.39
  6    辽宁       2.60     32.32      8.81
  7    吉林       2.15     26.31     10.49
  8    黑龙       2.14     28.46     10.87
  9    上海       6.53     31.59     11.04
 10    江苏       1.47     26.43     17.23
 11    浙江       1.17     23.74     17.46
 12    安徽       0.88     19.97     24.43
 13    福建       1.23     16.87     15.63
 14    江西       0.99     18.84     16.22
 15    山东       0.98     25.18     16.87
 16    河南       0.85     26.55     16.15
 17    湖北       1.57     23.16     15.79
 18    湖南       1.14     22.57     12.10
 19    广东       1.34     23.04     10.45
 20    广西       0.79     19.14     10.61
 21    海南       1.24     22.53     13.97
 22    西藏       0.96     21.65     16.24
 23    贵州       0.78     14.65     24.27
 24    云南       0.81     13.85     25.44
 25    西藏       0.57      3.85     44.43
 26    陕西       1.67     24.36     17.62
 27    甘肃       1.10     16.85     27.93
 28    青海       1.49     17.76     27.70
 29    宁夏       1.61     20.27     22.06
 30    新疆       1.85     20.66     12.75
 ;

proc transpose data=d633 out=dd633 prefix=y;
  var x1-x3;
run;
proc print data=dd633;

proc corr data=dd633 out=corr633;
   var y1-y30;
run;

data dc633(type=distance);
  set corr633;
  array y{30} y1-y30;
  if _type_='CORR';
  nn=_n_-3;
  do i=1 to 30;
    y{i}=1-y{i}*y{i};
    y{i}=sqrt(y{i});
  end;
  drop i;
run;

proc print;
run;
proc cluster data=dc633 method=com nonorm outtree=b633a;
   var y1-y30 ; copy nn;
run;
proc tree data=b633a  horizontal graphics ;
  title 'complete--例6.3.3';
  id nn;
run;
title;

proc cluster data=dc633 method=cen  nonorm outtree=b633b;
   var y1-y30; copy nn;
run;
proc tree data=b633b  horizontal graphics ;
  title 'centroid--例6.3.3'; id nn;
run;


proc cluster data=dc633 method=ward   nonorm outtree=b633c;
   var y1-y30; copy nn;
run;
proc tree data=b633c  horizontal graphics n=4 out=c633 ;
  title 'Ward--例6.3.3'; id nn;
run;

proc sort data=c633;
  by cluster;
run;

proc  print data=c633;
run;