会城及计划单列城综合经济实力评价指标统计分析

1、基于城市综合经济实力评价指标的统计分析张三摘要本文采用中国2010年省会城市和计划单列市主要经济指标统计数据，首先利用相关分析对地区生产总值各产业间进行了简单相关分析，以及控制了年末总人口对客运和货运量等进行了偏相关分析；接着利用各种经济指标对地区生产总值进行多元线性回归，利用逐步回归方法获知社会商品零售总额、剧场影剧院、普通高等学校在校学生数、固定资产投资总额、货物进出口总额、三废综合利用产品产值、客运量等和地区生产总值有显著关系；最后利用因子分析对各地区进行因子分析，通过计算综合得分，获得各地区综合排名，其中北京、上海一枝独秀，重庆、深圳、天津与北京、上海构成前五强。关键词地区生产总值；相

2、关分析；回归分析；因子分析1研究背景及目的2009年10月17日，第六届中国城市论坛北京峰会在朝阳规划艺术馆召开。这次峰会不仅吸引了城市发展领域内几百位专家的专注和参与，跟让来自全国各个城市的会议代表们受益匪浅。会议指出，“十二五”期间既是全球经济复苏的关键时期，也是我国加快城市化进程的关键时期。以前我国采取的城市外延式扩张战略导致城市发展中出现了资源浪费、环境污染、不注重保护城市历史文脉和特点等各种各样的问题。城市发展将呈现出五个新变化：一是城市发展开始从外延式扩张向内涵式发展转变；二是城市软实力成为城市发展的核心竞争力；三是城乡统筹和城乡一体化成为城市发展的新格局；四是综合配套改革实验区的

3、示范意义进一步凸显；五是城市群对城市建设与发展的作用日益增强。在这种大背景下对我国各城市的综合经济实力进行研究，不论是对于促进我国城市本身又好又快地发展，还是对于充分发挥城市在社会经济生活中所起的主导作用，都有着极为重要的意义。本研究的目的如下：通过对描述我国各城市综合经济实力的各种指标进行各种分析，一方面找出用来衡量我国城市综合经济实力的各个指标之间的内在联系，另一方面找出各城市综合经济实力的差异，为城市的发展提供政策建议1。2研究方法对城市综合经济实力的概念，中国城市经济发展研究中心指出：城市综合经济实力是指城市所拥有的全部实力、潜力及其在国内外经济社会中的地位和影响力。据此概念我们可以看

4、出，评价城市综合经济实力应该包括人口、地区生产总值、拥有的交通运输以及通讯能力、地方财政预算内收支、固定资产投资总额、城乡居民工资水平及储蓄水平、环境污染治理投资总额、商贸市场水平、人才状况及社会医疗保障水平等方面。所以我们采用的数据指标有：年末总人口、地区生产总值、第一产业、第二产业、第三产业、客运量、货运量、地方财政预算内收入、地方财政预算内支出、固定资产投资总额、城乡居民储蓄年末余额、在岗职工平均工资、年末邮政局数、年末固定电话用户数、社会商品零售总额、货物进出口总额、年末实有公共（汽）电车运营车辆数、剧场、影剧院数、普通高等学校在校学生数、医院、卫生院个数、执业（助理）医师、三废综合利

5、用产品总值等22个指标。本研究采用的数据是中国2010年省会城市和计划单列市主要经济指标统计（包括直辖市）»,数据摘自中国统计年鉴201111-3。采用数据分析方法主要有回归分析2、相关分析、因子分析等。基本思路是：首先使用回归分析、相关分析等分析方法研究构成城市综合经济实力的各个变量之间的关系；然后使用因子分析对构成城市综合经济实力的各个变量提取公因子；最后使用一些简单的SPSS数据处理技巧依据提出的公因子对各城市进行分类及排序。3实证分析政预算内收入”、“地方财政预算内支出”由于从国家统计局网站下载的数据为EXCEL格式，可以将数据导入成SPSS数据,我们设置了共23个变量，分别

