第七章主成分分析.pptx

《第七章主成分分析.pptx》由会员分享，可在线阅读，更多相关《第七章主成分分析.pptx（54页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、寻找主成分的正交旋转 旋转公式：第1页/共54页7.2 总体的主成分一、主成分的定义及导出二、主成分的性质三、从相关阵出发求主成分第2页/共54页一、主成分的定义及导出设 为一个p维随机向量，E(x)=，V(x)=。考虑如下的线性变换希望在约束条件 下寻求向量a1，使得 达到最大，y1就称为第一主成分。设1 2p0为的特征值，i=1,2,p为相应的单位特征向量，且相互正交。则可求得第一主成分为它的方差具有最大值1。第3页/共54页如果第一主成分所含信息不够多，还不足以代表原始的p个变量，则需考虑再使用一个综合变量 ，为使y2所含的信息与y1不重叠，应要求Cov(y1,y2)=0我们在此条件和约

2、束条件 下寻求向量a2，使得 达到最大，所求的 称为第二主成分。求得的第二主成分为其方差为2。一般来说，x的第i主成分是指：在约束条件 和 Cov(yk,yi)=0,k=1,2,i1下寻求ai，使得达到最大。第i主成分为第4页/共54页主成分的几何意义在几何上，ti表明了第i主成分的方向，yi是x在ti上的投影值（其绝对值即为投影长度），i是这些值的方差，它反映了在ti上投影点的分散程度。第5页/共54页x投影到ti上的值其中i是ti与x的夹角。第6页/共54页主成分向量与原始向量之间的关系式第7页/共54页第8页/共54页主成分与原始变量之间的关系式矩阵y1y2ypx1t11t12t1px2

3、t21t22t2pxptp1tp 2tpp第9页/共54页正交变换 的几何意义正交变换 的几何意义是将Rp中由x1,x2,xp构成的原p维坐标轴作一正交旋转，一组正交单位向量t1,t2,tp表明了p个新坐标轴的方向，这些新坐标轴彼此仍保持正交（或说垂直）。第10页/共54页二、主成分的性质1.主成分向量的协方差矩阵2.主成分的总方差 3.原始变量xi与主成分yk之间的相关系数4.m个主成分对原始变量的贡献率5.原始变量对主成分的影响第11页/共54页1.主成分向量的协方差矩阵V(y)=其中=diag(1,2,p)，即V(yi)=i,i=1,2,p，且y1,y2,yp互不相关。第12页/共54页

4、2.主成分的总方差 由于所以 或第13页/共54页总方差中属于第i主成分yi（或被yi所解释）的比例为称为主成分yi的贡献率。第一主成分y1的贡献率最大，表明它解释原始变量 x1,x2,xp的能力最强，而y2,y3,yp的解释能力依次递减。主成分分析的目的就是为了减少变量的个数，因而一般是不会使用所有p个主成分的，忽略一些带有较小方差的主成分将不会给总方差带来大的影响。第14页/共54页前m个主成分的贡献率之和称为主成分y1,y2,ym的累计贡献率，它表明y1,y2,ym解释x1,x2,xp的能力。通常取（相对于p）较小的m，使得累计贡献达到一个较高的百分比（如8090）。此时，y1,y2,y

5、m可用来代替x1,x2,xp，从而达到降维的目的，而信息的损失却不多。第15页/共54页3.原始变量xi与主成分yk之间的相关系数 x=Ty即xi=ti1y1+ti2y2+tipyp所以Cov(xi,yk)=Cov(tikyk，yk)=tikk在实际应用中，通常我们只对xi(i=1,2,p)与yk(k=1,2,m)的相关系数感兴趣。第16页/共54页4.m个主成分对原始变量的贡献率m个主成分y1,y2,ym从原始变量x1,x2,xp中提取的信息量，可度量为xi与y1,y2,ym的复相关系数的平方，称为m个主成分y1,y2,ym对原始变量xi的贡献率。其值为当m=p时，第17页/共54页例7.2

6、.1 设x=(x1,x2,x3)的协方差矩阵为其特征值为1=5.83，2=2.00，3=0.17相应的特征向量为若只取一个主成分，则贡献率为5.83/(5.83+2.00+0.17)=0.72875=72.875%第18页/共54页可见，y1对第三个变量的贡献率为零，这是因为x3与x1和x2都不相关，在y1中未包含一点有关x3的信息，这时仅取一个主成分就显得不够了，故应再取y2，此时累计贡献率为(5.83+2.00)/8=97.875%(y1，y2)对每个变量xi的贡献率分别为 ，都比较高。表7.2.1 y1及(y1,y2)对每个原始变量的贡献率i(y1,xi)(y2,xi)10.9250.8

