多元统计分析课程设计

主成分分析法在我国居民生活质量状况

综合评价中的应用内容摘要:改革开放以来，我国各地区间的经济发展速度有着明显差别，而人民的生活质量也因此产生了不同，本文用主成分分析法，选取多个指标，对全国31个省市居民的生活质量进行了简单的分析。关键词：数据选取数据分析主成分分析使用软件：SPSS一主成分分析主成分分析定义主成分分析也称主分量分析，旨在利用降维的思想，把多指标转化为少数几个综合指标。在实证问题研究中，为了全面、系统地分析问题，我们必须考虑众多影响因素。这些涉及的因素一般称为指标，在多元统计分析中也称为变量。因为每个变量都在不同程度上反映了所研究问题的某些信息，并且指标之间彼此有一定的相关性，因而所得的统计数据反映的信息在一定程度上有重叠。在用统计方法研究多变量问题时，变量太多会增加计算量和增加分析问题的复杂性，人们希望在进行定量分析的过程中，涉及的变量较少，得到的信息量较多。主成分分析法方法简介主成分分析法是一种数学变换的方法，它把给定的一组相关变量通过线性变换转成另一组不相关的变量，这些新的变量按照方差依次递减的顺序排列。在数学变换中保持变量的总方差不变，使第一变量具有最大的方差，称为第一主成分，第二变量的方差次大，并且和第一变量不相关，称为第二主成分。依次类推，I个变量就有I个主成分。其中Li为p维正交化向量(Li*Li=1)，Zi之间互不相关且按照方差由大到小排列，则称Zi为X的第I个主成分。设X的协方差矩阵为S，则S必为半正定对称矩阵，求特征值入i(按从大到小排序)及其特征向量，可以证明，入i所对应的正交化特征向量，即为第I个主成分Zi所对应的系数向量Li，而Zi的方差贡献率定义为Xi/2Xj,通常要求提取的主成分的数量k满足2Xk/SXj>0.85o主成分分析主要目的是希望用较少的变量去解释原来资料中的大部分变异，将我们手中许多相关性很高的变量转化成彼此相互独立或不相关的变量。通常是选出比原始变量个数少，能解释大部分资料中的变异的几个新变量，即所谓主成分，并用以解释资料的综合性指标。由此可见，主成分分析实际上是一种降维方法。二问题背景及数据问题及背景背景：随着生产力水平的不断提高，我国居民生活水平不断提高，生活质量也在不断改善。但是，受各地生产力发展水平不平衡的影响，我国各地居民的生活质量也表现为不平衡。

问题：利用主成分分析法对我国31个省市、自治区居民的生活状况进行评价分析。为全面分析各地居民生活状况，可选取如下指标体系进行反应：职工人均工资、人均居住面积、城市人均用水普及量、城市煤气普及量、人均拥有道路面积、人均绿地公共面积、批发零售贸易商品销售总额、旅游外汇收入。数据地区居住面理人均用水煤旦昔茂人均道路商品急质施溢外汇北京56323节i70nnrnirrt=21258324044S天津41"28.31^00101CC.ZCU.339900401ZC1河北247363071995797.1114.499.4939765027』LJ3艺二E2C5293.27F375390ZC1闰炭口归1U2L4782.037425127E11.1U2127.9057717宁27729M酒耶「日9；-「9"092700r%n功2348GJIU4洗上3W.二1-.JJJ.JU304040呈职2324B21.728J2479.459.2B346227640870二溥532CCC230^00101CC.ZC4.】7.0229712.504：72江节3IEE7110699.E8ga.2322.23I3.1I20643.203EEC浙江34何GO4399.7037.7212.239.6010270000发侬^36329.8896.11a?.EO14.159.293765.40心福莲JR397397-1z1511F574340E21-LL3'5bUb/ULL.1t11.JoIJ.bU1340301.1*行「£qq-q□F=rPnlM?n177hRH|-cl河南24Z1C31CD65.56EE二1二刀3.204403302273939.0497.as90.9013.D39.40E183.E0JJ3混市24870Ml7?叫「7「』汗1zn79G263030717厂?JJ11L2789ii士11..lb22348.80厂匹2SEE0317592.a?日40411E3B.61199860ECQ海南21ZC^22C4C3.C772.Z112.J5100704CO二1』重&^0536.0393.20go.879.499.622891.20460四JI2jZ2l349466.09l1.ZS107Bj.74g701口员刊246OE25.2788.6967.826.22B.161076.50117云南27449n；77Fir1/V53075801TF西抻4FJzLJJU；86.59a.zL1z.4o□,b4b4'Llclnn96.66F匚w17F"IMFFzn=F「甘而24217190787"EU侦G001C3333319.78inn.oogj.7811.IBB.63186.9010一点IT19qnnT1dpnnn48930「』土/2278tt.zl7.91236330let图1数据来源：《中国统计年鉴2009》三详细分析过程第一步：录入数据，有以下变量：职工人均工资，人均居住面积，城市人口用水普及量，城市煤气普及量，人均拥有道路面积，人均绿地公共面积，批发零售贸易商品销售总额，旅游外汇收入，见图2

