yyexa610.sas

来自「高老师的经典多元统计分析教材源代码非常经典使用」· SAS 代码 · 共 84 行

/*  yyexa610.sas  */


data data610;
  input x1-x7;
  cards;
0.05798    5.5150    347.10    21.910     8586    1742     61.69
0.08441    3.9700    347.20    19.710     7947    2000   2440.00
0.07217    1.1530     54.85     3.052     3860    1445   9497.00
0.15010    1.7020    307.50    15.030    12290    1461   6380.00
5.74400    2.8540    229.60     9.657     8099    1266   12520.0
0.21300    0.7058    240.30    13.910     8980    2820   4135.00
;
proc print;
run;


* 6-10(1) ;

proc cluster data=data610 method=ave outtree=o610;
  var x1-x7;
run;
proc tree data=o610 horizontal graphics;
run;
proc cluster data=data610 method=ward outtree=o610;
  var x1-x7;
run;
proc tree data=o610 horizontal graphics;
run;

proc corr data=data610 outp=oc610;
   var x1-x7;
run;

data ob610(drop=i);
  set oc610(firstobs=4);
  array  x(7) x1-x7;
  do i=1 to 7;
    x(i)=1-abs(x(i));
  end;
run;
proc print data=ob610;
run;
proc cluster data=ob610 method=ave outtree=o610;
  var x1-x7;
run;
proc tree data=o610 horizontal graphics;
run;
proc cluster data=ob610 method=ward outtree=o610;
  var x1-x7;
run;
proc tree data=o610 horizontal graphics;
run;


data od610(drop=i);
  set oc610(firstobs=4);
  array  x(7) x1-x7;
  do i=1 to 7;
    x(i)=sqrt(1-abs(x(i)*x(i)));
  end;
run;
proc print data=od610;
run;
proc cluster data=od610 method=ave outtree=o610;
  var x1-x7;
run;
proc tree data=o610 horizontal graphics;
run;
proc cluster data=od610 method=ward outtree=o610;
  var x1-x7;
run;
proc tree data=o610 horizontal graphics;
run;


* 6-10 (2) ;
proc varclus data=data610;
run;
proc varclus data=data610 maxc=3 summary outtree=tree;
run;
proc tree data=tree horizontal graphics;
run;