《第三章主成分分析new课件.ppt》由会员分享,可在线阅读,更多相关《第三章主成分分析new课件.ppt(113页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第三章主成分分析第三章主成分分析new第1页,此课件共113页哦主成分分析主成分分析主成分回归主成分回归立体数据表的主成分分析立体数据表的主成分分析第2页,此课件共113页哦 一项十分著名的工作是美国的统计学家斯通(stone)在1947年关于国民经济的研究。他曾利用美国1929一1938年各年的数据,得到了17个反映国民收入与支出的变量要素,例如雇主补贴、消费资料和生产资料、纯公共支出、净增库存、股息、利息外贸平衡等等。1基本思想第3页,此课件共113页哦 在进行主成分分析后,竟以97.4的精度,用三个新变量就取代了原17个变量。根据经济学知识,斯通给这三个新变量分别命名为总收入F1、总收入
2、变化率F2和经济发展或衰退的趋势F3。第4页,此课件共113页哦 主主成成分分分分析析是把各变量之间互相关联的复杂关系进行简化分析的方法。在社会经济的研究中,为了全面系统的分析和研究问题,必须考虑许多经济指标,这些指标能从不同的侧面反映我们所研究的对象的特征,但在某种程度上存在信息的重叠,具有一定的相关性。第5页,此课件共113页哦 主成分分析试图在力保数据信息丢失最少的原则下,对这种多变量的截面数据进行最佳综合简化,也就是说,对高维变量空间进行降维处理。很显然,识辨系统在一个低维空间要比在一个高维空间容易得多。第6页,此课件共113页哦 在力求数据信息丢失最少的原则下,对高维的变量空间降维,
3、即研研究究指指标标体体系系的的少少数数几几个个线线性性组组合合,并并且且这这几几个个线线性性组组合合所所构构成成的的综综合合指指标标将将尽尽可可能能多多地地保保留留原原来来指指标标变异方面的信息。这些综合指标就称为主成分变异方面的信息。这些综合指标就称为主成分。第7页,此课件共113页哦要讨论的问题:(1)基于相关系数矩阵还是基于协方差矩阵做主成分分析。当分析中所选择的经济变量具有不同的量纲,变量水平差异很大,应该选择基于相关系数矩阵的主成分分析。第8页,此课件共113页哦 (2)选择几个主成分。主成分分析的目的是简化变量,一般情况下主成分的个数应该小于原始变量的个数。关于保留几个主成分,应该
4、权衡主成分个数和保留的信息。(3)如何解释主成分所包含的经济意义。第9页,此课件共113页哦2数学模型与几何解释 假设我们所讨论的实际问题中,有p个指标,我们把这p个指标看作p个随机变量,记为X1,X2,Xp,主成分分析就是要把这主成分分析就是要把这p p个指标的问题,转变为讨个指标的问题,转变为讨论论p p个指标的线性组合的问题个指标的线性组合的问题,而这些新的指标F1,F2,Fk(kp),要求保留主要信息量的原则(即充分反映原指标的信息),并且相互独立。第10页,此课件共113页哦这种由讨论多个指标降为少数几个综合指标多个指标降为少数几个综合指标的过程在数学上就叫做降维的过程在数学上就叫做
5、降维。主成分分析通常的主成分分析通常的做法是,寻求原指标的线性组合做法是,寻求原指标的线性组合Fi。第11页,此课件共113页哦满足如下的条件:主成分之间相互独立,即无重叠的信息。即主成分之间相互独立,即无重叠的信息。即主成分的方差依次递减,重要性依次递减,即主成分的方差依次递减,重要性依次递减,即每个主成分的系数平方和为每个主成分的系数平方和为1。即。即第12页,此课件共113页哦主成分分析的几何解释平移、旋转坐标轴第13页,此课件共113页哦主成分分析的几何解释平移、旋转坐标轴第14页,此课件共113页哦主成分分析的几何解释平移、旋转坐标轴第15页,此课件共113页哦主成分分析的几何解释平
6、移、旋转坐标轴第16页,此课件共113页哦为了方便,我们在二维空间中讨论主成分的几何意义。设有n个样品,每个样品有两个观测变量xl和x2,在由变量xl和x2所确定的二维平面中,n个样本点所散布的情况如椭圆状。由图可以看出这n个样本点无论是沿着xl轴方向或x2轴方向都具有较大的离散性,其离散的程度可以分别用观测变量xl的方差和x2的方差定量地表示。显然,如果只考虑xl和x2中的任何一个,那么包含在原始数据中的经济信息将会有较大的损失。第17页,此课件共113页哦 如果我们将xl 轴和x2轴先平移,再同时按逆时针方向旋转角度,得到新坐标轴Fl和F2。Fl和F2是两个新变量。第18页,此课件共113
7、页哦 根据旋转变换的公式:第19页,此课件共113页哦 旋转变换的目的是为了使得n个样品点在Fl轴方向上的离轴方向上的离 散程度最大,即散程度最大,即Fl的方差最大。的方差最大。变量变量Fl代表了原始数据的绝大代表了原始数据的绝大 部分信息,部分信息,在研究某经济问题时,即使不考虑变量F2也无损大局。经过上述旋转变换原始数据的大部分信息集中到Fl轴上,对数据中包含的信息起到了浓缩作用。第20页,此课件共113页哦 Fl,F2除了可以对包含在Xl,X2中的信息起着浓缩作用之外,还具有不相关的性质,这就使得在研究不相关的性质,这就使得在研究复杂的问题时避免了信息重叠所带来的虚假性。二维平复杂的问题
8、时避免了信息重叠所带来的虚假性。二维平面上的样本点的方差大部分都归结在面上的样本点的方差大部分都归结在Fl轴上,而轴上,而F2轴上轴上的方差很小的方差很小。Fl和F2称为原始变量x1和x2的综合变量。F简化了系统结构,抓住了主要矛盾。第21页,此课件共113页哦3 3 主成分的推导及性质主成分的推导及性质 一、两个线性代数的结论一、两个线性代数的结论 1、若A是p阶实对称阵,则一定可以找到正交阵U,使其中 是A A的特征根。第22页,此课件共113页哦 2、若上述矩阵的特征根所对应的单位特征向量为 则实对称阵 属于不同特征根所对应的特征向量是正交的,即有令第23页,此课件共113页哦 二、主成
9、分的推导 (一)(一)第一主成分第一主成分设X的协方差阵为由于x为非负定的对称阵,则有利用线性代数的知识可得,必存在正交阵U,使得第24页,此课件共113页哦 其中1,2,p为x的特征根,不妨假设1 2 p。而U恰好是由特征根相对应的特征向量所组成的正交阵。下面我们来看,是否由U的第一列元素所构成为原始变量的线性组合是否有最大的方差。第25页,此课件共113页哦设有P维正交向量第26页,此课件共113页哦当且仅当 时,即 时,有最大的方差 。因为第27页,此课件共113页哦如果第一主成分的信息不够,则需要寻找第二主成分。第28页,此课件共113页哦(二)(二)第二主成分第二主成分如果第一主成分
10、的信息不够,则寻找第二主成分 因为第29页,此课件共113页哦 所以如果取线性变换:则 的方差次大。类推 第30页,此课件共113页哦写为矩阵形式:第31页,此课件共113页哦4 4 主成分的性质主成分的性质一、方差为所有特征根之和一、方差为所有特征根之和 说明主成分分析把P个随机变量的总方差分解成为P个不相关的随机变量的方差之和。协方差矩阵的对角线上的元素之和等于特征根之和。第32页,此课件共113页哦 二、精度分析二、精度分析 1)贡献率贡献率:第i个主成分的方差在全部方差中所占比重 ,称为贡献率,因为,第一主成分的方差最大,其贡献率也最大,说明它综合原来P个指标的信息的能力最强,其它主成
11、分依次渐弱。2)累积贡献率累积贡献率:前k个主成分共有多大的综合能力,用这k个主成分的方差和在全部方差中所占比重来描述,称为累积贡献率。第33页,此课件共113页哦 我们进行主成分分析的目的之一是希望用尽可能少的主成分F1,F2,Fk(kp)代替原来的P个指标。到底应该选择多少个主成分,在实际工作中,主成分个数的多少取决于能够反映原来变量80%以上的信息量为依据,即当累累积积贡贡献献率率80%80%时时的的主主成成分分的的个个数数就就足足够够了了。最常见的情况是主成分为2到3个。第34页,此课件共113页哦三、原始变量与主成分之间的相关系数变量与主成分之间的相关系数第35页,此课件共113页哦
12、 可见,和 的相关的密切程度取决于对应线性组合系数的大小。第36页,此课件共113页哦第37页,此课件共113页哦四、原始变量被主成分的提取率原始变量被主成分的提取率 前面我们讨论了主成分的贡献率和累计贡献率,它度度量量了了F F1 1,F F2 2,F Fm m分分别别从从原原始始变变量量X X1 1,X X2 2,XXP P中中提提取取了了多多少少信信息息。那那么么X X1 1,X X2 2,XXP P各各有有多多少少信信息息分分别别被被F F1 1,F F2 2,F Fm m提提取取了了。应该用什么指标来度量?我们考虑到当讨论F1分别与X1,X2,XP的关系时,可以讨论F1分别与X1,X
13、2,XP的相关系数,但是由于相关系数有正有负,所以只有考虑相关系数的平方。第38页,此课件共113页哦 如果我们仅仅提出了m个主成分,则第i 原始变量信息的被提取率为:是Fj 能说明的第i 个原始变量的方差是Fj 提取的第i 个原始变量信息的比重第39页,此课件共113页哦 例例:设 的协方差矩阵为 解得特征根为 ,第 一 个 主 成 分 的 贡 献 率 为5.83/(5.83+2.00+0.17)=72.875%,尽管第一个主成分的贡献率并不小,但在本题中第一主成分不含第三个原始变量的信息,所以应该取两个主成分。第40页,此课件共113页哦Xi与F1的相关系数平方Xi与F2的相关系数平方信息
14、提取率xi10.9250.855000.8552-0.9980.996000.996300110第41页,此课件共113页哦 定定义义:如果一个主成分仅仅对某一个原始变量有作用,则称为特殊成分。如果一个主成分所有的原始变量都起作用称为公共成分。(该题无公共因子)第42页,此课件共113页哦5 5 主成分分析的步骤主成分分析的步骤 在实际问题中,X的协方差阵通常是未知的,样品有 则样本协差阵为:第 一 步:由X的 协 方 差 阵x,求 出 其 特 征 根,即 解 方 程 ,可得特征根 。一、基于协方差矩阵第43页,此课件共113页哦 第二步:求出特征根分别对应的特征向量U1,U2,Up,第三步:
15、计算累积贡献率,给出恰当的主成分个数。第四步:综合评价。计算所选出的k个主成分的得分。将原始数据代入前k个主成分的表达式,分别计算出各单位k个主成分的得分,以方差贡献率为权数,求得k个主成分的得分的加权平均数,并按得分值的大小排队评价。第44页,此课件共113页哦 二、基于相关系数矩阵 如果变量有不同的量纲,则必须基于相关系数矩阵进行主成分分析。不同的是计算得分时应采用标准化后的数据。第45页,此课件共113页哦 例一例一 应收账款是指企业因对外销售产品、材料、提供劳务及其它原因,应向购货单位或接受劳务的单位收取的款项,包括应收销货款、其它应收款和应收票据等。出于扩大销售的竞争需要,企业不得不
16、以赊销或其它优惠的方式招揽顾客,由于销售和收款的时间差,于是产生了应收款项。应收款赊销的效果的好坏,不仅依赖于企业的信用政策,还依赖于顾客的信用程度。由此,评价顾客的信用等级,了解顾客的综合信用程度,做到“知己知彼,百战不殆”,对加强企业的应收账款管理大有帮助。某企业为了了解其客户的信用程度,采用西方银行信用评估常用的5C方法,5C的目的是说明顾客违约的可能性。第46页,此课件共113页哦 1、品格(用X1表示),指顾客的信誉,履行偿还义务的可能性。企业可以通过过去的付款记录得到此项。2、能力(用X2表示),指顾客的偿还能力。即其流动资产的数量和质量以及流动负载的比率。顾客的流动资产越多,其转
17、化为现金支付款项的能力越强。同时,还应注意顾客流动资产的质量,看其是否会出现存货过多过时质量下降,影响其变现能力和支付能力。3、资本(用X3表示),指顾客的财务实力和财务状况,表明顾客可能偿还债务的背景。4、附带的担保品(用X4表示),指借款人以容易出售的资产做抵押。5 5、环境条件(用X5表示),指企业的外部因素,即指非企业本身能控制或操纵的因素。第47页,此课件共113页哦 首先抽取了10家具有可比性的同类企业作为样本,又请8位专家分别给10个企业的5个指标打分,然后分别计算企业5个指标的平均值,如表。76.581.57675.871.78579.280.384.476.570.67367
18、.668.178.5949487.589.59290.787.39181.58084.666.968.864.866.477.573.670.969.874.857.760.457.460.86585.668.57062.276.57069.271.764.968.9;第48页,此课件共113页哦 Total Variance=485.31477778 Eigenvalues of the Covariance Matrix Eigenvalue Proportion Cumulative PRIN1 410.506 0.845854 0.84585 PRIN2 43.264 0.089146
19、 0.93500 PRIN3 20.670 0.042591 0.97759 PRIN4 8.071 0.016630 0.99422 PRIN5 2.805 0.005779 1.00000 Eigenvectors PRIN1 PRIN2 PRIN3 PRIN4 PRIN5 X1 0.468814 -.830612 0.021406 0.254654 -.158081 X2 0.484876 0.329916 0.014801 -.287720 -.757000 X3 0.472744 -.021174 -.412719 -.588582 0.509213 X4 0.461747 0.43
20、0904 -.240845 0.706283 0.210403 X5 0.329259 0.122930 0.878054 -.084286 0.313677第49页,此课件共113页哦 第一主成份的贡献率为84.6%,第一主成份 Z1=0.469X1+0.485X2+0.473X3+0.462X4+0.329X5 的各项系数大致相等,且均为正数,说明第一主成份对所有的信用评价指标都有近似的载荷,是对所有指标的一个综合测度,可以作为综合的信用等级指标。可以用来排序。将原始数据的值标准化后,代入第一主成份Z1的表示式,计算各企业的得分,并按分值大小排序:在正确评估了顾客的信用等级后,就能正确制定
21、出信用期、收帐政策等,这对于加强应收帐款的管理大有帮助。序号序号1 12 23 34 45 56 67 78 89 91010得分得分3.163.1613.613.6-9.01-9.0135.935.925.125.1-10.3-10.3-4.364.36-33.8-33.8-6.416.41-13.8-13.8排序排序4 43 37 71 12 28 85 510106 69 9第50页,此课件共113页哦 根据主成分分析的定义及性质,我们已大体上能看出主成分分析的一些应用。概括起来说,主成分分析主要有以下几方面的应用。1主成分分析能降低所研究的数据空间的维数主成分分析能降低所研究的数据空间
22、的维数。即用研究m维的Y空间代替p维的X空间(mp),而低维的Y空间代替高维的x空间所损失的信息很少。即使只有一个主成分Yl(即 m1)时,这个Yl仍是使用全部X变量(p个)得到的。例如要计算Yl的均值也得使用全部x的均值。6主成分分析主要有以下几方面的应用第51页,此课件共113页哦 2有时可通过因子负荷aij的结构,弄清X变量间的某些关系。3.多维数据的一种图形表示方法。我们知道当维数大于3时便不能画出几何图形,多元统计研究的问题大都多于3个变量。要把研究的问题用图形表示出来是不可能的。然而,经过主成分分析后,我们可以选取前两个主成分,根据主成分的得分,画出n个样品在二维平面上的分布况,由
23、图形可直观地看出各样品在主分量中的地位。第52页,此课件共113页哦 4由主成分分析法构造回归模型。即把各主成分作为新自变量代替原来自变量x做回归分析。5用主成分分析筛选回归变量。回归变量的选择有着重要的实际意义,为了使模型本身易于做结构分析、控制和预报,好从原始变量所构成的子集合中选择最佳变量,构成最佳变量集合。用主成分分析筛选变量,可以用较少的计算量来选择量,获得选择最佳变量子集合的效果。第53页,此课件共113页哦主成分回归介绍主成分回归介绍第54页,此课件共113页哦 国际旅游外汇收入是国民收入是国民经济发展的国际旅游外汇收入是国民收入是国民经济发展的重要组成部分,影响一个国家或地区旅
24、游收入的因素重要组成部分,影响一个国家或地区旅游收入的因素包括自然、文化、社会、经济、交通等多方面的因素。包括自然、文化、社会、经济、交通等多方面的因素。中国统计年鉴中国统计年鉴把第三次产业划分为把第三次产业划分为12个组成部个组成部分,分别为:分,分别为:一、提出问题一、提出问题第55页,此课件共113页哦x1:农林牧渔服务业:农林牧渔服务业 x2:地质勘查水利管理业:地质勘查水利管理业x3:交通运输仓储和邮电通讯业:交通运输仓储和邮电通讯业 x4:批发零售贸易和餐食业:批发零售贸易和餐食业x5:金融保险业:金融保险业 x6:房地产业:房地产业 x7:社会服务业:社会服务业 x8:卫生体育和
25、社会福利业:卫生体育和社会福利业 x9:教育文艺和广播:教育文艺和广播 x10:科学研究和综合艺术:科学研究和综合艺术x11:党政机关:党政机关 x12:其他行业:其他行业 选自选自1998年我国年我国31个省、市、自治区的数据。以旅游外汇收入个省、市、自治区的数据。以旅游外汇收入(百万美圆)为因变量。自变量的单位为亿元人民币。数据略。(百万美圆)为因变量。自变量的单位为亿元人民币。数据略。第56页,此课件共113页哦InterceptCoefficients-205.236116.8459-1.756460.096008标准误差tStatP-valueXVariable1-1.4004522
26、.8676-0.061240.951842XVariable22.67500118.575080.144010.887092XVariable33.3008772.4645561.3393390.197128XVariable4-0.944021.296117-0.728340.475774XVariable5-5.50164.508593-1.220250.238117XVariable64.0544343.9537451.0254670.318728XVariable74.1425.0699840.8169650.42463XVariable8-15.364910.82589-1.4192
27、70.172905XVariable917.367668.353372.0791210.052178XVariable109.07888310.147280.8947110.38275XVariable11-10.585.610696-1.885690.075582XVariable121.3507095.0015040.270060.790186这个模型是不理想的,一个最严重的问题是多重共线性的问题。第57页,此课件共113页哦线性回归模型的方差分析表方差来源自由度离差平方和方差F统计量显著性水平回归分析1211690140974178.310.51335 8.15025E-06残差1816
28、6789992661.04总计3113358039 利用主成分的互不相关性来建立因变量与主成分的回归,在理论上可以达到消除多重共线性。第58页,此课件共113页哦 二、主成分回归方法二、主成分回归方法第59页,此课件共113页哦原始数据观测矩阵主成分系数矩阵第60页,此课件共113页哦主成分得分矩阵第61页,此课件共113页哦根据最小二乘估计,则基于协方差矩阵的主成分回归基于相关系数矩阵的主成分回归第62页,此课件共113页哦主成分回归系数的协方差矩阵第63页,此课件共113页哦1、经济分析数据Y:进口总额X1:GDPX2:积累总额X3:消费总额求进口总额与GDP、积累总额和消费总额之间的回归
29、方程。三、主成分回归的实例第64页,此课件共113页哦dataa;inputx1-x3y;cards;149.34.2108.115.9161.24.1114.816.4171.53.1123.219.0175.53.1126.919.1180.81.1132.118.8190.72.2137.720.4202.12.1146.022.7212.45.6154.126.5226.15.0162.328.1231.95.1164.327.6239.00.7167.626.3;procregoutest=b;modely=x1-x3/pcomit=1,2outvif;procprintdata=b
30、;procstandarddata=aout=cmean=0std=1;varx1-x3y;procprincompdata=cout=dprefix=z;varx1-x3;procregdata=d;modely=z1z2/noint;run;第65页,此课件共113页哦Analysis of Variance Sum of Mean Source DF Squares Square F Value Pr F Model 3 204.77614 68.25871 285.61|t|Intercept 1 -10.12799 1.21216 -8.36 F Model 2 9.88278 4.
31、94139 379.38|t|F1 1 0.68998 0.02552 27.03 F Model 6 498504 83084 47.22|t|Intercept 1 386505 122516 3.15 0.0116 x1 1 13.71162 11.68424 1.17 0.2707 x2 1 0.00846 0.00461 1.84 0.0995 x3 1 0.09405 0.06720 1.40 0.1952 x4 1 0.20562 0.02948 6.97 .0001 x5 1 -0.00435 0.03111 -0.14 0.8918 x6 1 -199.20213 62.67
32、100 -3.18 0.0112第77页,此课件共113页哦 浅谈时序立体数据的主成分分析第78页,此课件共113页哦前面介绍的主成分分析方法,成功地实现了截面数据的最佳综合和简化。然而,在现实生活中,随着时间的发展与数据的积累,人们开始拥有大量按时间顺序排列的平面数据表序列,这样一组按时间顺序排放的数据表序列就像一个数据匣,被称为时序立体数据表。本章将介绍如何对这种多维动态数据系统进行立体式的综合简化,并在此基础上,迅速提取立体数据表中的重要信息,充分发掘其中的丰富内涵,从而简化扼要地把握系统的动态规律。第79页,此课件共113页哦第一节全局分析的概念时序立体数据表时一个按时间顺序排放的数据
33、表序列。如果对每一张数据表分别进行主成分分析,则不同的数据表有完全不同的简化空间,就无法保证系统分析的统一性、整体性和可比性。因此,对这种数据表进行主成分分析,得到一个统一的简化子空间。一、全局概念假设有个样本,个指标,时间的跨度为。时序立体数据表,第80页,此课件共113页哦若以为变量的指标,在时刻数据表中对上列数据的分析称为全局分析。第81页,此课件共113页哦二、全局变量 全局群点在j指标上的取值分布被称为全局变量,表示为三、全局重心 全局数据表的重心为 第82页,此课件共113页哦 权数应该根据不同时刻的重要性来决定,也可以等权,等权时,均值为:时刻t的数据表重心为四、全局方差全局变量
34、的方差:第83页,此课件共113页哦五、全局协方差全局变量的协方差为:全局协方差矩阵:第84页,此课件共113页哦第二节第二节 全局主成分分析全局主成分分析 一、全局主成分分析的步骤为(1)求全局相关系数矩阵 (2)求 的特征根不妨假设 和对应的特征向量:第85页,此课件共113页哦第86页,此课件共113页哦第三节对经典主成分分析的继承性一、全局主成分一定对应于数据变易最大的方向二、全局主成分是对原始变量系统的最佳综合在全局主成分分析中,还可以证明,若全局数据表种有p个变量,如果想以一个综合变量来取代原来所有的全局变量,则第一个主成分F1就是最好的选择。第87页,此课件共113页哦这个结论可
35、以推广到m维空间:第88页,此课件共113页哦三、全局分析与单张数据表分析的联系设j(j=1,2,m)是全局特征值 (j=1,2,m)是第t时刻的数据表所计算的特征值第89页,此课件共113页哦上式反映了全局第h个主成分与单张数据表个主成分之间的数量关系。特别当h=1时:因此,如果各年数据表的重心在第一主成分上的投影不发生改变,则说明,第一主成分与单张数据表的主成分之间最相关。第90页,此课件共113页哦第四节精度分析 一、全局精度 以数据变异的大小来恒量数据中的信息量如果变量已经被标准化,则精度为:第91页,此课件共113页哦二、数据表Xt的表现精度数据表Xt的表现精度是指群点在全局主成分上
36、的近似精度。令是第t张表中的第i个样本在全局第h个主成分的得分。第92页,此课件共113页哦第93页,此课件共113页哦第五节数据主要特征的动态分析为了迅速把握多维动态数据群种的主要信息,还应该对数据系统的主要特征进行动态分析研究。数据群点有如下特征:(1)的总体水平(2)的主轴(3)的主轴上的分布偏差(4)中各样本点间的相对位置和排列顺序。第94页,此课件共113页哦一、总体水平第t年数据群点的总体水平为。可以从三个方面研究其动态数据信息。(1)的时序轨迹(2)对于1一p个变量指标,研究哪一个指标在1一T年间发生的变化最大。首先,j指标在1一T年间的变化可以用aj表示,有第95页,此课件共1
37、13页哦所有指标在1T年的变化为a表示,有使cj最大的指标xj,在1T年发生的变化最大,在经济系统分析中,过大过小的cj都应是分析人员关注的对象。(3)从1T年,研究在哪一年发生了较大的变化。这是比值,比cj更加深入的分析。第96页,此课件共113页哦则说明j指标在tt+1年间的变化比其它年间更大。第97页,此课件共113页哦二、主轴对第t年的数据表xt做平面主成分分析,可以得到一组主轴,对应的有特征值,分析是如何随时间变化的,可以了解数据的主要特征发展变化的历史过程。第98页,此课件共113页哦从前面的分析可以知道,是第t年数据变异最大的方向,数据在这个方向被拉得最长。如果研究国民生活水平的
38、话,则在这一方向人们生活水平的差距最大,所以,是最能反映国民生活水平的主要特征。与对应的是主成分。数据的主要特征随时间的发展会发生变化,这个变化可以通过的变化过程来观察。特别对于第一、第二主轴(即h1,2),以及后续含数据信息量较大的那些主轴,更应给予重点研究。第99页,此课件共113页哦三、方差的变化在数据表由x1,x2,xT的变化过程中,除了需要研究数据的主要特征随时间的变化以外,还要分析数据在主轴上的分布方差是否发生了较大的变化。分别从以下三个指标来观察数据在主轴散布范围发生的变化。第100页,此课件共113页哦(1)在h轴上,数据的分散程度的差分(2)比较在t+1年,哪个主轴的散布范围
39、较大(3)比较1T年间,哪个主轴的分散范围较大第101页,此课件共113页哦四、样本点间相对位置和排列顺序的变化随着时间的发展,群点在某一方向上的相对位置和排列顺序也会发生变化。例如,改革开放以来,我国沿海城市经济发展速度较其他地区的城市要快,特别在对外贸易方面,其发展更为显著。如果第一主轴反映了城市经济的综合实力,则在这个轴上可以看出,在不同的年份上,各城市由于发展速度不一,因此,相对位置和顺序都有变化,沿海城市的经济实力显然日趋向前。第102页,此课件共113页哦如何反映样本点间位置和顺序的变化呢?有一个要点必须注意,这就是必须在同一的轴上比较样本点的位置和顺序,因此,取全局主成分分析的第
40、h主轴,它对所有时刻的数据表都是同一的。在其上的投影为1、在上的投影坐标是否有明显移动第103页,此课件共113页哦2、样本点排列顺序的改变第104页,此课件共113页哦下例是我国1998年和1999年城镇居民分地区的消费支出资料:X1:食品支出X2:衣着支出X3:家庭设备用品及服务支出X4:医疗保健支出X5:交通和通讯支出X6:娱乐教育文化支出X7:居住支出X8:杂项商品支出进行主成分分析,并比较全局主成分分析和单张数据表主成分分析的结果。第105页,此课件共113页哦EigenvaluesoftheCorrelationMatrix(全局主成分特征根)EigenvalueDifferenc
41、eProportionCumulativeA16.991256.443290.8739060.87391A20.547960.395310.0684950.94240A30.152660.030190.0190820.96148A40.122470.039720.0153090.97679A50.082750.020420.0103440.98714A60.062330.021900.0077920.99493A70.040440.040300.0050550.99998A80.00014.0.0000181.00000第106页,此课件共113页哦 全局主成分特征向量全局主成分特征向量 A1
42、 A2 A3 A4 A5 A6 A7 A8 X1 0.374493 -.172257 0.030143 0.136213 0.076849 0.062345 0.005073 -.894875 X2 0.346007 -.445411 0.024956 0.532852 0.438070 0.136731 -.241623 0.358262 X3 0.311984 0.710728 0.411674 0.164345 0.360232 -.055544 0.253882 0.061138 X4 0.362343 -.194425 0.293868 0.105955 -.623604 0.310
43、654 0.461967 0.185796 X5 0.360705 -.096981 0.484438 -.536079 -.081534 -.154719 -.548523 0.083447 X6 0.345751 0.425463 -.516114 -.016448 -.203576 0.488600 -.386444 0.057272 X7 0.364743 0.060889 -.332168 0.215494 -.291244 -.784080 0.002187 0.082907 X8 0.358775 -.186733 -.362278 -.570254 0.388105 0.019
44、425 0.462062 0.124385 第107页,此课件共113页哦EigenvaluesoftheCorrelationMatrix98年数据表的主成分分析EigenvalueDifferenceProportionCumulativeB17.105926.589490.8882400.88824B20.516430.391980.0645530.95279B30.124440.024300.0155550.96835B40.100140.023200.0125170.98087B50.076940.020310.0096170.99048B60.056620.037210.00707
45、80.99756B70.019420.019320.0024270.99999B80.00010.0.0000121.00000第108页,此课件共113页哦 Eigenvectors98年数据表的主成分分析 B1 B2 B3 B4 B5 B6 B7 B8 X1 0.372150 -.159966 -.071551 -.057458 0.102394 0.118105 -.006808 -.896111 X2 0.349028 -.418593 -.335008 -.152609 0.517027 0.310746 0.279572 0.354117 X3 0.312789 0.729505
46、0.271016 -.034411 0.496477 0.150524 -.151145 0.057958 X4 0.365701 -.120455 0.051476 -.389966 -.416759 0.320595 -.618933 0.193606 X5 0.361312 -.092321 0.626102 -.266618 -.242140 -.169833 0.551419 0.079394 X6 0.347155 0.403651 -.479912 0.280454 -.487565 0.195188 0.358038 0.059746 X7 0.364365 0.038843
47、-.325100 -.154922 0.069968 -.833079 -.175368 0.079806 X8 0.352541 -.282022 0.280252 0.803997 0.027822 -.055566 -.229918 0.120426第109页,此课件共113页哦EigenvaluesoftheCorrelationMatrix(99年数据表的主成分分析)EigenvalueDifferenceProportionCumulativeC16.943786.340700.8679730.86797C20.603080.443010.0753850.94336C30.1600
48、80.042350.0200100.96337C40.117730.029250.0147160.97808C50.088480.035560.0110610.98914C60.052920.019150.0066150.99576C70.033770.033610.0042210.99998C80.00015.0.0000191.00000第110页,此课件共113页哦 Eigenvectors99年数据表的主成分分析 C1 C2 C3 C4 C5 C6 C7 C8 X1 0.375326 -.176598 -.025298 0.142244 0.043048 0.053525 0.0552
49、25 -.894037 X2 0.342405 -.456571 -.142611 0.634127 0.299171 0.174451 -.033774 0.361620 X3 0.311094 0.687727 0.332464 0.410055 0.013385 -.119860 0.364141 0.066550 X4 0.360191 -.263051 0.143639 -.224628 -.638715 0.343534 0.414770 0.178802 X5 0.359409 -.098416 0.703325 -.237833 0.141298 -.042767 -.5301
50、38 0.084078 X6 0.345183 0.435106 -.413152 -.254514 0.167365 0.591044 -.274788 0.056631 X7 0.365788 0.080977 -.388392 0.085724 -.476896 -.561740 -.389637 0.081534 X8 0.365029 -.112389 -.177901 -.479953 0.474450 -.409397 0.427557 0.128860第111页,此课件共113页哦第112页,此课件共113页哦F1F2F39898关于全局主成分关于全局主成分的解释方差的解释方差