多元统计分析主成分分析

1、主成分分析的基本思想主成分分析的基本思想主成分的计算主成分的计算主成分的性质主成分的性质主成分分析的应用主成分分析的应用主成分回归主成分回归 一项十分著名的工作是美国的统计学家斯通(stone)在1947年关于国民经济的研究。他曾利用美国1929一1938年各年的数据，得到了17个反映国民收入与支出的变量要素，例如雇主补贴、消费资料和生产资料、纯公共支出、净增库存、股息、利息外贸平衡等等。1 基本思想 在进行主成分分析后，竟以97.4的精度，用三新变量就取代了原17个变量。根据经济学知识，斯通给这三个新变量分别命名为总收入总收入f1、总收入总收入变化率变化率f2和经济发展趋势经济发展趋势f3。

2、更有意思的是，这三个变量其实都是可以直接测量的。斯通将他得到的主成分与实际测量的总收入i、总收入变化率i以及时间t因素做相关分析，得到下表：f1f1f2f2f3f3i ii it tf1f11 1f2f20 01 1f3f30 00 01 1i i0.9950.995-0.041-0.0410.0570.057l li i-0.056-0.0560.9480.948-0.124-0.124-0.102-0.102l lt t-0.369-0.369-0.282-0.282-0.836-0.836-0.414-0.414-0.112-0.1121 1主成分分析的基本思想 主成分分析就是把原有的多

3、个指标转化成少数几个代表性较好的综合指标，这少数几个指标能够反映原来指标大部分的信息（85%以上），并且各个指标之间保持独立，避免出现重叠信息。主成分分析主要起着降维和简化数据结构的作用。 主成分分析是把各变量之间互相关联的复杂关系进行简化分析的方法。 在社会经济的研究中，为了全面系统的分析和研究问题，必须考虑许多经济指标，这些指标能从不同的侧面反映我们所研究的对象的特征，但在某种程度上存在信息的重叠，具有一定的相关性。 主成分分析试图在力保数据信息丢失最少的原则下，对这种多变量的截面数据表进行最佳综合简化，也就是说，对高维变量空间进行降维处理。 很显然，识辨系统在一个低维空间要比在一个高维空

