主成分的求解方法
1求相关矩阵
2、求特征值与特征向量
3、确定主成分个数
4、计算主成分得分。
data p108;
input x$ x1-x8;
datalines;
北京1394.89 2505.00 519.01 8144 373.90 117.30 112.60 843.43 天津920.11 2720.00 345.46 6501 342.80 115.20 110.60 582.51 河北2849.52 1258.00 704.87 4839 2033.30 115.20 115.80 1234.85 山西1092.48 1250.00 290.90 4721 717.30 116.90 115.60 697.25 内蒙832.88 1387.00 250.23 4134 781.70 117.50 116.80 419.39 辽宁2793.37 2397.00 387.99 4911 1371.10 116.10 114.00 1840.55 吉林1129.20 1872.00 320.45 4430 497.40 115.20 114.20 762.47 黑龙江2014.53 2334.00 435.73 4145 824.80 116.10 114.30 1240.37 上海2462.57 5343.00 996.48 9279 207.40 118.70 113.00 1642.95 江苏5155.25 1926.00 1434.95 5943 1025.50 115.80 114.30 2026.64 浙江3524.79 2249.00 1006.39 6619 754.40 116.60 113.50 916.59 安徽2003.58 1254.00 474.00 4609 908.30 114.80 112.70 824.14 福建2160.52 2320.00 553.97 5857 609.30 115.20 114.40 433.67 江西1205.11 1182.00 282.84 4211 411.70 116.90 115.90 571.84 山东5002.34 1527.00 1229.55 5145 1196.60 117.60 114.20 2207.69 河南3002.74 1034.00 670.35 4344 1574.40 116.50 114.90 1367.92 湖北2391.42 1527.00 571.68 4685 849.00 120.00 116.60 1220.72 湖南2195.70 1408.00 422.61 4797 1011.80 119.00 115.50 843.83 广东5381.72 2699.00 1639.83 8250 656.50 114.00 111.60 1396.35 广西1606.15 1314.00 382.59 5105 556.00 118.40 116.40 554.97 海南364.17 1814.00 198.35 5340 232.10 113.50 111.30 64.33 四川3534.00 1261.00 822.54 4645 902.30 118.50 117.00 1431.81 贵州630.07 942.00 150.84 4475 301.10 121.40 117.20 324.72 云南1206.68 1261.00 334.00 5149 310.40 121.30 118.10 716.65 西藏55.98 1110.00 17.87 7382 4.20 117.30 114.90 5.57 陕西1000.03 1208.00 300.27 4396 500.90 119.00 117.00 600.98 甘肃553.35 1007.00 114.81 5493 507.00 119.80 116.50 468.79 青海165.31 1445.00 47.76 5753 61.60 118.00 116.30 105.80
宁夏169.75 1355.00 61.98 5079 121.80 117.10 115.30 114.40 新疆834.57 1469.00 376.95 5348 339.00 119.70 116.70 428.76 ;
proc princomp;
var x1-x8;
run;
proc princomp out=score n=3;
var x1-x8;
run;
proc print data=score;
run;
Proc plot data=score;
Plot prin2*prin1=x;
Run;
例:进行主成分分析
1、用程序求样本相关矩阵的特征值和特征向量。
2、求两个主成分的累计贡献率,写出贡献率计算表达式。
3、写出两个主成分的模型表达式。
4、计算主成分得分,并按主成分得分名次
例:对下表的指标进行主成分分析,
5、用程序求样本相关矩阵的特征值和特征向量。
6、求两个主成分的累计贡献率,写出贡献率计算表达式。
7、写出两个主成分的模型表达式。
8、计算主成分得分,并按主成分得分名次
9、按两个主成分得分作散点图
10、对主成分得分聚类
例:
学生身体各指标的主成分分析.
随机抽取30名某年级中学生,测量其身高(X1)、体重(X2)、胸围(X3)和坐高(X4),数据见
试对中学生身体指标数据做主成分分析.。