yydy731b.sas

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

/* yydy731.sas   */


data d731(type=corr);
  input _name_$  x1-x16;
  _type_='corr';
  cards;
x1 1.0 0.79 0.36 0.96 0.89 0.79 0.76 0.26 0.21 0.26 0.07 0.52 0.77 0.25 0.51 0.21
x2  .  1.00 0.31 0.74 0.58 0.58 0.55 0.19 0.07 0.16 0.21 0.41 0.47 0.17 0.35 0.16
x3  .   .   1.00 0.38 0.31 0.30 0.35 0.58 0.28 0.33 0.38 0.35 0.41 0.64 0.58 0.51
x4  .   .   .    1.00 0.90 0.78 0.75 0.25 0.20 0.22 0.08 0.53 0.79 0.27 0.57 0.26
x5  .   .   .    .    1.0 0.79 0.74 0.25 0.18 0.23 -.02 0.48 0.79 .27 .51 .23
x6  .   .   .    . .  1 .73 .18 .18 .23 .00 .38 .69 .14 .26 .00
x7  .   .   .    . . . 1 .24 .29 .25 .10 .44 .67 .16 .38 .12
x8  .   .   .    . . . . 1 -.04 .49 .44 .30 .32 .51 .51 .38
x9  .   .   .    . . . . . 1 -.34 -.16 -.05 .23 .21 .15 .18
x10 .   .   .    . . . . . . 1 .23 .50 .31 .15 .29 .14
x11 .   .   .    . . . . . . . 1 .24 .10 .31 .28 .31
x12 .   .   .    . . . . . . . . 1 .62 .17 .41 .18
x13 .   .   .    . . . . . . . . . 1 .26 .50 .24
x14 .   .   .    . . . . . . . . . . 1 .63 .50
x15 .   .   .    . . . . . . . . . . . 1 .65
x16 .   .   .    . . . . . . . . . . . . 1
;
proc print;
run;
proc princomp data=d731(type=corr) n=3 outstat=s731 ;
   var x1-x16;
run;

proc print data=s731;
run;

data t731;
  set s731;
  if _type_='SCORE';
run;
proc transpose data=t731 out=tt731 prefix=fac;
   var x1-x16;
run;
data t731;
  set tt731;
  p1=sqrt(7.03648)*fac1;
  p2=sqrt(2.61403)*fac2;
run;

options ps=35 ls=80;

proc plot data=t731;
   plot p2*p1 $ _name_='*'/vref=0 ;
run;
proc gplot data=t731;
   plot p2*p1 ;
run;