6、是“城市名称”、“年末总人口”、“地区生产总值”、“第一产业”、“第二产业”、“第三产业”、“客运量”、“货运量”、“地方财、“固定资产投资总额”、“城乡居民储蓄年末余额”、“在岗职工平均工资”、“年末邮政局数”、“年末固定电话用户数”、“社会商品零售总额”、“货物进出口总额”、“年末实有公共（汽）电车运营车辆数”、“剧场、影剧院数”、“普通高等学校在校学生数”、“医院、卫生院数”、“执业（助理）医师”、“三废综合利用产品总值”等。样本是中国2010年省会城市和计划单列市主要经济指标统计的相关数据。数据参见附表1。3.1 相关分析对于相关分析，主要有以下几个部分：第一，对“地区生产总值”的三个

7、组成部分：“第一产业”、“第二产业”和“第三产业”进行简单相关分析；第二，在控制“地区生产总值”的前提下，对“第一产业”、“第二产业”、“第三产业”进行偏相关分析；第三，在控制“年末总人口”的前提下，对“客运量”和“货运量”进行偏相关分析；第四，在控制“年末总人口”的前提下，对“地方政府预算内收入”和“地方政府预算内支出”进行偏相关分析；第五，对“年末总人口”、“地区生产总值”、“三废综合利用产品总值”这三个变量进行简单相关分析。3.1.1 地区生产总值的简单相关分析通过使用主菜单分析下的相关工具条，可以获得分析结果，见表1。表1产业相关分析结果相关性第一产业第二产业第三产业Pearson相关

8、性1.295.063第一产业显著性（双侧）.085.718N353535Pearson相关性.2951*.820第二产业显著性（双侧）.085.000N353535Pearson相关性.063*.8201第三产业显著性（双侧）.718.000N353535*.在.01水平（双侧）上显著相关。易见只有“第二产业”和“第三产业”之间具有很强的相关关系并且在0.01的显著性水平上显著。这是因为“第二产业”和“第三产业”本身就具有很大的相关性。3.1.2 第一产业、第二产业和第三产业的偏相关分析表2偏相关分析结果相关性控制变量第一产业第二产业第三产业相关性1.000.321-.517第一产业显著性（双

9、侧）.064.002df03232相关性.3211.000-.977地区生产总值第二产业显著性（双侧）.064.000df32032相关性-.517-.9771.000第三产业显著性（双侧）.002.000.df32320从分析结果可以看出，在控制“地区生产总值”的前提下，构成“地区生产总值”的3个组成部分中，“第二产业”和“第三产业”的相关关系依然显著。但相关系数却变成了负值，这是由于总额是一定。3.1.3 客运量和货运量的偏相关分析表3偏相关分析结果相关性控制变量客运量货运量相关性1.000.127客运量显著性（双侧）.474df032年末总人口相关性.1271.000货运量显著性（双侧）

10、.474.df320可见在控制“年末总人口”的前提下，“客运量”和“货运量”之间的相关关系不显著3.1.4 地方财政预算内收入和地方财政预算内支出的偏相关分析分析结果如表4所示:表4偏相关分析结果相关性控制变量地方财政预算内收入地方财政预算内支出相关性1.000.994地方财政预算内收入显著性（双侧）.000df032年末总人口相关性.9941.000地方财政预算内支出显著性（双侧）.000.df320可见，在控制“年末总人口”的前提下，“地方财政预算内收入”和“地方财政预算内支出”之间的相关系数很大，而且相关关系非常显著。3.1.5 年末总人口、地区生产总值、三废综合利用产品总值的相关分析表

11、5简单相关分析结果相关性年末总人口地区生产总值三废综合利用值Pearson相关性1*.524.128年末总人口显著性（双侧）.001.479N353533Pearson相关性*.5241.153地区生产总值显著性（双侧）.001.394N353533Pearson相关性.128.1531三废综合利用值显著性（双侧）.479.394N333333*.在.01水平（双侧）上显著相关。分析结果如表5所示，可见，“年末总人口”与“地区生产总值”正相关而且这种相关关系十分显著；“地区生产总值”与“三废综合利用产品总值”之间也存在显著的相关关系；“年末总人口”与“三废综合利用产品总值”之间的相关关系不显著

