《主成分分析法 (2)讲稿.ppt》由会员分享,可在线阅读,更多相关《主成分分析法 (2)讲稿.ppt(32页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、关于主成分分析法 (2)第一页,讲稿共三十二页哦概念:概念:把原来多个变量划为少数几个综合指标的一种把原来多个变量划为少数几个综合指标的一种统计分析方法,是一种降维处理技术统计分析方法,是一种降维处理技术. 主成分分析法主成分分析法 一个研究对象,往往是多要素的复杂系统。变量太多无疑会增加分析问一个研究对象,往往是多要素的复杂系统。变量太多无疑会增加分析问题的难度和复杂性,利用原变量之间的相关关系,用较少的新变量代替原题的难度和复杂性,利用原变量之间的相关关系,用较少的新变量代替原来较多的变量,并使这些少数变量尽可能多的保留原来较多的变量所反应来较多的变量,并使这些少数变量尽可能多的保留原来较
2、多的变量所反应的信息,这样问题就简单化了。的信息,这样问题就简单化了。第二页,讲稿共三十二页哦研究对象指 标x1 x2 xj xp12inx11 x12 x1j x1px21 x22 x2j x2p xi1 xi2 xij xip xn1 xn2 xnj xnp 假设有假设有n个对象,每一个对象都有个对象,每一个对象都有x1,x2,xp个指标构成,它们所个指标构成,它们所对应的要素数据用下表给出:对应的要素数据用下表给出:一、基本原理一、基本原理第三页,讲稿共三十二页哦 原变量为原变量为x1,x2,xp,降维处理后,设它们的综合指标,降维处理后,设它们的综合指标,即新变量为即新变量为z1,z2
3、,z3,zm(mp),则则pmpmmmppppxlxlxlzxlxlxlzxlxlxlz.22112222121212121111第四页,讲稿共三十二页哦z1称为原变量称为原变量x1,x2,xp的第的第一一主成分主成分Z2称为原变量称为原变量x1,x2,xp的第的第二二主成分主成分 .zm称为原变量称为原变量x1,x2,xp的第的第m主成分主成分 找找主成分主成分zi就是要就是要确定系数确定系数lij。从数学上知道,它们从数学上知道,它们分别是分别是x1,x2,xp的的相关系数矩阵相关系数矩阵的的m个较大的个较大的特征特征值值所对应的所对应的特征向量特征向量。第五页,讲稿共三十二页哦二、二、
4、主成分分析的计算步骤主成分分析的计算步骤nknkjkjikinkjkjikiijxxxxxxxxr11221相关系数计算公式相关系数计算公式1、计算相关系数、计算相关系数据公式得这据公式得这p个变量之间的相关系数矩阵为个变量之间的相关系数矩阵为pppppprrrrrrrrrR.212222111211第六页,讲稿共三十二页哦主成分主成分zi的贡献率为的贡献率为 piQi,.21p1kki前前i个主成分的累计贡献率为个主成分的累计贡献率为 piQpkkikk,.21113、计算主成分贡献率及累计贡献率、计算主成分贡献率及累计贡献率 当前当前 i个主成分累计贡献率达到个主成分累计贡献率达到85%9
5、5%,就取前就取前i个主成分作个主成分作为新变量。为新变量。第七页,讲稿共三十二页哦 pjielijiij,.214、计算主成分载荷、计算主成分载荷计算公式为计算公式为得前得前i个主成分在原变量上的载荷个主成分在原变量上的载荷原变量xi主成分Z1Z2Zix1x2xpl11l12l1pl21l22l2pli1li2lip第八页,讲稿共三十二页哦三、主成分分析方法的三、主成分分析方法的SPSS实现实现对45个城市7项经济指标进行主成分分析 1、数据导入到数据窗口中,定义各变量,确保各变量均为数值型。第九页,讲稿共三十二页哦2、激活Analysis 菜单选Data Reduction 的Factor
6、.命令项,弹出Factor Analysis 对话框。在对话框左侧的变量列表中选变量X1 至X7,点击钮使之进入Variables 框。第十页,讲稿共三十二页哦Descriptives.钮变量的均数变量的均数与标准差与标准差相关系数矩阵相关系数矩阵第十一页,讲稿共三十二页哦Extraction.钮:因子提取方法钮:因子提取方法主成分分析法主成分分析法未加权最小平方法未加权最小平方法极大似然估计法极大似然估计法第十二页,讲稿共三十二页哦指定分析矩阵指定分析矩阵提取结果提取结果与因子提取有关的与因子提取有关的输出项输出项迭代次数迭代次数第十三页,讲稿共三十二页哦Rotation.钮:因子旋转方法钮
7、:因子旋转方法旋转的目的是为了获得简单结构,以帮助我们解释因子不作因子旋转不作因子旋转正交旋转正交旋转第十四页,讲稿共三十二页哦Scores.钮:估计因子得分系数的方法钮:估计因子得分系数的方法第十五页,讲稿共三十二页哦D De es sc cr ri ip pt ti iv ve e S St ta at ti is st ti ic cs s614.8580530.2189735.4800.1649435121.157494.6763835840.73171017.738233566.760684.6623635646.6091622.4701335118.4957127.3909335x
8、1x2x3x4x5x6x7MeanStd. DeviationAnalysis N结果解释结果解释均数与标准差均数与标准差第十六页,讲稿共三十二页哦C Co or rr re el la at ti io on n M Ma at tr ri ix x1.000-.341.843.360.404.497.469-.3411.000-.475.309.358.261.311.843-.4751.000.336.324.446.374.360.309.3361.000.941.848.861.404.358.324.9411.000.923.953.497.261.446.848.9231.000
9、.973.469.311.374.861.953.9731.000 x1x2x3x4x5x6x7Correlationx1x2x3x4x5x6x7相关系数矩阵相关系数矩阵K KM MO O a an nd d B Ba ar rt tl le et tt t s s T Te es st t.739313.41721.000Kaiser-Meyer-Olkin Measure of SamplingAdequacy.Approx. Chi-SquaredfSig.Bartletts Test ofSphericityBartlett 检验Bartlett 值值= 313.417,P0.0001
10、,即相关矩阵不,即相关矩阵不是一个单位矩阵,故考虑是一个单位矩阵,故考虑进行因子分析。进行因子分析。第十七页,讲稿共三十二页哦T To ot ta al l V Va ar ri ia an nc ce e E Ex xp pl la ai in ne ed d4.31561.63861.6384.31561.63861.6381.95427.91789.5541.95427.91789.554.3605.13894.692.1852.64497.335.1381.97899.313.033.47399.786.015.214100.000Component1234567Total % of
11、VarianceCumulative %Total % of VarianceCumulative %Initial EigenvaluesExtraction Sums of Squared LoadingsExtraction Method: Principal Component Analysis.特征值、贡献率及累积贡献率特征值、贡献率及累积贡献率第十八页,讲稿共三十二页哦C Co om mp po on ne en nt t M Ma at tr ri ix xa a.613-.686.212.866.549-.775.910.210.954.234.963.096.962.162x
12、1x2x3x4x5x6x712ComponentExtraction Method: Principal Component Analysis.2 components extracted.a. 主成分载荷主成分载荷第十九页,讲稿共三十二页哦因子得分因子得分第一主成分第一主成分第二主成分第二主成分第二十页,讲稿共三十二页哦三、实例三、实例1:对中国大陆对中国大陆31个省(市、区)的第三产业综合发展水个省(市、区)的第三产业综合发展水 平平进行主成分分析。进行主成分分析。1、构建评价第三产业综合发展水平的指标体系,并对原始数据、构建评价第三产业综合发展水平的指标体系,并对原始数据进行标准化处理。
13、这里用于评价第三产业综合发展水平的指标有进行标准化处理。这里用于评价第三产业综合发展水平的指标有人均人均GDP x1 人均第三产业增加值人均第三产业增加值 x2 第二产业增加值比重第二产业增加值比重 x3 第三产业增加值比重第三产业增加值比重 x4 第三产业从业人员比重第三产业从业人员比重 x5 第三产业固定资产投资比重第三产业固定资产投资比重 x6 城市化水平城市化水平 x7第二十一页,讲稿共三十二页哦2、计算各指标之间的相关系数矩阵、计算各指标之间的相关系数矩阵1219. 0839. 0648. 0175. 0808. 0783. 0219. 01003. 0046. 007. 0003.
14、 0001. 0839. 0003. 01712. 0218. 0854. 0822. 0648. 0046. 0712. 01271. 0752. 0657. 0175. 007. 0218. 0271. 01241. 0339. 0808. 0003. 0854. 0752. 0241. 01988. 0783. 001. 0822. 0657. 0339. 0988. 0176543217654321xxxxxxxxxxxxxxRij第二十二页,讲稿共三十二页哦3、计算出相关系数矩阵的特征值、计算出相关系数矩阵的特征值主成分zi特征值i贡献率(%)累计贡献率(%)Z1Z2Z3Z4Z5Z6
15、Z74.1981.2261.0360.2680.1740.09610.0028759.97217.50714.7943.8322.4821.3730.041159.97277.47992.27396.10598.58699.959100从上表知:前三个主成分累计贡献率达从上表知:前三个主成分累计贡献率达92.273%,因此,这三个主成分,因此,这三个主成分Z1、Z2、Z3能够充分反映能够充分反映31个区域第三产业发展的综合水平个区域第三产业发展的综合水平 。第二十三页,讲稿共三十二页哦4、计算主成分载荷、计算主成分载荷 原变量xi主成分载荷lij第一主成分l1i第二主成分l2i第三主成分l3j
16、x1x2x3x4x5x6x70.9460.9710.2200.7950.930-0.07630.8990.1390.02250.945-0.5140.00086-0.220.03320.1230.10.162-0.0910.07020.967-0.187第二十四页,讲稿共三十二页哦 3、第三主成分第三主成分在在第三产业固定资产投资比重第三产业固定资产投资比重x6上有较大的载荷上有较大的载荷,这说明对第三产业的,这说明对第三产业的投资投资是第三产业是第三产业可持续发展可持续发展的重要支撑。的重要支撑。从上表可知:从上表可知: 1、第一主成分第一主成分在在人均人均GDPx1、人均第三产业增加值人均
17、第三产业增加值x2、第三产第三产业从业人员比重业从业人员比重x5、城市化水平城市化水平x7上具有很大的载荷,这些原变量几乎上具有很大的载荷,这些原变量几乎包含了第三产业的主要指标。因此,包含了第三产业的主要指标。因此,第一主成分第一主成分在一定程度上就代表了在一定程度上就代表了第三产业的综合发展水平。第三产业的综合发展水平。 2、第二主成分第二主成分在在第二产业增加值比重第二产业增加值比重x3上有较大的载荷,表明上有较大的载荷,表明第二主成分第二主成分在一定程度上代表在一定程度上代表工业化程度工业化程度。第二十五页,讲稿共三十二页哦5 计算各省区在一二三主成分上的综合得分计算各省区在一二三主成
18、分上的综合得分序号 得分区域第一主成分 L1第二主成分 L2第三主成分L3综合得分L位次1北京3.055-1.1640.8861.75922天津1.6950.418-1.5670.85833上海3.2570.878-0.07322.096123四川-0.7450.1430.571-0.33723第二十六页,讲稿共三十二页哦队员科学成绩智力水平 动手能力写作能力外语水平 协作能力其他特长A8.698.287.99.56B8.28.88.16.57.79.12C88.68.58.59.29.68D8.68.98.39.69.79.78E8.88.48.57.78.69.29F9.29.28.27.
19、9996G9.29.697.29.19.29H789.86.28.79.76I7.78.28.46.59.69.35J8.38.18.66.98.59.44K98.287.899.55L9.69.18.19.98.79.76M9.59.68.38.199.37N8.68.38.28.1995O9.18.78.88.48.89.45P9.38.48.68.88.69.56Q8.489.49.28.49.17R8.78.39.29.18.79.28S7.88.19.67.699.69T98.89.57.97.796第二十七页,讲稿共三十二页哦第二十八页,讲稿共三十二页哦第二十九页,讲稿共三十二页哦第三十页,讲稿共三十二页哦第三十一页,讲稿共三十二页哦感谢大家观看第三十二页,讲稿共三十二页哦