用SPSS进行主成分分析

1、一、引言主成分分析和因子分析在社会经济统计综合评价中是两个常被使用的统计分析方法。现在SPSSSAS等统计软件使用越来越普遍，但SPSS并未像SAS一样，将主成分分析与因子分析作为两个独立的方法并列处理注：主成分分析与因子分析二者是又有着区别与联系，最主要的不同在于它们的数学模型的构建上，具体区别请见参考文献2,而是根据二者之间的关系有机地将主成分分析嵌入到因子分析之中，这样虽然简化了分析程序，却为主成分分析的计算带来不便。且国内许多SPSS教程并没有详细讲解如果应用SPSS进行主成分分析，如何使用SPSS对主成分分析进行计算呢？为使读者能够正确使用SPSS软件进行主成分分析，本文将通过一个实

2、例来详细介绍如何用SPSS做主成分分析。接下来先简单介绍主成分分析原理与模型，以便读者对主成分分析有个大致的了解。二、主成分分析原理和模型1（一）主分成分析原理主成分分析是设法将原来众多具有一定相关性（比如P个指标），重新组合成一组新的互相无关的综合指标来代替原来的指标。通常数学上的处理就是将原来P个指标作线性组合，作为新的综合指标。最经典的做法就是用Fl（选取的第一个线性组合，即第一个综合指标）的方差来表达，即Var（F1）越大，表示Fi包含的信息越多。因此在所有的线性组合中选取的Fi应该是方差最大的，故称Fi为第一主成分。如果第一主成分不足以代表原来P个指标的信息，再考虑选取F2即选第二个

3、线性组合，为了有效地反映原来信息，Fi已有的信息就不需要再出现再F2中，用数学语言表达就是要求Cov（Fi,F2）=0,则称F2为第二主成分，依此类推可以构造出第三、第四，第P个主成分。（二）主成分分析数学模型F2=ai2ZX+a22Z%,+a4Fp=aimZXi+a2mZX2+amZXp其中aii,a2i,pi（jFi,沏X的协方差阵2的特征值多对应的特征向量，ZXi,ZX2,ZXp原始变量经过标准化处理的值，因为在实际应用中，往往存在指标的量纲不同，所以在计算之前须先消除量纲的影响，而将原始数据标准化，本文所采用的数据就存在量纲影响注：本文指的数据标准化是指Z标准化。A=（aij）pxm=

4、（ai,a2,知,），Rai=入iai,R为相关系数矩阵，入、ai是相应的特征值和单位特征向量，入n2>->p>0。进行主成分分析主要步骤如下：i.指标数据标准化（SPSS软件自动执行）；2.指标之间的相关性判定；3. 确定主成分个数m;4. 主成分Fi表达式；5. 主成分Fi命名；主成分与综合主成分（评彳）值。三、对沿海10个省市经济综合指标进行主成分分析（一）指标选取原则本文所选取的数据来自中国统计年鉴2003»中2002年的统计数据，在沿海10省市经济状况主要指标体系中选取了10个指标：X1GDPX2人士匀GDPX3农业增加值X4工业增加值X5第二产业增加值X

5、6固定资产投资X7基本建设投资X8国内生产总值占全国比重（）X9海关出口总额X10地方财政收入图表1沿海10个省市经济数据地区GDP人均GDP农业增加值工业增加值第二产业增加值固定资产投资基本建设投资社会消费品零售总额海关出口总额地方财政收入辽宁5458.21300014883.31376.22258.41315.95292258.4123.7399.7山东105501164313903502.538512288.71070.73181.9211.1610.2河北6076.69047950.21406.72092.61161.6597.11968.345.9302.3天津2022.622068

6、83.9822.8960703.7361.9941.4115.7171.8江苏10636143971122.63536.33967.223201141.33215.8384.7643.7上海5408.84062786.22196.22755.81970.2779.32035.2320.5709浙江7670165706802356.530652296.61180.62877.5294.2566.9福建4682135106631047.11859964.5397.91663.3173.7272.911770150301023.94224.64793.63022.91275.55013.61843.

7、71202广西2437.25062591.4367995.7542.2352.71025.515.1186.7（二）主成分分析在SPSS中的具体操作步骤运用SPSS统计分析软件Factor过程对沿海10个省市经济综合指标进行主成分分析。具体操作步骤如下：1. AnalyzeaDataReductionFFactorAnB竭isFactorAnalysis对话框2. 把XiXio选入Variables框3. Descriptives:CorrelationMatrix框组中选中Coefficients,然后点击Continue,返回FactorAnalysis对话框4. 点击“OK”图表2Fac

8、torAnalyze对话中1与Descriptives子对话框.r峪田£的3匕”广日毕血£4dQdniWWWr婕TMQE厂y*0vdBfldfeiYI-汕ec即SPSS在调用FactorAnalyze过程进行分析时，SPSS会自动对原始数据进行标准化处理，所以在得到计算结果后指的变量都是指经过标准化处理后的变量，但SPSS不会直接给出标准化后的数据，如需要得到标准化数据，则需调用Descriptives过程进行计算。图表3相关系数矩阵人为砌耸Pw如惊工业事m也国引雷产枚寝S*4ttttl社曾清中国手笆自如雄吴*口由一ttA0WI.W09MJO幻W$及923WIOT且M人电面

