yyexa809.sas

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

/*  YYexa809.sas  */

data data809;
  input x1-x4;
  n=_n_;
  cards;
1.611    10.59    0.69    1.67
1.429     9.44    0.61    1.50
1.447     5.97    0.24    1.25
1.572    10.72    0.75    1.71
1.483    10.99    0.75    1.44
1.371     6.46    0.41    1.31
1.665    10.51    0.53    1.52
1.403     6.11    0.17    1.32
2.620    21.51    1.40    2.59
2.033    24.15    1.80    1.89
2.015    26.86    1.93    2.02
1.501     9.74    0.87    1.48
1.578    14.52    1.12    1.47
1.735    14.64    1.21    1.91
1.453    12.88    0.87    1.52
1.765    17.94    0.89    1.40
1.532    29.42    2.52    1.80
1.488     9.23    0.81    1.45
2.586    16.07    0.82    1.83
1.992    21.63    1.01    1.89
;
proc print data=data809;
run;


* 8-9(1) ;

proc factor data=data809 p=0.8 method=prin simple;
  var x1-x4;
run;

proc factor data=data809 rotate=varimax  n=2 ;
  var x1-x4;
run;
proc factor data=data809 score n=2 out=o89;
   var x1-x4;
run;

* 7-8 (2) ;
proc sort data=o89;
  by factor1;
run;

proc print data=o89;
   var n factor1 factor2 x1-x4;
run;
proc plot data=o89;
   plot factor2*factor1 $ n='*'/href=0 vref=0;
run;

proc cluster data=data809   method=ave outtree=oc89;
  var x1-x4;
run;
proc tree data=oc89 horizontal graphics;
run;
proc cluster data=data809 method=ward outtree=oc89;
  var x1-x4;
run;
proc tree data=oc89 horizontal graphics;
run;

quit;