7、550.0000.85520.9980.9960.0000.99630.0000.0001.0001.000第19页/共54页5.原始变量对主成分的影响 yk=t1kx1+t2kx2+tpkxp称tik为第k主成分yk在第i个原始变量xi上的载荷，它度量了xi对yk的重要程度。在解释主成分时，我们需要考察载荷，同时也应考察一下相关系数。方差大的那些变量与具有大特征值的主成分有较密切的联系，而方差小的另一些变量与具有小特征值的主成分有较强的联系。通常我们取前几个主成分，因此所取主成分会过于照顾方差大的变量，而对方差小的变量却照顾得不够。第20页/共54页例7.2.2 设x=(x1,x2,x3)的

8、协方差矩阵为经计算，的特征值及特征向量为1=109.793，2=6.469，3=0.738 相应的主成分分别为第21页/共54页y1=0.305x1+0.041x2+0.951x3y2=0.944x1+0.120 x20.308x3y3=0.127x1+0.992x20.002x3可见，方差大的原始变量x3在很大程度上控制了第一主成分y1，方差小的原始变量x2几乎完全控制了第三主成分y3，方差介于中间的x1则基本控制了第二主成分y2。y1的贡献率为这么高的贡献率首先归因于x3的方差比x1和x2的方差大得多，其次是x1,x2,x3相互之间存在着一定的相关性。y3的特征值相对很小，表明x1,x2,

9、x3之间有这样一个线性依赖关系：0.127x1+0.992x20.002x3c其中c=0.1271+0.99220.0023为一常数。第22页/共54页三、从相关阵出发求主成分最常用的标准化变换是令 。显然，的协方差矩阵正是x的相关矩阵R。从R出发求主成分，主成分分析将均等地对待每一个原始变量。从R出发求得主成分的方法与从出发是完全类似的，并且主成分的一些性质具有更简洁的数学形式。设 为R的p个特征值，为相应的单位特征向量，且相互正交，则p个主成分为 。记 ，于是y*=T*x*第23页/共54页从R出发的主成分性质(1)E(y*)=0，V(y*)=*，其中 (2)。(3)变量 与主成分 之间的

10、相关系数即有第24页/共54页 因此，在解释主成分 时，由相关矩阵R求得的载荷 和相关系数 所起的作用是完全相同的，只需选其一用来作主成分解释即可。(4)主成分 对变量 的贡献率(5)。第25页/共54页例7.2.3 在例7.2.2中，x的相关矩阵R的特征值及特征向量为相应的主成分分别为第26页/共54页 的贡献率为 和 累计贡献率为现比较本例中从R出发和例7.2.2中从 出发的主成分计算结果。从R出发的 的贡献率0.705明显小于从出发的y1的贡献率0.938，事实上，原始变量方差之间的差异越大，这一点也就倾向于越明显。可用标准化前的原变量表达如下：第27页/共54页可见，在原变量x1,x2

11、,x3上的载荷相对大小与例7.2.2中yi在x1,x2,x3上的载荷相对大小之间有着非常大的差异。这说明，标准化后的结论完全可能会发生很大的变化，因此标准化不是无关紧要的。第28页/共54页7.3 样本的主成分我们可以从协差阵或相关阵R出发求得主成分。但在实际问题中，或R一般都是未知的，需要通过样本来进行估计。设数据矩阵为则样本协差阵和样本相关阵分别为第29页/共54页7.3 样本的主成分一、样本主成分的定义二、从S出发求主成分三、从 出发求主成分四、主成分分析的应用五、若干补充及应用中需注意的问题第30页/共54页一、样本主成分的定义若向量a1在约束条件 下，使得的样本方差 达到最大，则称线

