应用统计学因子分析与主成分分析案例解析SPSS操作分析

1、因子分析与主成分分析学号：2010209583 姓名：钱志雷 专业：信息管理与信息系统 班级：四班一、 问题概述现希望对30个省市自治区经济发展基本情况的八项指标进行分析。具体采用的指标只有：gdp、居民消费水平、固定资产投资、职工平均工资、货物周转量、居民消费价格指数、商品零售价格指数、工业总产值。这是一个综合分析问题，八项指标较多，用主成分分析法进行综合。二、 数据由于样本数比较多，这里不再给出，可参见factor1.sav文件三、数据处理与分析1. 因子分析打开数据后，在spss中进行因子分析的步骤如下：选择“分析-降维-因子分析”，在弹出的对话框里（1）描述-系数、kmo与bartle

2、tt的球形度检验（2）抽取-碎石图、未旋转的因子解（3）旋转-最大方差法、旋转解、载荷图（4）得分-保存为变量、显示因子得分系数矩阵（5）选项-按大小排序点击确定得到如下各图： 图3-1相关矩阵gdp居民消费水平固定资产投资职工平均工资货物周转量居民消费价格指数商品价格指数工业总产值相关gdp1.000.267.951.187.617-.273-.264.874居民消费水平.2671.000.426.716-.151-.235-.593.363固定资产投资.951.4261.000.396.431-.280-.359.792职工平均工资.187.716.3961.000-.357-.145-.

3、543.099货物周转量.617-.151.431-.3571.000-.253.022.659居民消费价格指数-.273-.235-.280-.145-.2531.000.763-.125商品价格指数-.264-.593-.359-.543.022.7631.000-.192工业总产值.874.363.792.099.659-.125-.1921.000图3-2kmo 和 bartlett 的检验取样足够度的 kaiser-meyer-olkin 度量。.620bartlett 的球形度检验近似卡方231.285df28sig.000图3-3公因子方差初始提取gdp1.000.945居民消费

4、水平1.000.799固定资产投资1.000.902职工平均工资1.000.873货物周转量1.000.857居民消费价格指数1.000.957商品价格指数1.000.928工业总产值1.000.904提取方法：主成份分析。 图3-4解释的总方差成份初始特征值提取平方和载入旋转平方和载入合计方差的 %累积 %合计方差的 %累积 %合计方差的 %累积 %13.75446.92446.9243.75446.92446.9243.20740.09240.09222.20327.53274.4562.20327.53274.4562.21727.70867.80031.20815.09689.5511

5、.20815.09689.5511.74021.75289.5514.4035.04294.5935.2142.67397.2666.1381.72298.9887.066.82999.8178.015.183100.000提取方法：主成份分析。图3-5图3-6成份矩阵a成份123固定资产投资.911.163.213gdp.884.385.120工业总产值.822.429.210商品价格指数-.621.596.433居民消费水平.606-.596.277货物周转量.486.737-.279职工平均工资.465-.725.362居民消费价格指数-.510.257.794提取方法 :主成份。a.

6、已提取了 3 个成份。图3-7旋转成份矩阵a成份123gdp.955.124-.131工业总产值.944.109-.014固定资产投资.872.351-.137货物周转量.751-.507-.192职工平均工资.048.925-.121居民消费水平.219.841-.209居民消费价格指数-.135-.013.969商品价格指数-.104-.496.819提取方法 :主成份。 旋转法 :具有 kaiser 标准化的正交旋转法。a. 旋转在 5 次迭代后收敛。图3-8成份转换矩阵成份1231.817.407-.4082.548-.769.3313.179.494.851提取方法 :主成份。 旋转

7、法 :具有 kaiser 标准化的正交旋转法。 图3-9图3-10成份得分系数矩阵成份123gdp.306.011.047居民消费水平.025.387.040固定资产投资.270.129.075职工平均工资-.025.451.096货物周转量.248-.319-.139居民消费价格指数.070.180.653商品价格指数.077-.098.462工业总产值.317.026.123提取方法 :主成份。 旋转法 :具有 kaiser 标准化的正交旋转法。 构成得分。图3-11成份得分协方差矩阵成份12311.000.000.0002.0001.000.0003.000.0001.000提取方法 :

8、主成份。 旋转法 :具有 kaiser 标准化的正交旋转法。 构成得分。 （2）因子模型中各统计量的意义 a）因子载荷aij：因子载荷aij为第i个变量在第j个因子上的载荷，实际上就是xi与fj的相关系数，表示变量xi依赖因子fj的程度，反应了第i个变量xi对于第j个因子fj的重要性。 b）变量xi的变量共同度：k个公因子对第i个变量方差的贡献，也称为公因子方差比，记为hi2，公式为：hi2=i=1paij2 （j=1,2,.,k）表示全部公因子对变量xi的总方差所做出的贡献，也即是变量xi的信息能够被k个公因子所描述的程度。 c）公因子fj的方差贡献率：在因子载荷矩阵a中，各列元素aj的平方

9、和记为gi2，表示第j个公因子fj对于x所提供方差的总和，它是衡量公因子相对重要性的指标。方差贡献率越大，表明公因子对x的贡献越大。（3）基本分析结果a）kmo和球形bartlett检验用于因子分析的适用性检验。kmo检验变量间的偏相关是否较小，bartlett球形检验是判断相关矩阵是否是单位阵，参见图3-2。由bartlett检验可以看出，应拒绝个变量独立的假设，即变量间具有较强的相关性，但是kmo的统计量为0.620，小于0.7，说明个变量间信息的重叠程度可能不是特别的高，有可能做出的因子分析模型不是很完善，但还是值得尝试的。b）变量共同度communalities是表示各变量中所含原始信