12、。3.2回归分析对于回归分析，以“地区生产总值”为被解释变量，“年末总人口”、“客运量”、“货运量”、“地方财政预算内收入”、“地方财政预算内支出”、“固定资产投资总额”、“城乡居民储蓄年末余额”、“在岗职工平均工资”、“年末邮政局数”、“年末固定电话用户数”、“社会商品零售总额”、“货物进出口总额”、“年末实有公共（汽）电车营运车辆数”、“剧场、影剧院数”、“普通高等学校在校学生数”、“医院、卫生院数”、“执业（助理）医师”、“三废综合利用产品总值”等为解释变量，进行多元线性回归，采用逐步回归法。得到输出结果参见表6-表8所示。表6模型拟合情况表模型汇总模型RR方调整R方标准估计的误差Dur

13、bin-Watson1.975a.951.949837.32417.998g.995.994288.45561.595a.预测变量：（常量），社会商品零售总额。g.预测变量：（常量），社会商品零售总额，剧场、影剧院，高等学校学生数，国内资产投资总额，货物进出口，客运量，三废综合利用值。h.因变量：地区生产总值a表7方差分析表模型平方和df均方FSig.回归417501904.5831417501904.583595.486.000b1残差21734459.78231701111.606总计439236364.36432回归437156198.448762450885.493750.552.00

14、0h7残差2080165.9172583206.637总计439236364.36432Anovaa.因变量：地区生产总值b.预测变量：（常量），社会商品零售总额。h.预测变量：（常量），社会商品零售总额，剧场、影剧院，高等学校学生数，国内资产投资总额，货物进出口，客运量，三废综合利用值。表8回归系数表模型非标准化系数标准系数tSig.B标准误差试用版(常量)1-107.524238.711-.450.656社会商品零售总额.000.000.97524.403.000(常量)16.834113.567.148.883社会商品零售总额.000.000.7238.873.000剧场、影剧院-14.

15、7683.067-.137-4.815.000高等学校学生数7-.001.000-.058-2.430.023国内资产投资总额5.243E-005.000.2417.273.000货物进出口.000.000.3026.098.000客运量-.007.002-.071-3.952.001三废综合利用值.001.000.0523.666.001a.因变量：地区生产总值综上所述，我们可以得到:(1)最终模型的表达式为：地区生产总值=16.834+0.000*社会商品零售总额-14.768*剧场、影剧院-0.001*高等学校学生数+0.00005*固定资产投资总额+0.000*货物进出口+0.001*

16、三废综合利用产品产值-0.007*客运量；(2)最终模型的拟合优度较好，修正后可决系数近乎1;(3)模型中各自变量系数的显著性p值都小于0.05,回归方程的线性关系显著；(4)分析结论：经过以上多元线性回归分析，可以发现我国城市的地区生产总值与社会商品零售总额、剧场、影剧院、普通高等学校在校学生数、固定资产投资总额、货物进出口总额、三废综合利用产品产值、客运量有显著关系，与其他变量之间的关系并不显著。表9残差统计表残差统计量极小值极大值均值标准偏差N预测值419.71017058.1454505.6193696.096833残差-767.9445547.0401.0000254.961133标

17、准预测值-1.1053.396.0001.00033标准残差-2.6621.896.000.88433a.因变量：地区生产总值表9给出了回归分析的残差统计结果，可以看出预测值及标准化的预测值、残差及残差预测值的最小值、最大值、均值、标准差和样本数。这些数据中无离群值，可以认为模型是健康的。除了分析残差统计外，还可以直接做出标准残差的直方图和正态P-P图来观察其是否服从正态分布。从图1和图2可以看出残差具有正态分布的趋势，因此可以认为回归模型是恰当的。直方图图1标准化残差直方图标腐化残雉的标制P-P图期并晌俎和概率图2标准化残差正态P-P图上一部分我们利用线性回归来对各城市地区生产总值的影响因素

18、进行相关和回归分析，获知各种影响指标，以及影响程度和方向，但是由于各种变量的多重共线性问题，不太利于具体经济意义解释，因此我们可以利用因子分析对构成城市综合经济实力的各个变量提取公因子，对数据进行降维处理。首先利用主菜单分析中的降维工具条进行因子分析。在对话框中，依次将“年末总人口”、“地区生产总值”、“客运量”、“货运量”、“地方财政预算内收入”、“地方财政预算内支出”、“固定资产投资总额”、“城乡居民储蓄年末余额”、“在岗职工平均工资”、“年末邮政局数”、“年末固定电话用户数”、“社会商品零售总额”、“货物进出口总额”、“年末实有公共(汽)电车运营车辆数”、“剧场、影剧院数”、“普通高等学

