cluster1.sas

SAS for Windows统计分析系统教程,里面有源码

options nodate nonumber;
filename cluster1 'd:\spss9data\body1.dbf';
proc dbf db3=cluster1 out=c;
label x1='Body Height'
      x2='Sitting Height'
      x3='Body Weight'
      x4='Chest Circumference'
      x5='Shoulder Width'
      x6='Pelvis Width';
/* Code:
   OB1: Beijing Shi              OB15: Shanghai Shi
   OB2: Tianjin Shi              OB16: Jiangsu Sheng
   OB3: Hebei Sheng              OB17: Zhejiang Sheng
   OB4: Shanxi Sheng             OB18: Anhui Sheng
   OB5: Nei Mongol Zizhiqu       OB19: Fujian Sheng
   OB6: Liaoning Sheng           OB20: Jiangxi Sheng
   OB7: Jilin Sheng              OB21: Henan Sheng
   OB8: Heilongjiang Sheng       OB22: Hubei Sheng
   OB9: Shandong Sheng           OB23: Hunan Sheng
   OB10:Shaanxi Sheng            OB24: Guangdong Sheng
   OB11:Gansu Sheng              OB25: Guangxi Zhuangzu Zizhiqu
   OB12:Qinghai Sheng            OB26: Sichuan Sheng
   OB13:Ningxia Huizu Zizhiqu    OB27: Guizhou Sheng
   OB14:Xinjiang Uygur Zizhiqu   OB28: Yunnan Sheng*/;
title '1. Cluster Analysis, Average Method';
proc cluster method=average pseudo data=c;
    var x1-x6;
run;
title '2. Cluster Analysis, Centroid Method';
proc cluster method=centroid pseudo data=c;
    var x1-x6;
run;
title '3. Cluster Analysis, Complete Method';
proc cluster method=complete data=c;
    var x1-x6;
run;
title '4. Cluster Analysis, Density Method';
proc cluster method=density k=3 data=c;
    var x1-x6;
run;
title '5. Cluster Analysis, Eml Method';
proc cluster method=eml data=c;
    var x1-x6;
run;
title '6. Cluster Analysis, Flexible Method';
proc cluster method=flexible data=c;
    var x1-x6;
run;
title '7. Cluster Analysis, Mcquitty Method';
proc cluster method=mcquitty data=c;
    var x1-x6;
run;
title '8. Cluster Analysis, Median Method';
proc cluster method=median data=c;
    var x1-x6;
run;
title '9. Cluster Analysis, Single Method';
proc cluster method=single data=c;
    var x1-x6;
run;
title '10. Cluster Analysis, Twostage Method';
proc cluster method=twostage k=3 data=c;
    var x1-x6;
run;
title '11. Cluster Analysis, Ward Method';
proc cluster method=Ward pseudo data=c;
    var x1-x6;
run;