NameT憎WidihDecimalsLabelValues-MiaBingColumnsAlign地区String80NoneNoneBSLeft净5.JuJ..3I工住Numerica0FMoneNoneB=Left疔，血Numerica2NoneNoneB=Righl/Scale人上用木Numeric82NoneNoneB=Right#Scale:miNumenca2NoneNoneB=Right疫Scale人切道路Numeric82NoneNoneB=Right疥Sc;lc.j蜡地NumEriE:NoneNoneBMRighi:/Scale商品孔矛Numeric3zNoneNoneE三Fi；l/5.Ju旅涛外FNumericaDNoneNoneBSRight#ScaleV/raO：O2Numeric82NoneNoneB三Riqhl疳3.Ju第二步：选择功能模块MewQshIrensformAnalyzeGraphsLHittir?Add-oniHelph■f□Em®[Retirris►件：tDescriptKtStatistics►地区Tsbfe,煤气普及人均值暗人岫弛商品总辙旅泳外汇北京5心CorripareMesns>l.UUJ•一£UGJ2E632.404心r老J17JBLBTfeT：日LJn^HTMDCltd>Gener：日l歹1LinearModels,M咎dModeis►Correlarte►Regresson►Lo^inear►1QDDO14.39ESDD.401QD1可北24766iM.d?9■!?397E.5017125B2S:2b4:；b4:；5:.jU:Ul内奈古26114叫mI；兀1'-1:-KF丁.二=277^9923B9.95937ES27.B0152623436NeuralTJetELorks►34J321D.399.2D丑I4H和2115：_23MGCbMify►出k9464U:,'U上海569G5QalaReducllon►\Eaclor..70229712.0049讶-玉.31BG7scale>&CorresponderceAna^sla..13.112C643.20^80淅江<?微汕纯;i;:::Nanpsrmetnclests*|T[meSeries>fcjj□pllmalSeeling..~s?cn伺1二UbJil.'Jl:2,'...UJ.245串:Surviv-ai*=7W|;-=、in42-74=JO2-A4•二亏21TOOMissingVeiusAtsI^se...90.1B11.C610.6D134D.30252±KJS4_4M或祖日RBBfMnse►.a5jIj.zUUJJllr,'7.：UIzUl问南CorrplDKSornplK►■A91顷n?i御-FF缶.L22739fiudHyContrd►90.9D13.039.4DB183.E0J43湖南24B70\FROCCurjic..日4展12.01798353B.30617第三步：将变量添加到Varicribles中第四步：输入信息地区]SignificancelevelsReproduced[KMOandBartlett'stestolsphericityCorrelationMatrixScores...Ojjtions...DeterrriinantIAnti-irriageDeacriptives]CoefficierrtsInverseCancelUnivariatedeacriptives叵Lnrtialscl'.-rtionExtrsution..Rotation...ContinuerE3<tractnAnalyzeMaxirriumItersrtiQnaforConwergence:CorrelationmatrixC：Covariancemalrix回UnrutateclfactorsolLJtion「巳巳尸|匚计<*；：Eigenvaluesover:QNurriber-o-f-factors:CancelContinueMetnod:PrincipalccmiFjonentaDe^criptives...Extraction...Rut_ation...Scores...Cyjtions...PromaxKappa[4DirectOblirTiinDisplay0Rotatedsoh-rtion回Lpadicgi'*)^arimax•,J'EquarnaxMaximumIterationsforConvergence:Cancel歹地区MethodNoneQuartirnaxContinue~SacIoi-AjL^lysie:Rotation图8图9第五步：单击“OK”按钮，完成运算。图10三.结果分析1.CommunalitiesInitialExtraction人均工资1.000.730居住面积1.000.598人均用水1.000.636煤气普及1.000.794人均道路1.000.776人均绿地1.000.771商品总额1.000.883旅游外汇1.000.653分析：第一列是列出八个原始变量，第二列是根据主成分分析初始解计算出变量共同度，第三列是是根据主成分分析最终解计算出变量共同度，这时由于因子变量个数少于原始变量个数，因此每个变量的共同度必然小于1。例如，第一行中0.730表示m个因子变量共同解释掉原始变量“人均工资”方差72.2%。2.TotalVarianceExplainedComponentInitialEigenvaluesExtractionSumsofSquaredLoadingsRotationSumsofSquaredLoadingsTotal%ofVarianceCumulative%Total%ofVarianceCumulative%Total%ofVarianceCumulative%13.95549.43849.4383.95549.43849.4383.87448.41948.41921.88623.58173.0191.88623.58173.0191.96824.60073.0193.88311.03384.0514.4635.78389.8345.4425.53095.3646.2302.87498.2397.0931.16499.4038.048.597100.000ExtractionMethod:PrincipalComponentAnalysis.分析：上表为spss输出结果，从上表可以看出特征值和和贡献率。从上表可以看出公共因子对原变量总体的描述情况。可以看出前2个公共因子的的贡献