4、间容易得多。2 数学模型与几何解释 假设我们所讨论的实际问题中，有p个指标，我们把这p个指标看作p个随机变量，记为x1，x2，xp，主成分分析就是要把这p个指标的问题，转变为讨论 m 个新的指标f1，f2，fm(mp），按照保留主要信息量的原则充分反映原指标的信息，并且相互独立。npnnppxxxxxxxxxx212222111211niiiixxxx21其中pxxx21xxaxaxafxxaxaxafxxaxaxafppppppppppp2211222221122112211111 这种由讨论多个指标降为少数几个综合指标的过程在数学上就叫做降维。主成分分析通常的做法是，寻求原指标的线性组合f

5、i。)()(121xvarkxkvar所以如果不对 加以限制，问题就变得无意义。1最大最大因此限制 为单位向量。1)()(1piiixavarxvarijjpjijiipiiiisaasa1,122piijjpjisaa11)(xvar满足如下的条件：122221piiiaaapjijiffcovji，），（210)()(21pfvarfvarfvar）（主成分之间相互独立，即无重叠的信息。即主成分之间相互独立，即无重叠的信息。即主成分的方差依次递减，重要性依次递减，即主成分的方差依次递减，重要性依次递减，即每个主成分的系数平方和为每个主成分的系数平方和为1。即。即2x1x1f2f主成分分析的

6、几何解释平移、旋转坐标轴2x1x1f2f主成分分析的几何解释平移、旋转坐标轴2x1x1f2f 主成分分析的几何解释平移、旋转坐标轴 旋转变换的目的是为了使得n个样品点在fl轴方向上的离散程度最大，即fl的方差最大。变量fl代表了原始数据的绝大 部分信息，在研究某经济问题时，即使不考虑变量f2也无损大局。经过上述旋转变换原始数据的大部分信息集中到fl轴上，对数据中包含的信息起到了浓缩作用。 fl，f2除了可以对包含在xl，x2中的信息起着浓缩作用之外，还具有不相关的性质，这就使得在研究复杂的问题时避免了信息重叠所带来的虚假性。二维平面上的n个点的方差大部分都归结在fl轴上，而f2轴上的方差很小。

7、fl和f2称为原始变量x1和x2的综合变量。f简化了系统结构，抓住了主要矛盾。 2x1x1f2f主成分分析的几何解释平移、旋转坐标轴主成分的计算先讨论二维情形212122211211xxxxxxxxxnn求主成分f1和f2。21,xx观察图，我们已经把主成分f1和f2 的坐标原点放在平均值 所在处，从而使得f1和f2 成为中心化的变量，即f1和f2 的样本均值都为零。因此f1可以表示为)()(222111111xxaxxaf),(2111aa关键关键是，寻找合适的单位向量 ，使f1的方差最大。122111222211121112)(saasasafvar2111222112112111)(aa

8、ssssaa最大最大1问题的答案问题的答案是：x的协方差矩阵s 的最大特征根 所对应的单位特征向量即为 。并且 就是f1的方差。2111,aa1推导同样，f2可以表示为)()(222211122xxaxxaf),(2212aa寻找合适的单位向量 ，使f2与f1独立，且使f2的方差（除f1之外）最大。2问题的答案问题的答案是：x的协方差矩阵s 的第二大特征根 所对应的单位特征向量即为 。并且 就是f2的方差。2212,aa2推导求解主成分的步骤：求解主成分的步骤：1. 求样本均值 和样本协方差矩阵s;),(21xxx 2. 求s的特征根求解特征方程 ，其中i是单位矩阵，解得2个特征根 0 is2

9、121,3. 求特征根所对应的单位特征向量4. 写出主成分的表达式例1 下面是8 个学生两门课程的成绩表 6585709065455565数学10090707085555545语文1x2x对此进行主成分分析。1. 求样本均值和样本协方差矩阵求样本均值和样本协方差矩阵5 .6725.7121xxx5 .1871 .1034 .323sis2. 求解特征方程 0 05 .1871 .1031 .1034 .32301 .103)5 .187)(4 .323(2 化简得：09 .500079 .5102 解得：132, 9 .37821 5 .1871 .1034 .323s3.求特征值所对应的单位

10、特征向量 1所对应的单位特征向量 ， 0)(11s 其中21111aa0)9 .3785 .187(1 .10301 .103)9 .3784 .323(21112111aaaa1221211aa解得 （2111,aa）=)47. 0 ,88. 0( 2所对应的单位特征向量 0)(22s ，其中221220)1325 .187(1 .10301 .103)1324 .323(22122212aaaa1222212aa解得： )88. 0 ,47. 0(),(2212aa5 .1871 .1034 .323s4. 得到主成分的表达式 )5 .67(47. 0)25.71(88. 0211xxf第

11、二主成分：)5 .67(88. 0)25.71(47. 0212xxf第一主成分：5.主成分的含义通过分析主成分的表达式中原变量前的系数来解释各主成分的含义。 第一主成分f1是 和 的加权和，表示该生成绩的好坏。 1x2x第二主成分f2表示学生两科成绩的均衡性 6. 比较主成分重要性比较主成分重要性 第一主成分f1的方差为9 .3781第二主成分f2的方差为1322方差贡献率 )()()(211211fvarfvarfvar%16.741329 .3789 .378%84.251329 .378132212方差贡献率为 主成分f1和f2的方差总和为219 .5101329 .378原变量和1x

12、2x的方差总和为9 .5105 .1874 .3232211ss总方差保持不变总方差保持不变身高x1(cm)胸围x2(cm)体重x3(kg)149.5162.5162.7162.2156.5156.1172.0173.2159.5157.769.577.078.587.574.574.576.581.574.579.038.555.550.865.549.045.551.059.543.553.5例2 下表是10位学生的身高1x、胸围2x、体重3x的数据。对此进行主成分分析。 1. 求样本均值和样本协方差矩阵 2 .513 .772 .161321xxx53.5558.3200.3011.21

13、12.1767.46s 2. 求解协方差矩阵的特征方程 0 is 053.5558.3200.3058.3211.2112.1700.3012.1767.463.解得三个特征值 15.98160.23256. 13)71. 0 ,42. 0 ,56. 0(),(312111aaa)48. 0,33. 0,81. 0(),(322212aaa)53. 0,85. 0 ,03. 0(),(332313aaa和对应的单位特征向量：4. 由此我们可以写出三个主成分的表达式： )2 .51(71. 0)3 .77(42. 0)2 .161(56. 03211xxxf)2 .51(48. 0)3 .77(

14、33. 0)2 .161(81. 03212xxxf)2 .51(53. 0)3 .77(85. 0)2 .161(03. 03213xxxf5. 主成分的含义f1表示学生身材大小。 f2反映学生的体形特征 三个主成分的方差贡献率分别为：%6 .7931.12315.9856. 160.2315.9815.98311ii%1 .1931.12360.23312ii%3 . 131.12356. 1313ii前两个主成分的累积方差贡献率为： %7 .9831.12375.1213121ii5 . 155. 030. 008. 021. 075. 0543212xxxxxf8 .1960. 060

15、. 015. 042. 030. 0543213xxxxxf1 .1118. 052. 000. 078. 030. 0543214xxxxxf9 .1315. 029. 092. 019. 008. 0543215xxxxxf例3 对88个学生5 门不同课程的考试成绩进行分析，要求用合适的方法对这5 门课程成绩进行平均，以对88个学生的成绩进行评比。这5门课程是：mechanics vectors (闭)，algebra analysis statistics (开)。1x2x4x3x5x经计算，得到5个主成分的表达式如下： 这5个主成分的方差分别为679.2，199.8，102.6, 83

16、.7和31.8。前两个主成分各自的贡献率和累积贡献率为%91.611 .10972 .679511ii%21.181 .10978 .199512ii%12.80%21.18%91.615121ii在一般情况下，设有n个样品，每个样品观测p个指 标，将原始数据排成如下矩阵： npnnppxxxxxxxxx.212222111211),.,(21pxxxx1.求样本均值和样本协方差矩阵s;2.求解特征方程is=0, 其中i是单位矩阵0.212222111211ppppppsssssssss,解得p个特征根p,.,21).(21p3. 求k所对应的单位特征向量k),.,2 , 1(pk 即需求解方

17、程组0)(kkis其中),.,(21pkkkkaaa 0.21212222111211pkkkkpppppkpkaaasssssssss 再加上单位向量的条件 1.22221pkkkaaa解得 ),.,(21pkkkkaaa4. 写出主成分的表达式 )(.)()(222111pppkkkkxxaxxaxxaf根据累积贡献率的大小取前面m 个(m|t|interceptx1x2x31111-10.12799-0.051400.586950.286851.212160.070280.094620.10221-8.36-0.736.202.810.00010.48830.00040.0263para

18、meterestimatesdependentmean21.89091r-square0.9919rootmse0.48887adjr-sq0.9884summaryoffitf1f2f3x1x2x30.70630.04350.7065-0.03570.9990-0.02580.70700.0070-0.7072eigenvectorseigenvaluedifferenceproportioncumulativepcr1pcr2pcr31.99920.99820.00261.00100.99550.66640.33270.00090.66640.99911.0000eigenvaluesof

19、thecorrelationmatrixf1=0.7063x1+0.0435x2+0.7065x3f2=-0.0357x1+0.9990 x2-0.0258x3 sourcedfsumofsquaresmeansquaref值值probfmodelerrortotal28109.88280.117210.00004.94140.0147337.23020.0001analysisofvariancevariabledfestimatestandarderrort值值prob|t|f1f2110.69000.19130.02710.038325.48594.99300.00010.0011par

20、ameterestimates2119130. 068998. 0ffy*3*2*14825. 02211. 04804. 0 xxxy3211062. 06091. 00727. 0130. 9xxxy标准化后的变量63.2073.1394826. 065. 13 . 322. 03059.1944805. 05437. 489.21321xxxy把标准化变量还原，代入得：影响人们外出旅游的因素有居民收入、交通、闲影响人们外出旅游的因素有居民收入、交通、闲暇时间、旅游目的地治安状况、旅游目的地的环暇时间、旅游目的地治安状况、旅游目的地的环境卫生以及接待能力等等。境卫生以及接待能力等等。 由于

21、资料的可得性和代表性，选择以下变量由于资料的可得性和代表性，选择以下变量。 国内旅游人数（百万人）y农村居民人均纯收入（元）1x城镇居民人均可支配收入（元）2x公路线路里程（万公里） 3x数据见sasuser.tourmx例例2 国内旅游人数模型国内旅游人数模型variabledfestimatestandarderrort值值prob|t|interceptincomeonincomeochighway1111417.8201-0.13810.1737-3.000974.02300.06990.03020.81925.6445-1.97595.7589-3.66330.00050.08360

22、.00040.0064parameterestimatesdependentmean558.1017r-square0.9920rootmse19.2003adjr-sq0.9890summaryoffitf1f2f3x1x2x30.58100.59180.5588-0.5167-0.26230.81500.6289-0.76220.1533eigenvectorseigenvaluedifferenceproportioncumulativepcr1pcr2pcr32.80880.18500.00622.62380.17880.93630.06170.00210.93630.99791.00

23、00eigenvaluesofthecorrelationmatrixf1=0.5810 x1+0.5918x2+0.5588x3f2=-0.5167x1-0.2623x2+0.8150 x3 sourcedfsumofsquaresmeansquaref值值probfmodelerrortotal291110.71130.288711.00005.35560.0321166.93280.0001analysisofvariancevariabledfestimatestandarderrort值值prob|t|f1f2110.5767-0.46200.03220.125617.8977-3.

24、67940.00010.0051parameterestimates214620. 05767. 0ffy*3*2*11741. 03231. 03213. 0 xxxy标准化后的变量36.1990.1211741. 084.186566.41673231. 026.67063.15753213. 091.18210.558321xxxy把标准化变量还原，代入得：3216448. 103167. 008768. 08482.286xxxy 地区经济发展现状及潜力分析 长江三角洲经济发展状况分析 长江三角洲产业发展状况分析 城市竞争力评价指标体系 区域智力资本的测度 区域创新能力对经济增长的影响分析 区域智力资本对经济增长的影响分析 区域软实力评价体系研究主成分的推导主成分的推导 （一）（一）第一主成分第一主成分xxaxafpp1111111111)()( xvarfvarmax寻找合适的单位向量 ，使f1的方差最大。1) 1(1111 q022111q022111q11表明： 应为 的特征值，而 为与 对应的单位特征向量。 1 11111)(fvar而且可见 应取 的最大特征根。如果第一主成分的信息不够，则需要寻找第二主成分。xxaxafpp221122（