yyexa711.sas

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

/*  yyexa711.sas  */

data data711;
   input name : $18.  x1-x8;
   cards;
1(冶金)          90342  52455 101091 19272  82.000 16.100 197435 0.172
2(电力)           4903   1973   2035 10313  34.200  7.100 592077 0.003
3(煤炭)           6735  21139   3767  1780  36.100  8.200 726396 0.003
4(化学)          49454  36241  81557 22504  98.100 25.900 348226 0.985
5(机械)         139190 203505 215898 10609  93.200 12.600 139572 0.628
6(建材)          12215  16219  10351  6382  62.500  8.700 145818 0.066
7(森工)           2372   6572   8103 12329 184.400 22.200  20921 0.152
8(食品)          11062  23078  54935 23804 370.400 41.000  65486 0.263
9(纺织)          17111  23907  52108 21796 221.500 21.500  63806 0.276
10(缝纫)          1206   3930   6126 15586 330.400 29.500   1840 0.437
11(皮革)          2150   5704   6200 10870 184.200 12.000   8913 0.274
12(造纸)          5251   6155  10383 16875 146.400 27.500  78796 0.151
13(文教艺术用品)  14341  13203  19396 14691  94.600 17.800   6354 1.574
 ;
proc print data=data711;
run;

* 7-11(1) ;
proc princomp data=data711  prefix=z out=o711;
  var x1-x8;
run;

* 7-11(2) ;
proc sort data=o711;
  by z1;
run;
proc print data=o711;
  var name z1 z2 x1-x8;
run;

proc plot data=o711;
  plot z2*z1 $ name='*';
run;

proc cluster data=o711 method=ave outtree=b711;
  var z1-z3;
run;
proc tree data=b711 horizontal graphics;
run;