因子分析
- 格式:doc
- 大小:512.00 KB
- 文档页数:10
现欲对我国18个省市自治区2008年经济发展基本情况的八项指标进行因子分析。
具体采用的指标:GDP (单位:亿元)x1,居民消费水平(单位:元)x2、固定资产投资(单位:亿元)x3、职工平均工资(单位:元)x4、货物周转量(单位:亿吨公里)x5、居民消费价格指数(单位:%)x6、商品零售价指数(单位:%)x7、工业总产值(单位:亿元)x8。
部分数据见表如下:2008年18个省市自治区八项经济指标的部分数据序号 地区 x1x2x3x4x5x6x7x81 北京 10488.03 20346.00 3814.70 56328.00 758.90 105.10 104.40 10413.092 天津 6354.3814000.00 3389.80 41748.02703.40 105.40 105.10 12503.253 河北 16188.61 6570.00 8866.60 24756.00 5925.50 106.20 106.70 23030.734 山西 6938.73 6187.00 3531.20 25828.00 2562.00 107.20 107.20 10023.875 内蒙 7761.80 8108.00 5475.40 26114.00 3658.70 105.70 104.70 8740.186 辽宁 13461.57 9625.00 10019.10 27729.00 7033.90 104.60 105.30 24769.097 吉林6424.06 7591.00 5038.90 23486.01157.80 105.10 106.20 8406.858 黑龙江 8310.00 7039.00 3656.00 23046.00 1690.90 105.60 105.80 7624.549上海13698.15 27343.00 4823.10 56565.00 16029.80 105.80 105.30 25120.9210 江苏 30312.61 11013.00 15300.60 31667.00 4300.90 105.40 104.90 67798.6811 浙江 21486.92 13893.00 9323.00 34146.00 4974.90 105.00 106.30 40832.1012 安徽 8874.17 6377.00 6747.00 26363.00 5843.20 106.20 106.30 11162.1613 福建 10823.11 10361.00 5207.70 25702.00 2396.20 104.60 105.70 15212.8114 江西 6480.33 5753.00 4745.40 21000.00 2285.50 106.00 106.10 8499.5815 山东 31072.06 9573.00 15435.90 26404.90 10107.80 105.30 104.90 62958.5316 河南 18407.78 5877.00 10490.60 24816.00 5165.10 107.00 107.50 26028.4117 湖北 11330.38 7406.00 5647.00 22739.00 2526.40 106.30 106.30 13454.9418 湖南11156.64 7145.00 5534.00 24870.02349.80 106.00 105.60 11553.31具体实验步骤:1、按照顺序:分析——降维——因子分析进入到因子分析的对话框中,将八个变量选入到变量框中。
(2)、单击“描述”按钮,弹出“因子分析:描述统计”对话框,选中“KMO和Bartlett的球形度检验”即可。
如下是KMO检验和Bartlett球体检验的结果。
其中KMO的值为0.440,Bartlett球体检验的p值为0.000<0.05,故认为适合进行因子分析。
KMO 和Bartlett 的检验取样足够度的Kaiser-Meyer-Olkin 度量。
.440Bartlett 的球形度检验近似卡方152.055df 28Sig. .000(3)、点击“抽取”按钮,选中“碎石图”。
此对话框中,表明萃取公因子的方法为主成分分析法,选取公因子的准则为特征值大于1,并显示未经旋转的因子载荷矩阵,运行后的结果如下:公因子方差初始提取GDPX1 1.000 .958居民消费水平x2 1.000 .970固定资产投资x3 1.000 .972职工平均工资x4 1.000 .909货物周转量x5 1.000 .765居民消费价格指数x6 1.000 .839商品零售价格指数x7 1.000 .822工业总产值x8 1.000 .962提取方法:主成份分析。
结论:此表最后一列给出了提取三个公因子后变量的共同度,除了“货物周转量”对应的共同度为0.765外,其他的值都在0.8以上,说明每个变量被提取的三个公因子说明的程度都比较高,损失的信息较少。
解释的总方差成份初始特征值提取平方和载入合计方差的% 累积% 合计方差的% 累积%1 3.432 42.903 42.903 3.432 42.903 42.9032 2.508 31.345 74.247 2.508 31.345 74.2473 1.257 15.718 89.965 1.257 15.718 89.9654 .418 5.227 95.1935 .299 3.736 98.9296 .056 .699 99.6287 .022 .276 99.9048 .008 .096 100.000提取方法:主成份分析。
结论:由此表知,SPSS提取了三个公因子(其特征值都大于1 )前三个因子的累积贡献率达到89.965%,提取三个公因子基本已经足够。
结论:此图也说明提取三个公因子比较合适,到第四个公因子时,特征值已开始趋于平稳。
成份矩阵a成份1 2 3GDPX1 .849.488 .007居民消费水平x2 .546 -.771.278固定资产投资x3 .742.643 -.087职工平均工资x4 .451 -.790.285货物周转量x5 .647-.064 .585居民消费价格指数x6 -.443 .383 .705商品零售价格指数x7 -.522 .546.501工业总产值x8 .879.434 -.031提取方法:主成份。
a. 已提取了3 个成份。
结论:此表为因子载荷矩阵,它是因子命名的主要依据。
从此表可以看出,第一个公因子和第二个公因子在商品零售价格变量上的因子载荷相差不大,第一个公因子个第三个公因子在货物周转量变量上的载荷相差也不大,所以无法明确解释各公因子的含义,也就无法进行因子命名,因此,必须进行因子旋转。
(4)、为了使公因子的含义更容易解释。
,点击“旋转”,并选择“最大方差法”和“旋转解”、“载荷图”得到如下结果:解释的总方差成份初始特征值提取平方和载入旋转平方和载入合计方差的% 累积% 合计方差的% 累积% 合计方差的% 累积%1 3.432 42.903 42.903 3.432 42.903 42.903 3.102 38.778 38.7782 2.508 31.345 74.247 2.508 31.345 74.247 2.309 28.868 67.6463 1.257 15.718 89.965 1.257 15.718 89.965 1.786 22.319 89.9654 .418 5.227 95.1935 .299 3.736 98.9296 .056 .699 99.6287 .022 .276 99.9048 .008 .096 100.000提取方法:主成份分析。
结论:此表与前面的表“解释的总方差”类似,只不过多出旋转后的特征根,方差贡献率以及累计方差贡献率,对比旋转前后的“方差贡献率”,可以发现:前三个公因子的累计贡献率相同(都是89.965%),但每个公因子的特征根发生了变化,以及它们的贡献率也发生了变化。
也就是说,因子旋转相当于在确定的公因子数目下,将相同的累积贡献率在这些公因子上重新分配。
旋转成份矩阵a成份1 2 3GDPX1 .972.084 -.079居民消费水平x2 -.014 .934-.312固定资产投资x3 .977-.124 -.049职工平均工资x4 -.103 .905-.281货物周转量x5 .496 .686.220居民消费价格指数x6 -.110 -.064 .907商品零售价格指数x7 -.081 -.327 .842工业总产值x8 .963.113 -.145提取方法:主成份。
旋转法:具有Kaiser 标准化的正交旋转法。
a. 旋转在5 次迭代后收敛。
成份矩阵a成份1 2 3GDPX1 .849.488 .007居民消费水平x2 .546 -.771.278固定资产投资x3 .742.643 -.087职工平均工资x4 .451 -.790.285货物周转量x5 .647-.064 .585居民消费价格指数x6 -.443 .383 .705商品零售价格指数x7 -.522 .546.501工业总产值x8 .879.434 -.031提取方法:主成份。
a. 已提取了3 个成份。
结论:此表为旋转后的因子载荷矩阵,将其与前面得出的表进行对比,可以发现,个因子在八个变量上的载荷趋于两极分化。
结论:此表为因子旋转的变换矩阵。
结论:此图为旋转后的因子散点图,很清晰的描绘了3个公因子与8个原始变量之间的关系及变量间的关系。
可以看出:工业总产值、GDP、固定资产投资三个变量距离较近,且在第一个公因子轴上的投影坐标较大;居民消费水平、职工平均工资、货物周转量很近,且在第二个公因子轴上的投影坐标很大;但居民消费价格指数、商品零售价格指数很近,但到底是在第三公因子轴还是在第二公因子轴上的投影坐标更大一些,这个比较难以确定。
故,还需要借助于因子载荷矩阵进行因子命名。
(5)、由于本案例涉及变量较多,若要更好更准确进行因子命名,直接通过成分矩阵那个表还不够清晰,所以可以考虑将因子载荷进行排序。
因此,单击“选项”按钮,并选中“按大小排序”运行结果如下:旋转成份矩阵a成份1 2 3固定资产投资x3 .977 -.124 -.049GDPX1 .972 .084 -.079工业总产值x8 .963 .113 -.145居民消费水平x2 -.014 .934 -.312职工平均工资x4 -.103 .905 -.281货物周转量x5 .496 .686 .220居民消费价格指数x6 -.110 -.064 .907商品零售价格指数x7 -.081 -.327 .842提取方法:主成份。