10、息能被提取的公因子所表示的程度，由图3-3所示的变量共同度可知：几乎所有变量的共同度都在80%以上，因此提取出的这几个公因子对各变量的解释能力是较强的。c）碎石图用于显示各因子的重要程度，横轴为因子序号，纵轴表示特征根大小，从中可以非常直观的了解到哪些是最主要的因子，参见图3-5。本例中可见前三个因子的散点位于陡坡之上，而后五个因子散点成了平台，且特征根均小于1，因此至多考虑前三个公因子即可。d）图3-4给出的是各成分的方差贡献率和累计贡献率，以及进行因子旋转后的方差贡献率和累计贡献率，前者将在主成分分析中进行说明。e）图3-6为因子载荷矩阵，在前面已经直接按列的方向将其解释为个成分的系数，实

11、际上严格讲因子载荷矩阵应该是各因子在各变量上的载荷，即是各因子对各变量的影响度。表示如下：zx1=0.884f1+0.385f2+0.120f3+1zx2=0.606f1-0.596f2-0.277f3+2.zx8=0.822f1+0.429f2-0.210f3+8在表达式中各变量已经不是原始变量，而是标准化变量。i表示特殊因子，是除了这3个公因子之外影响该变量的其他因素。原来设计了8个指标来表示经济发展水平，但是经过因子分析后，只需要三个因子即可描述影响地区经济发展状况。f）为了使因子载荷矩阵中系数更加显著，可以对初始因子载荷矩阵进行转换，使因子和原始变量间的关系进行重新分配，相关系数向0-

12、1分化，从而更加容易解释。图3-9是进行因子旋转的空间示意图，值得注意的是旋转前后各变量散点的相对位置保持不变，即旋转并不改变因子分析的整体结果，只是影响各因子在各变量上的载荷分布，并影响各因子的贡献率。本例中采用的是方差最大正交旋转法进行因子旋转，输出的结果参见图3-4.，由图可知，只有前三个特征根大于1，因此spss只提取了前三个公因子。在旋转后三个公因子的方差累计贡献率均发生了变化，但仍然会保持从大到小的顺序，而且前三个因子的方差贡献率仍为89.55%，和旋转前完全相同，因此选前三个因子已足够描述经济发展的水平。g）进行方差最大旋转后，旋转后的因子载荷矩阵如图3-7所示，由图可以看出，第

13、一公因子在x1、x3、x5、x8有较大的载荷，主要从gdp、固定资产投资、货物周转量和工业总产值反映经济发展状况，可以命名为总量因子。第二公因子在x2、x4上有较大载荷，从居民消费水平和职工平均工资方面反映经济发展水平，因此命名为消费因子。第三公因子在x6和x7上有较大载荷，表现为居民消费价格指数和水平价格指数方面，因此命名为价格因子。与未旋转前相比较，旋转后各公因子的意义显然更加明确合理。h）因子得分：前面得到了因子结构表达式， 可以将各变量表示为公因子的线性形式，但是更多的时候需要将公因子表达为各变量的线性形式。公因子的得分系数函数不能通过矩阵变换的方法由因子载荷阵得到，只能采用估计的方法

14、求得，本例采用的是回归法。因子得分系数矩阵如图3-10所示，据此可以直接写出各公因子的得分表达式：f1=0.306zx1+0.025zx2+0.270zx3-0.025zx4+0.248zx5+0.070zx6+0.077zx7+0.317zx8f2=0.011zx1+0.387zx2+0.129zx3+0.451zx4-0.319zx5+0.180zx6-0.098zx7+0.026zx8f3=0.047zx1+0.040zx2+0.075zx3+0.096zx4-0.139zx5+0.653zx6+0.462zx7+0.123zx8spss已经给出三个公因子的得分，保存在fac_1fac_

15、3中，按各因子对应的方差贡献率为权数计算如下综合统计量： f=11+2+3f1+21+2+3f2+31+2+3f3 =0.730f1+0.141f2+0.129f3在spss中用程序计算综合因子得分：comp score=0.73* fac1_1+0.141* fac2_1+0.129* fac3_12. 主成分分析a)由图3-1（各变量相关系数矩阵）可以看出，许多变量之间直接的相关性比较强，的确存在信息上的重叠。b)由图3-4（具体不再阐述）可知，只有前三个特征根大于1，因此spss只提取了前三个主成分，前三个主成分的方差贡献率达到89.55%，因此选前三个主成分已足够描述经济发展的水平。c

16、)图3-6输出为主成分系数矩阵，从而得到各主成分的表达式，在表达式中各变量已经不是原始变量，而是标准化变量。f1=0.884zx1+0.606zx2+0.911zx3+0.465zx4+0.486zx5-0.51zx6-0.621zx7+0.822zx8f2=0.385zx1-0.596zx2+0.163zx3-0.725zx4+0.737zx5+0.257zx6-0.596zx7+0.429zx8f3=0.120zx1+0.277zx2+0.213zx3+0.362zx4-0.279zx5+0.794x6-0.433zx7+0.210zx8因为各自变量已经过标准化，因此以上三个主成分的均数均为0。可以证明，各主成分的方差应当为前述特征