yydy632.sas

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

/* yydy632.sas */

title;
 goptions  ftext="宋体";

 data d632;
   input num @6 Country $8. x1-x6;
   cards;
  1  美国          631.60    161.90   0.36     403     26073   35.34
  2  日本          498.40    143.20   3.57     176     10223    6.26
  3  德国          557.60     70.60   2.18     199     11571    9.48
  4  瑞典          684.10    281.80   1.40     286     16660   29.39
  5  瑞士          644.00     93.50   1.98     234     13621   22.68
  6  丹麦          620.30    248.60   2.56     296     17210   21.84
  7  新加坡        498.40    147.50   2.50     284     13578   13.49
  8  中国台湾      469.40     56.10   3.68     119      6911    1.72
  9  韩国          434.50     73.00   3.36      99      5795    1.66
 10  巴西           81.90     16.30   3.02      19       876    0.52
 11  智利          138.60      8.20   1.40      31      1411    1.28
 12  墨西哥         92.20      9.80   2.61      31      1751    0.35
 13  俄罗斯        174.90      5.00   5.12      24      1101    0.48
 14  波兰          169.00      6.50   3.68      40      1796    1.45
 15  匈牙利        262.20     49.40   2.66      68      3067    3.09
 16  马来西亚      195.50     88.40   4.19      53      2734    1.25
 17  泰国           78.60     27.80   4.95      22      1662    0.11
 18  印度           13.60      0.30   6.28       2       101    0.01
 19  法国          559.10     42.90   1.27     201     11702    4.76
 20  英国          521.10    122.50   0.98     248     14461   11.91
 ;


proc cluster data=d632 method=com std  nonorm pseudo ccc
             outtree=b632;
   var x1-x6;
   id num;
proc tree data=b632  horizontal graphics ;
  title   'Average--例6.3.2';
run;
title;

proc cluster data=d632 method=cen std  nonorm  pseudo ccc
             outtree=b632;
   var x1-x6;
   id num;
proc tree data=b632  horizontal graphics ;
  title  'centroid--例6.3.2';
run;


proc cluster data=d632 method=ward std  nonorm  pseudo ccc
             outtree=b632;
   var x1-x6;
   id num;
proc tree data=b632  horizontal graphics n=5 out=c632 ;
  copy num x1-x6;
  title 'Ward--例6.3.2';
run;
title;

proc sort data=c632;
  by cluster;
run;

proc  print data=c632;
    var cluster num x1-x6;
run;
proc means data=c632 ;
  by cluster;
  var x1-x6;
run;