12、性组合 为第一样本主成分。若向量a2在约束条件 和的样本协方差第31页/共54页下，使得 的样本方差 达到最大，则称线性组合 为第二样本主成分。一般地，若向量ai 在约束条件 和 的样本协方差第32页/共54页下，使得的样本方差达到最大，则称线性组合 为第i样本主成分，i=1,2,p。需要指出的是，样本主成分是使样本方差而非方差达到最大，是使样本协方差而非协方差为零。第33页/共54页二、从S出发求主成分用类似于上一节的方法，以S代替即可求得样本主成分。设 为S的特征值，为相应的单位特征向量，且彼此正交。则第i样本主成分为 ，它具有样本方差,i=1,2,p，各主成分之间的样本协方差为零。在几何

13、上，p个样本主成分的方向为 所在的方向，且彼此垂直。n个样品点在 上的投影点最为分散，在其余 上投影点的分散程度依次递减。第34页/共54页总样本方差xi与 的样本相关系数其中 ，k=1,2,p。第35页/共54页主成分得分在实际应用中，我们常常让xj 减去 ，使样本数据中心化。这不影响样本协差阵S，在前面的论述中惟一需要变化的是，将第i主成分改写成中心化的形式，即若将各观测值xj代替上式中的观测值向量x，则第i主成分的值称之为观测值xj的第i主成分得分。所有观测值的平均主成分得分第36页/共54页三、从 出发求主成分设样本相关阵 的p个特征值为 ，为相应的正交单位特征向量，则第i样本主成分其

14、中x*是各分量经（样本）标准化了的向量，即第37页/共54页标准化后的主成分得分令这是xj的各分量数据经标准化后的数据向量，将其代替上述样本主成分公式中的x*，即得观测值xj在第i主成分上的得分所有观测值的平均主成分得分第38页/共54页四、主成分分析的应用在主成分分析中，我们首先应保证所提取的前几个主成分的累计贡献率达到一个较高的水平，其次对这些被提取的主成分必须都能够给出符合实际背景和意义的解释。主成分的解释其含义一般多少带有点模糊性，不像原始变量的含义那么清楚、确切，这是变量降维过程中不得不付出的代价。因此，提取的主成分个数m通常应明显小于原始变量个数p（除非p本身较小），否则维数降低的

15、“利”可能抵不过主成分含义不如原始变量清楚的“弊”。第39页/共54页如果原始变量之间具有较高的相关性，则前面少数几个主成分的累计贡献率通常就能达到一个较高水平，也就是说，此时的累计贡献率通常较易得到满足。主成分分析的困难之处主要在于要能够给出主成分的较好解释，所提取的主成分中如有一个主成分解释不了，整个主成分分析也就失败了。主成分分析是变量降维的一种重要、常用的方法，简单的说，该方法要应用得成功，一是靠原始变量的合理选取，二是靠“运气”。第40页/共54页例7.3.1 在制定服装标准的过程中，对128名成年男子的身材进行了测量，每人测得的指标中含有这样六项：身高（x1）、坐高（x2）、胸围（

16、x3）、手臂长（x4）、肋围（x5）和腰围（x6）。所得样本相关矩阵列于表7.3.1。x1x2x3x4x5x6x11.000 x20.791.000 x30.360.311.000 x40.760.550.351.000 x50.250.170.640.161.000 x60.510.350.580.380.631.000表7.3.1 男子身材六项指标的样本相关矩阵第41页/共54页经计算，相关阵 的前三个特征值、相应的特征向量以及贡献率列于表7.3.2。表7.3.2 的前三个特征值、特征向量以及贡献率特征向量 ：身高0.4690.3650.092 ：坐高0.4040.3970.613 ：胸围

17、0.3940.3970.279 ：手臂长0.4080.3650.705 ：肋围0.3370.5690.164 ：腰围0.4270.3080.119特征值3.2871.4060.459贡献率0.5480.2340.077累计贡献率0.5480.7820.859第42页/共54页前三个主成分分别为从表7.3.2中可以看到，前两个主成分的累计贡献率已达78.2，前三个主成分的累计贡献率达85.9，因此可以考虑只取前面两个或三个主成分，它们能够很好地概括原始变量。第一主成分 对所有（标准化）原始变量都有近似相等的正载荷，故称第一主成分为（身材）大小成分。第43页/共54页第二主成分 在 上有中等程度的