19、校在校学生数”、“医院、卫生院数”、“执业(助理)医师”、“三废综合利用产品总值”添加到“变量”列表中。接着在相应的按钮中选中输出相关系数矩阵和KMO和Bartlett的球形度检验的复选框，并且设置特征根大于1为提取公因子的要求，利用最大方差法进行因子旋转，最后保存因子得分变量，并显示因子得分系数矩阵，输出结果参见表10-表16。(1) KMO和Bartlett的检验，如表10所示,KMO的取值为0.844,表明变量间有较强的相关性，数据很适合做因子分析。Bartlett检验的Sig.值为0.000,说明数据来自正态总体，适合进一步分析。表10KMO和Bartlett的检验结果KMOandBa

20、rtlett'sTestKaiser-Meyer-OlkinMeasureofSamplingAdequacy.Bartlett'sTestofApprox.Chi-SquareSphericitydfSig.8441123.563171.000(2)变量共同度，指的是按照所选标准提取相应数量主成分后，各变量中信息分别被提取的比例。如表11所示，除了“三废综合利用值”的共同度较低之外，大多数变量共同度都在90犯上，所以提取这几个公因子对各变量的解释力还可以。表11变量共同度CommunalitiesInitialExtraction年末总人口1.000.948地区生产总值1.0

21、00.962客运量1.000.697货运量1.000.735地方财政预算内收入1.000.957地方财政预算内支出1.000.951国内资产投资总额1.000.835城乡居民储蓄1.000.969在岗职工平均工资1.000.876年末邮政局数1.000.865年末固定电话用户1.000.970社会商品零售总额1.000.974货物进出口1.000.962年末实有公共汽1.000.859剧场、影剧院1.000.673高等学校学生数1.000.633医院、卫生院1.000.922执业助理医生1.000.844三废综合利用值1.000.341ExtractionMethod:PrincipalCom

22、ponentAnalysis.(3)解释的总方差，由表12中可以看出，“初始特征值”一栏显示只有前三个特征值大于1,所以只选取了前三个公因子；”提取平方和载入”一栏显示第一公因子的方差贡献率是63.754%,前三个公因子的方差总和占所有主成分方差的84.067%,可见选取前三个因子已足够替代原来的变量。表12解释总方差TotalVarianceExplainedComponentInitialEigenvaluesExtractionSumsofSquaredLoadingsRotationSumsofSquaredLoadingsTotal%ofVarianceCumulative%Tota

23、l%ofVarianceCumulative%Total%ofVarianceCumulative%112.11363.75463.75412.11363.75463.7548.96747.19347.19322.65913.99477.7482.65913.99477.7485.54029.15876.35131.2016.31984.0671.2016.31984.0671.4667.71684.0674.9334.91088.9775.7754.07793.0546.5462.87495.9287.2311.21697.1448.137.72297.8679.118.62398.4901

24、0.083.43798.92711.071.37399.30012.052.27399.57313.026.13799.71014.021.10999.81915.015.07799.89516.010.05599.95017.006.03399.98318.002.01199.99419.001.006100.000ExtractionMethod:PrincipalComponentAnalysis.(4)碎石图，是按照特征根大小排列的主成分散点图。如图3所示，我们可以看到有三个成分的特征值超过1。STHO7=81234K47ft91011131415ift171019成分数图3碎石图(5

25、)成分矩阵，如表13可见，在所有变量中打*号的三个变量归为哪一因子不是很明确，因此有必要进行因子旋转。我们利用最大方差法进行因子旋转得到的旋转成分矩阵如表14所示。表13成分矩阵ComponentMatrixComponent123年末固定电话用户.980-.093.026社会商品零售总额.970-.137.123地方财政预算内支出.967-.120-.029地区生产总值.965-.165.056城乡居民储蓄.949-.248.078地方财政预算内收入.944-.254.035执业助理医生.909.129.036年末实有公共汽.810-.387-.231国内资产投资总额.810.409.113

