求高人指导SAS总主成分分析综合得分程序

我按照你说的做了，去掉后面3步程序虽然是按降序排列，仔细看一下你截图的结果，0.7y应该在第三位，结果却排在负数后面，其他不走都是正确的，我需要的结果是按y的值把这14个城市按降序排列，麻烦你再帮我看看，我就是这个排序不会做。

明月独醉清风 发表于 2013-4-7 08:07
我按照你说的做了，去掉后面3步程序虽然是按降序排列，仔细看一下你截图的结果，0.7y应该在第三位，结果却排 ...

这个要冒泡排序，还得把数据跟着排序，很麻烦。

那就不用了，谢谢你！

proc sort data=gszjia;把这步中的sort改成order后就可以按降序排列了，但问题是排序中有两个数据不合适，再麻烦您看看。

proc sort data=gszjia;把这步中的sort改成order后就可以按降序排列了，但问题是排序中有两个数据不合适，再麻烦您看看。

明月独醉清风 发表于 2013-4-7 13:39
proc sort data=gszjia;把这步中的sort改成order后就可以按降序排列了，但问题是排序中有两个数据不合适，再 ...

我和你一样，。没注意细节，sort过程有个equals选项，可以按数值大小排序。

1. proc sort data=gszjia equals;
   
2.     by descending y ;
   
3. run;

复制代码

这样就对了。

1. data gs2;
   
2.     input name $ x1-x11;
   
3. cards;
   
4. 
   
5. 兰州市 1.43 1.02 1.36 1.04 2.92 5.90 5.50 1.11 1.01 4.90 0.70
   
6. 嘉峪关 3.84 1.73 0.47 2.40 0.22 0.55 0.26 1.32 0.92 0.29 0.02
   
7. 金昌市 1.90 1.65 0.46 2.09 0.42 0.47 0.37 1.40 0.96 0.38 0.11
   
8. 白银市 0.83 1.21 0.88 1.33 0.74 0.76 0.88 1.11 0.92 0.73 1.06
   
9. 天水市 0.42 0.84 1.16 0.63 0.98 1.01 1.29 0.91 1.00 1.21 1.81
   
10. 武威市 0.57 0.89 0.93 0.67 0.96 0.34 0.76 0.93 1.09 0.88 0.88
    
11. 张掖市 0.81 0.79 0.97 0.68 0.49 0.42 0.69 0.87 0.85 0.67 0.48
    
12. 平凉市 0.51 1.02 0.87 0.98 0.97 0.61 0.89 0.93 1.07 0.80 0.93
    
13. 酒泉市 1.66 1.10 0.99 1.00 1.81 1.01 0.90 1.20 0.88 1.05 0.58
    
14. 庆阳市 0.78 1.34 0.67 1.57 1.73 1.78 0.96 1.00 0.93 0.92 1.39
    
15. 定西市 0.26 0.58 1.23 0.26 0.91 0.30 0.54 0.86 0.93 0.70 1.68
    
16. 陇南市 0.29 0.65 1.23 0.48 0.97 0.46 0.42 0.85 1.20 0.73 1.81
    
17. 临夏州 0.25 0.65 1.35 0.37 0.47 0.23 0.35 0.68 1.24 0.50 1.89
    
18. 甘南州 0.45 0.53 1.44 0.38 0.41 0.16 0.20 0.82 0.99 0.24 0.68
    
19. ;
    
20. 
    
21. proc princomp data=gs2 prefix=z n=11 out=result(drop=name x1-x11);
    
22.     var x1-x11;
    
23. run;
    
24. 
    
25. proc print data=result;
    
26.     var z1-z11;
    
27. run;
    
28. 
    
29. data gszjia;
    
30.     set gs2(keep=name);
    
31.     set result;
    
32.     y=0.4765*z1+0.3528*z2+0.0741*z3+0.0478*z4+0.0190*z5+0.0149*z6+0.0098*z7+0.0044*z8+0.0007*z9+0.0001*z10;
    
33. run;
    
34. 
    
35. proc sort data=gszjia equals;
    
36.     by descending y ;
    
37. run;
    
38. 
    
39. proc print data=gszjia;
    
40.     var name y z1-z11;
    
41. run;
    
42. 
    
43. proc datasets lib=work;
    
44.     delete result;
    
45. quit;

复制代码

嗯，看了一下，这回合适了，非常感谢你，

请问，主成分分析不是应该基于协方差或相关系数矩阵做吗？data=gs2 这个数据集是原始数据还是相关系数矩阵？

Kalet 发表于 2013-4-7 00:24
应该不是特征根。

特征值对应的特征向量