18、正载荷，而在 上有中等程度的负载荷，称第二主成分为形状成分（或胖瘦成分）。第三主成分 在 上有大的正载荷，在 上有大的负载荷，而在其余变量上的载荷都较小，可称第三主成分为臂长成分。由于第三主成分的贡献率不高（7.65）且实际意义也不太重要，因此我们一般可考虑取前两个主成分。由于 非常小，所以存在共线性关系：第44页/共54页例7.3.2 在习题6.5中，如下八项男子径赛运动记录：x1：100米（秒）x5：1500米（分）x2：200米（秒）x6：5000米（分）x3：400米（秒）x7：10000米（分）x4：800米（秒）x8：马拉松（分）表7.3.3 八项男子径赛运动记录的样本相关矩阵x1

19、x2x3x4x5x6x7x8x11.000 x20.9231.000 x30.8410.8511.000 x40.7560.8070.8701.000 x50.7000.7750.8350.9181.000 x60.6190.6950.7790.8640.9281.000 x70.6330.6970.7870.8690.9350.9751.000 x80.5200.5960.7050.8060.8660.9320.9431.000第45页/共54页表7.3.4 的前三个特征值、特征向量以及贡献率特征向量 ：100米0.3180.5670.332 ：200米0.3370.4620.361 ：40

20、0米0.3560.2480.560 ：800米0.3690.0120.532 ：1500米0.3730.1400.153 ：5000米0.3640.3120.190 ：10000米0.3670.3070.182 ：马拉松0.3420.4390.263特征值6.6220.8780.159贡献率0.8280.1100.020累计贡献率0.8280.9370.957第46页/共54页例7.3.3 对例6.3.3中的数据从相关矩阵出发进行主成分分析。经计算，x1,x2,x8的样本相关矩阵 列于表7.3.5。的前三个特征值、特征向量以及贡献率列于表7.3.6。表7.3.5 消费性支出八个变量的样本相关矩

21、阵x1x2x3x4x5x6x7x8x11.000 x20.2471.000 x30.6980.2581.000 x40.4680.4230.6211.000 x50.8280.0860.5850.5311.000 x60.7690.2550.8560.6840.7081.000 x70.6700.2010.5690.3140.8000.6471.000 x80.8770.3490.6670.6280.7760.7450.5251.000第47页/共54页表7.3.6 的前三个特征值、特征向量以及贡献率特征向量 ：食品0.4010.0770.415 ：衣着0.1320.7490.332 ：家庭设

22、备用品及服务0.3750.0650.442 ：医疗保健0.3200.3450.478 ：交通和通讯0.3880.2320.279 ：娱乐教育文化服务0.4060.0270.310 ：居住0.3260.4960.034 ：杂项商品和服务0.3960.0960.345特征值5.0981.3520.574贡献率0.6370.1690.072累计贡献率0.6370.8060.878第48页/共54页表7.3.7 按第一主成分排序的31个地区地区地区江西2.2341.867新疆0.6970.647河南1.9470.388四川0.5330.041黑龙江1.9270.636广西0.2512.058吉林1.8

23、590.151山东0.1470.983山西1.8480.404福建0.2011.337内蒙古1.8260.509湖南0.2190.203安徽1.7960.519江苏0.4070.311甘肃1.5490.526云南0.4350.479宁夏1.5010.906西藏0.4372.365辽宁1.3130.844重庆1.1150.409贵州1.2980.341天津2.0060.044海南1.1571.913浙江3.5830.531青海1.0450.426北京5.4262.466陕西0.8590.501广东5.5833.072河北0.7690.580上海5.8660.195湖北0.7170.247第49页

24、/共54页表7.3.8 按第二主成分排序的31个地区地区地区广东5.5833.072山西1.8480.404广西0.2512.058重庆1.1150.409海南1.1571.913青海1.0450.426江西2.2341.867云南0.4350.479福建0.2011.337内蒙古1.8260.509安徽1.7960.519甘肃1.5490.526陕西0.8590.501浙江3.5830.531河南1.9470.388河北0.7690.580贵州1.2980.341黑龙江1.9270.636江苏0.4070.311新疆0.6970.647湖北0.7170.247辽宁1.3130.844湖南0.2190.203宁夏1.5010.906上海5.8660.195山东0.1470.983四川0.5330.041西藏0.4372.365天津2.0060.044北京5.4262.466吉林1.8590.151第50页/共54页第51页/共54页第52页/共54页五、若干补充及应用中需注意的问题1.关于时间序列数据2.主成分用于聚类分析3.关于不同时期的主成分分析4.对综合得分方法的质疑第53页/共54页感谢您的观看！第54页/共54页