26、在岗职工平均工资.790-.469.177货物进出口.777-.570-.181货运量.762.392.022剧场、影剧院.759-.182.252年末邮政局数.711.475-.366客运量.669.054-.497医院、卫生院*.647.697-.130年末总人口*.676.697-.073三废综合利用值.154.084.557高等学校学生数*.420.438.514表14旋转成分矩阵Component123货物进出口.961.032-.192城乡居民储蓄.920.317.151在岗职工平均工资.918.029.178地方财政预算内收入.918.320.107地区生产总值.884.397.

27、149年末实有公共汽.882.207-.197社会商品零售总额.873.405.221地方财政预算内支出.859.455.077年末固定电话用户.855.470.138剧场、影剧院.729.218.308执业助理医生.668.602.185年末总人口.150.947.169医院、卫生院.124.945.111年末邮政局数.302.866-.156国内资产投资总额.427.748.306货运量.397.731.208客运量.508.543-.379高等学校学生数.097.451.648三废综合利用值.085.012.578aRotatedComponentMatrix这样每个因子就很明确了。第一

28、因子上包括地区生产总值、地方财政预算内收入、地方财政预算内支出、城乡居民储蓄年末余额、在岗职工平均工资、年末固定电话用户数、社会商品零售总额、货物进出口总额、年末实有公共（汽）电车营运车辆数、剧场、影剧院、执业（助理）医师，可命名为存量因子；第二因子包括年末总人口、客运量、货运量、固定资产投资总额、年末邮政局（所）数、医院、卫生院，可命名为流量因子；第三因子包括普通高等学校在校学生数和三废综合利用产品产值，可命名为可持续因子。（6）成分得分系数矩阵表15成分得分系数矩阵ComponentScoreCoefficientMatrixComponent123年末总人口-.106.253.005地区

29、生产总值.101-.016.043客运量.028.151-.389货运量-.033.147.058地方财政预算内收入.119-.039.018地方财政预算内支出.091.015-.022国内资产投资总额-.033.135.133城乡居民储蓄.119-.046.054在岗职工平均工资.157-.140.114年末邮政局数-.059.251-.249年末固定电话用户.087.012.025社会商品零售总额.096-.022.099货物进出口.174-.095-.183年末实有公共汽.136-.029-.208剧场、影剧院.094-.072.197高等学校学生数-.061.041.454医院、卫生院

30、-.109.264-.041执业助理医生.034.073.049三废综合利用值-.002-.085.457通过表15,可知F1（存量因子）=-0.106*年末总人口+0.101*地区生产总值+0.028*客运量-0.033*货运量+0.119*地方财政预算内收入+0.091*地方财政预算内支出-0.033*固定资产投资总额+0.119*城乡居民储蓄年末余额+0.157*在岗职工平均工资-0.059*年末邮政局数+0.087*年末固定电话用户数+0.096*社会商品零售总额+0.174*货物进出口总额+0.136*年末实有公共汽车数+0.094*剧场、影剧院-0.061*普通高等学校学生数-0.

31、109*医院、卫生院+0.034*执业（助理）医师-0.002*三废综合利用值。F2和F3的得分也可类似求出。通过表16可见，各因子之间是正交的，即相互之间彼此独立。表16成分得分协方差矩阵成份得分协方差矩阵份成12311.000.000.0002.0001.000.0003.000.0001.000提取方法：主成份。旋转法：具有Kaiser标准化的正交旋转法。3.4 因子分析后续分析当我们获得了各城市在三个公因子的得分后，就可以利用因子得分进行分类和排序。我们可以利用计算变量和排序来进行各城市的经济综合得分计算和排血，具体操作步骤如下：11打开数据文件，依次单击“转换”一“计算变量”命令，弹出如图4所示的对话框。图4“计算变量”对话框在“目标变量”一栏中输入“综合得分”，这一变量将最终代表各个城市的综合经济实力。综合得分=0.47193*FAC1_1+0.29158*FAC2_1+0.07716*FAC3_1其中FAC1_1FAC2_1ffiFAC3_1是在做因子分析时提取的公因子保存变量，前面的系数是各个公因子的方差贡献率。单击“确定”，返