率达到73.019%，所以提取2个公共因子就可以反映原变量的大部分信息。3.ScreePlotIIIIIIIIIII12345678ComponentNumber分析：上图为公共因子碎石图，它的横坐标为公共因子数，纵坐标为公共因子的特征值。可以看出前2个公共因子的特征值变化非常明显，到2个以后趋于平稳。所以得出提取2个公共因子可以对原变量的信息描述有显著作用。这与Communalities的结论也相符合。4.ComponentScoreCoefficientMatrixComponent12人均工资.216-.272居住面积.200-.010人均用水.173.134煤气普及.187.170人均道路-.068.460人均绿地.018.426商品总额.249-.064|旅游外汇|.213|-.040|分析：该表格是因子得分矩阵，这是根据回归年算法计算出来的因子得分函数的系数，根据这个表格可得下面的因子得分函数F1=0.216x1+0.200x2+0.173x3+0.187x4-0.068x5+0.018x6+0.249x7+0.213x8F2=-0.272x1-0.010x2+0.134x3+0.170x4+0.460x5+0.426x6-0.064x7-0.040x8SPSS将根据2个因子得分函数自动计算样本的2个因子得分，并且2个因子作为新变量，保存到SPSS窗口中第一主成分在人均拥有道路面积的系数上为负，其他为正，而且职工人居工资、人均居住面积、批发零售贸易商品销售总额、旅游外汇收入的系数绝对值比较大，说明第一主成分代表了我国居民生活质量状况针对职工人居工资、人均居住面积、批发零售贸易商品销售总额、旅游外汇收入和其他居民生活质量状况的反应指标之间的差异。第二主成分在职工人均工资、人均居住面积、批发零售贸易商品销售总额、旅游外汇收入的系数上为负，其他为正，而且人均拥有道路面积和人居绿地公共面积的系数的绝对值比较大，说明第二主成分代表了我国居民生活质量状况针对人均公共设施需求（人均拥有道路面积和人居绿地公共面积）和其他居民生活质量状况的反应指标之间的差异。主成分得分表分析，主成分得分表如下：北京2.27203-1.59501天津1.892352-0.862534河北0.0436240.949832山西-0.39209-0.59951内蒙古-1.024830.27295辽宁0.231639-0.009