9、T31蝴一，1U3侬前帕1稣我生曲灿玲-0)3-1?11.000-.1J3通知，.媾IW01J.12J必H业脚匐曲M7LB门骁1.000那购9i9.935.TO5脐&黑二产业用国总919m-J030硼1A00晒Wt例71491!固定,卢根费靖JH确W3帆10009fl9部MX外4事上健立恨傅g期H116卡器M0邮1ItM蒯d制合社金由受星零告百中HlW01?的咖】碗期9»岛养出口名澳J«WJtl-innw.714.7176N颜i.m雨心开时就收人釉5m叫wsVI3934州83100图表4方差分解主成分提取分析表TotalVarianceEjlainedMithl音中

10、而屿ExfndSnSumsofSgudLoadingsTsui%afVirimwCumulate惴T而惴ofVmmceCm皿囱如$L7.220n.su7.210理加572,2052L23512.34664.5511.2351234664,5513的8砂53.31S4,5475耶598,7865.0S5部4沔川6.021.211998501012,11999.9K8.002.01899.9889D01,013100.(W10.000.000100.000EstrajctiatLMe由第1:Pi讪甲虹Cob平on修正痴型?it.图表5初始因子载荷矩阵ComponentMatrix*Componen

11、t12GDP5491S5人均GDP112-824农业增加值-.109.577工业增加值,978-.005第三产业增加值.070固定资产投资,983-J68基本建设投资347-.024社会消费品零售总额,977.176海关出口总额.800-051地方财政收入.954-.126ExtrActianMethod;PmuripAlComponentAnalysis.32匚qequrn由extrafted.从图表3可知GD*工业增加值，第三产业增加值、固定资产投资、基本建设投资、社会消费品零售总额、地方财政收入这几个指标存在着极其显著的关系，与海关出口总额存在着显著关系。可见许多变量之间直接的相关性比较

12、强，证明他们存在信息上的重叠。主成分个数提取原则为主成分对应的特征值大于1的前m个主成分。注：特征值在某种程度上可以被看成是表示主成分影响力度大小的指标，如果特征值小于1,说明该主成分的解释力度还不如直接引入一个原变量的平均解释力度大，因此一般可以用特征值大于1作为纳入标准。通过图表4（方差分解主成分提取分析）可知，提取2个主成分，即m=2,从图表5（初始因子载荷矩阵）可知GDR工业增加值、第三产业增加值、固定资产投资、基本建设投资、社会消费品零售总额、海关出口总额、地方财政收入在第一主成分上有较高载荷，说明第一主成分基本反映了这些指标的信息；人均GD环口农业增加值指标在第二主成分上有较高载荷

13、，说明第二主成分基本反映了人均GD济口农业增加值两个指标的信息。所以提取两个主成分是可以基本反映全部指标的信息，所以决定用两个新变量来代替原来的十个变量。但这两个新变量的表达还不能从输出窗口中直接得到，因为"ComponentMatrix”是指初始因子载荷矩阵，每一个载荷量表示主成分与对应变量的相关系数。用图表5（主成分载荷矩阵）中的数据除以主成分相对应的特征值开平方根便得到两个主成分中每个指标所对应的系数2。将初始因子载荷矩阵中的两列数据输入（可用复制粘贴的方法）到数据编辑窗口（为变量B1、B2），然后利用aTransformaComputeVariable:在ComputeVar

14、iable对话框中输入“A1=B1/SQR（7.22）”注：第二主成分SQ所的括号中填1.235,即可得到特征向量A1（见图表6）。同理，可得到特征向量A2。将得到的特征向量与标准化后的数据相乘，然后就可以得出主成分表达式注：因本例只是为了说明如何在SPSS进行主成分分析，故在此不对提取的主成分进行命名，有兴趣的读者可自行命名:Fi=0.353ZXi+0.042ZX2-0.041ZX3+0.364ZX4+0.367ZX5+0.366ZX6+0.352ZX7+0.364ZX8+0.298ZX9+0.355ZXi0F2=0.175ZXi-0.741ZX2+0.609ZX3-0.004ZX4+0.0

15、63ZX5-0.061ZX6-0.022ZX7+0.158ZX8-0.046ZX9-0.115ZX10图表6ComputeVariable对话框前文提到SPS8自动对数据进行标准化，但不会直接给出，需要我们自己另外算，我们可以通过AnalyzeaDescriptiveStatisticsaDescriptives对话框来实现：弹出Descriptives对话框后,把XX10选入Variables框，在Savestandardizedvaluesasvariables前的方框打上钩，点击“OK,经标准化的数据会自动填入数据窗口中，并以Z开头命名。图表7Descriptives对话框以每个主成分所

16、对应的特征值占所提取主成分总的特征值之和的比例作为权重计算主成分综合模型：即可得到主成分综合模型：F=0.327ZXi-0.072ZX2+0.054ZX3+0.310ZX4+0.323ZX5+0.304ZX6+0.297ZX7+0.334ZX8+0.248ZX9+0.286ZXi0根据主成分综合模型即可计算综合主成分值，并对其按综合主成分值进行排序，即可对各地区进行综合评价比较，结果见图表8。图表8综合主成分值城巾第一主成分F1排名第二主成分F2排名综合主成分F排名5.2310.1164.481江苏2.2520.2351.962山东1.9630.5021.753浙江1.164-0.1980.964上海0.305-2.3610-0.095辽宁-1.2461.961-0.786河北-1.3570.414-1.107福建-1.978-0.077-1.708天津-3.049-1.019-2.749广西-3.29100.413-2.7510对得出的综合主成分(评价)值，我们可用实际结果、经验与原始数据做聚类分析进行检验，对有争议的结果，可用原始数据做判别分析解决争议，具体评价与检验本文