本文基于多元统计模型研究可持续发展竞争力评价问题，该模型是全局主成分分析和层次分析的有机结合，我们利用该模型对云南省2005~2015年社会经济可持续发展竞争力指标体系进行评价分析，并在此基础上得到相应的结论。 This paper studied the evaluation of sustainable development competitiveness based on the mul-tivariate statistical model, and it is an organic combination of global principal component analysis and analytic hierarchy process. We use this model to make comprehensive evaluation analysis on sustainable development competitiveness of social economy in Yunnan Province from 2005 to 2015. Finally the corresponding conclusions are given.
莫玉莲1，费宇2*
1云南师范大学商学院，云南 昆明
2云南财经大学统计与数学学院，云南 昆明
收稿日期：2018年3月4日；录用日期：2018年3月18日；发布日期：2018年4月3日
本文基于多元统计模型研究可持续发展竞争力评价问题，该模型是全局主成分分析和层次分析的有机结合，我们利用该模型对云南省2005~2015年社会经济可持续发展竞争力指标体系进行评价分析，并在此基础上得到相应的结论。
关键词 :可持续发展竞争力，指标体系，全局主成分分析
Copyright © 2018 by authors and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
中国的经济可持续发展是提高经济综合竞争力的基础，而经济可持续发展竞争力的提升是保持经济可持续发展的关键，故中国经济的发展方式必须是基于可持续发展竞争力提升的可持续发展。按照中央对云南省社会经济发展的要求，云南要做生态文明的排头兵和民族团结的进步的示范区 [
目前对社会经济可持续发展竞争力的评价采用的是多目标、多因素的综合评价方法，根据对各指标的相对重要性加以评价的定量评价方式，主要分为主观评价法和客观评价法 [
在对社会经济可持续发展竞争力客观评价方面，主要使用的是主成分分析综合评价，为了从一个主观和客观相结合的方面对社会经济可持续发展竞争力进行动态评价，本文采用主观评价方法中使用频率最高的层次分析法和全局主成分分析(AHP & GPCA)对云南省社会经济可持续发展竞争力实现系统层次和逻辑结构分明的量化动态评价 [
主成分分析是利用降维的思想，把多指标转化为少数几个综合指标(即主成分)，其中每个主成分都能够反映原始变量的大部分信息，且所含信息互不重复。这种方法在引进多方面变量(指标)的同时将复杂因素归结为几个主成分，使问题简单化，同时得到的结果更加科学有效的数据信息 [
因主成分分析一般用于研究静态的统计指标问题，其仅仅反映了该统计指标在某个静态时间上的发展状况，故此时有必要引入一种动态的统计分析方法来对我们所要研究的统计指标做分析，而全局主成分分析就是其中一种可以研究动态的统计指标的统计方法。
全局主成分分析与主成分分析的原理基本相同，它们都是在力保数据信息损失最少的原则下，把多个变量经过一定的变换，最终用少数几个具有代表性的变量来对所要研究的内容进行分析。
即若统计n个地区，使用相同的p个指标 X 1 , X 2 , ⋯ , X p 来描述，一共统计了t年的数据，经典主成分分析是对这p个静态的(一般指某一时间)某年的样本X (截面数据表)进行逐一评价，而全局主成分分析是将这t年的样本构成一张时序立体数据表K，然后再对整张时序立体数据表K执行与主成分分析一致评价的过程，其中
K = { X t ∈ R n × p , t = 1 , 2 , ⋯ , T }
层次分析赋权法(AHP)是主观赋权法的一种，即专业的相关研究人员根据自身所掌握的知识，如其常识和经验，评估评价指标对评价目标的相对重要程度并据此赋予评价指标相应的权重的一种赋权方法 [
层次分析赋权法的步骤：
1) 构造判断矩阵。将指标层每个指标因素下有联系的各具体指标进行两两比较，从而得到相对重要性程度(层次单排序)。用具体的标度值写出，组成矩阵形式，即判断矩阵。
设 A = ( b i j ) n × n = [ b 11 b 12 b 21 b 22 ⋯ b 1 j ⋯ b 2 j ⋮ ⋮ b i 1 b i 2 ⋱ ⋮ ⋯ b i j ] n × n = [ u 1 / u 1 u 1 / u 2 u 2 / u 1 u 2 / u 2 ⋯ u 1 / u j ⋯ u 2 / u j ⋮ ⋮ u i / u 1 u i / u 2 ⋱ ⋮ ⋯ u i / u j ] n × n
并满足 b j i = 1 / b i j ；i不等于j； i , j = 1 , 2 , 3 , ⋯ , n ； b i j > 0 。
式中：A为判断矩阵；n为两两比较的因素数目； b i j 为因素 u i 比 u j 项对某准则重要性的比例标度值。
本文标度值确定采用A. L. Sany (1980)提出的1-9标度法，见表1所示。
| 标度 | 定义 | 解释 |
|---|---|---|
| 1 | 同等重要元素 | i和j同样重要 |
| 3 | 略重要元素 | i和j略微重要 |
| 5 | 相当重要元素 | i比j重要 |
| 7 | 明显重要元素 | i比j明显重要 |
| 9 | 绝对重要元素 | i比j绝对重要 |
| 2,4,6,8 | 介于两相邻重要程度间 | |
| 倒数 | 元素i与j比较，得到bij，则元素j与i比较 | |
表1. 1-9标度法
2) 权重的计算。由判断矩阵计算各比较指标对某准则的相对权数。假设有n个指标 u 1 , u 2 , ⋯ , u n ，对于某准则的判断矩阵A，对于另一准则的相对权重为 W 1 , W 2 , ⋯ , W n 。相对权重的向量形式可以写成 W = ( W 1 , W 2 , ⋯ , W n ) T ，然后根据特征根求解判断矩阵的最大特征根与相应的特征向量，随后将特征向量正规化，便形成了相对权重向量。
计算相对权重向量的步骤如下：
1) 计算判断矩阵每一行元素的乘积 M i
M i = ∏ j = 1 n b i j , i = 1 , 2 , ⋯ , n
2) 计算 M i 的n次方根 W ¯ i
W ¯ i = M i n
3) 对向量 W ¯ = [ W ¯ 1 , W ¯ 2 , ⋯ , W ¯ n ] T 正规化
W i = W ¯ i ∑ j = 1 n W ¯ j
则 W = [ W 1 , W 2 , ⋯ , W n ] T 即为所求权重向量。
4) 一致性检验：以一致性比率 C R < 0.1 为判断准则，其中 C R = C I R I ， C I = λ max − n n − 1 ，n为判断矩阵A的阶数，RI为平均一致性随机指标， λ max = ∑ i = 1 n ( A W ) i n w i , ( i = 1 , 2 , 3 , ⋯ , n ) ， λ max 为判断矩阵的最大特征根，
W i 为因素i特征向量，即相对权重。特别值得注意的是因为考虑到专家在构造判断矩阵时，不可避免地会产生认识上的不一致，故需要对判断矩阵进行一致性检验。
首先设总体层为K，系统层的各指标为 A i , ( i = 1 , ⋯ , o ) ，状态层的各指标为 B j , ( j = 1 , ⋯ , m ) ，变量层的各指标为 C k , ( k = 1 , ⋯ , n ) ，要素层的各指标为 D l t , ( l = 1 , ⋯ , s , t = 1 , ⋯ , q ) ；其次将状态层的指标权重记为 ω 1 ( 1 ) , ω 2 ( 1 ) , ⋯ , ω m ( 1 ) ，变量层的指标权重记为 ω 1 ( 2 ) , ω 2 ( 2 ) , ⋯ , ω n ( 2 ) ，要素层的指标权重记为 ω 1 ( 3 ) , ω 2 ( 3 ) , ⋯ , ω s ( 3 ) ，然后将标准化的指标值 D l t 通过 D l ˜ t = ( 1 + ω l ( 3 ) ) D t , t = 1 , ⋯ , q , l = 1 , ⋯ , s 变换为新的指标值 D l ˜ t ，最后形成新的数据表 K ˜ 。
将得到的数据表 K ˜ 进行全局主成分分析，得到变量层的评价指数 f k , ( k = 1 , 2 , ⋯ , n ) ；再根据AHP赋权所得到的变量层的权重系数，对变量层的评价指数加权求和得到状态层的得分指数 B j = w k ( 2 ) f k , k = 1 , 2 , ⋯ , m ，以此类推得到系统层的得分指数 F i = w j ( 1 ) B j , j = 1 , ⋯ , o ，最后把系统层的得分指数求平均得到总体层的得
分指数 Z = 1 p ∑ i = 1 p F i 。
为验证以上模型的有效性和实用性以及其具备的实践意义，本文选取了前期工作中所构建的云南省可持续发展竞争力指标体系(见表2所示)作为评价对象 [
| 总体层 | 系统层 | 状态层 | 变量层 | 要素层 |
|---|---|---|---|---|
| 云南省社会经济可持续发展竞争力 | 生存支持系统 | 生存资源禀赋 | 耕地资源指数、水资源指数、气候资源指数 | 138个指标 |
| 农业生产投入 | 资金投入指数 | |||
| 资源转化效率 | 生物转化效率指数、经济转化效率指数 | |||
| 生存持续能力 | 生存稳定指数、生存潜力指数、生存持续指数 | |||
| 发展支持系统 | 基础设施水平 | 交通基础设施指数、信息基础设施指数 | ||
| 经济发展水平 | 经济规模指数、经济结构指数、经济市场化指数、经济国际化指数、经济要素投入指数 | |||
| 经济效益水平 | 工业企业经济效益指数、工业企业运营效益指数 | |||
| 资源环境绩效 | 区域经济资源环境绩效指数 | |||
| 环境支持系统 | 环境质量 | 空气污染指数、水污染指数、固体废弃物污染指数 | ||
| 生态水平 | 生态脆弱性指数、生态质量指数 | |||
| 环境保护力度 | 环境污染治理指数、资源节约与循环利用指数、能源结构调整指数、生态保护建设指数 | |||
| 社会支持系统 | 人口发展水平 | 人口增长指数、人口素质指数、人口结构指数 | ||
| 居民生活质量 | 消费水平指数、居住条件指数、出行便利度指数、配套设施水平指数 | |||
| 公共卫生服务 | 公共卫生投入指数、公共卫生设施指数、公共卫生服务能力指数 | |||
| 社会安全水平 | 社会公平指数、社会安全指数、社会福利保障指数 | |||
| 文化发展水平 | 文化事业投入指数、文化设施服务指数 | |||
| 智力支持系统 | 区域教育水平 | 教育投入指数、教育规模指数、教育成就指数 | ||
| 区域创新能力 | 创新资源指数、创新效率指数、创新效益指数 | |||
| 区域管理能力 | 政府工作绩效指数、经济调控绩效指数、社会管理绩效指数、资源环境管理绩效指数 |
表2. 云南省社会经济可持续发展竞争力指标体系
通过查阅云南省2006年到2016年的《云南统计年鉴》和《中国统计年鉴》以及其他相关的统计年鉴 [
为了消除量纲对评价过程和结果的影响，基于层次分析赋权的全局主成分分析采用了直线型 Z − Kore 法对原始数据进行标准化处理。
本文首先运用专家打分的方式，将指标层每个指标因素下有联系的各具体指标进行两两比较，从而得到相对重要性程度(层次单排序)。用具体1-9标度值(如表1所示)写出，组成判断矩阵。然后检验各判断矩阵的一致性，最后根据权重的计算方法算出各层指标下相应因素的权重大小。
云南省社会经济可持续发展竞争力指标中，通过MCE (Modern Comprehensive Evaluation)软件，分别构造了53个变量层、19个状态层和5个系统层的判断矩阵，并进行了相应的一致性检验，其中五个系统层的判断矩阵分别为：
1) 生存支持系统状态层的判断矩阵
A 1 = ( 1 3 2 1 1 / 3 1 1 / 2 1 / 3 1 / 2 2 1 1 / 2 1 3 2 1 )
其中 C R = 0.0038 , λ max = 4.0104 , C I = 0.035 , R I = 0.9 ，故判断矩阵A1是一致的，可以根据其进行生存支持系统状态层的指标权重计算。
2) 发展支持系统状态层的判断矩阵
A 2 = ( 1 1 / 3 1 / 2 1 / 2 3 1 2 2 2 1 / 2 1 1 2 1 / 2 1 1 )
其中 C R = 0.0038 , λ max = 4.0104 , C I = 0.035 , R I = 0.9 ，故判断矩阵A2是一致的，可以根据其进行发展支持系统状态层的指标权重计算。
3) 环境支持系统状态层的判断矩阵
A 3 = ( 1 1 1 1 1 1 1 1 1 )
其中 C R = 0 , λ max = 3 , C I = 0 , R I = 0.58 ，故判断矩阵A3是一致的，可以根据其进行环境支持系统状态层的指标权重计算。
4) 社会支持系统状态层的判断矩阵
A 4 = ( 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 )
其中 C R = 0 , λ max = 5 , C I = 0 , R I = 1.12 ，故判断矩阵A4是一致的，可以根据其进行社会支持系统状态层的指标权重计算。
5) 智力支持系统状态层的判断矩阵
A 5 = ( 1 1 2 1 1 2 1 / 2 1 / 2 1 )
其中 C R = 0 , λ max = 3 , C I = 0 , R I = 0.58 ，故判断矩阵A5是一致的，可以根据其进行智力支持系统状态层的指标权重计算。
此时总体层即云南省社会经济可持续发展竞争力为A，系统层为五大支持系统，其各指标为
A i , ( i = 1 , ⋯ , 5 ) ，状态层的各指标为 B j , ( j = 1 , ⋯ , 19 ) ，变量层的各指标为 C k , ( k = 1 , ⋯ , 53 ) ，要素层的各指标为 D l t , ( l = 1 , ⋯ , 138 , t = 1 , ⋯ , 10 ) ；其次将状态层的指标权重记为 ω 1 ( 1 ) , ω 2 ( 1 ) , ⋯ , ω 19 ( 1 ) ，变量层的指标权重记为 ω 1 ( 2 ) , ω 2 ( 2 ) , ⋯ , ω 53 ( 2 ) ；要素层的指标权重记为 ω 1 ( 3 ) , ω 2 ( 3 ) , ⋯ , ω 138 ( 3 ) ，然后将标准化的指标值 D l t 通过 D l ˜ t = ( 1 + ω l ( 3 ) ) D t , t = 1 , ⋯ , 10 , l = 1 , ⋯ , 138 变换为新的指标值 D l ˜ t ，最后形成新的数据表 A ˜ 。
用上面得到的数据表 A ˜ 进行全局主成分分析，得到变量层的评价指数 f k , ( k = 1 , 2 , ⋯ , 53 ) ；再根据AHP赋权所得到的变量层的权重系数，对变量层的评价指数加权求和得到状态层的得分指数
B j = w k ( 2 ) f k , k = 1 , 2 , ⋯ , 53 ，以此类推得到系统层的得分指数 F i = w j ( 1 ) B j , j = 1 , ⋯ , 19 ，最后把系统层的得分指数求平均得到总体层的得分指数 Z = 1 5 ∑ i = 1 5 F i 。
通过上诉过程可得云南省社会经济的五大支持系统和可持续发展竞争力的得分函数分别为：
1) 生存支持系统得分函数
F 1 = 0.3512 B 1 + 0.1089 B 2 + 0.1887 B 3 + 0.3512 B 4
2) 发展支持系统得分函数
F 2 = 0.1221 B 5 + 0.4231 B 6 + 0.2274 B 7 + 0.2274 B 8
3) 环境支持系统得分函数
F 3 = 1 3 B 9 + 1 3 4 B 10 + 1 3 B 11
4) 社会支持系统得分函数
F 4 = 0.2 B 12 + 0.2 B 13 + 0.2 B 14 + 0.2 B 15 + 0.2 B 16
5) 智力支持系统得分函数
F 5 = 0.4 B 17 + 0.4 B 18 + 0.2 B 19
6) 云南省社会经济可持续发展竞争力综合得分函数
Z = 0.2 F 1 + 0.2 F 2 + 0.2 F 3 + 0.2 F 4 + 0.2 F 5
根据以上得分函数，可以得到如表3的云南省社会经济的五大支持系统及可持续发展竞争力的2005年至2015年综合得分。
如表3所示，云南省的社会经济可持续发展竞争力的得分逐年提升，其五大支持系统在2005年至2015年的得分也呈现的是逐年增加，这与2005年至2015年的云南省社会经济发展状况是一致的。
如图1所示，不管是由主成分分析综合评价得到的结果，还是由基于层次分析的全局主成分分析综合评价方法得到的结果，都显示了2005年~2015年云南省的社会经济可持续发展竞争力呈现的是一个平稳的上升趋势。
本文使用主客观相结合的基于层次分析的全局主成分分析评价方法对云南省2005年~2015年的社会经济可持续发展竞争力指标体系做了综合评价，得出了相应的评价结果(如表3所示)，并与客观的主成分
| 年份 | 生存支持系统 | 发展支持系统 | 环境支持系统 | 社会支持系统 | 智力支持系统 | 可持续发展竞争力 |
|---|---|---|---|---|---|---|
| 2005年 | −0.442 | −0.770 | −1.123 | −0.448 | −1.754 | −0.907 |
| 2006年 | 0.251 | −0.937 | −0.738 | −0.699 | −0.778 | −0.580 |
| 2007年 | −0.307 | −0.555 | −0.620 | −0.620 | −1.083 | −0.637 |
| 2008年 | −0.037 | −0.653 | −0.733 | −0.275 | −0.384 | −0.416 |
| 2009年 | −0.785 | −0.228 | −0.373 | −0.152 | −0.273 | −0.362 |
| 2010年 | −0.204 | −0.065 | 0.044 | 0.160 | 0.247 | 0.036 |
| 2011年 | −0.170 | −0.038 | 0.981 | −0.026 | 0.668 | 0.283 |
| 2012年 | 0.142 | 0.505 | 0.554 | 0.302 | 0.896 | 0.480 |
| 2013年 | 0.498 | 0.853 | 0.524 | 0.182 | 0.749 | 0.561 |
| 2014年 | 0.490 | 0.868 | 0.735 | 0.562 | 0.830 | 0.697 |
| 2015年 | 0.563 | 1.021 | 0.750 | 0.725 | 0.882 | 0.788 |
表3. 层次分析赋权的全局主成分分析的评价结果
图1. 两种方法的云南省社会经济可持续发展竞争力变化趋势
分析评价方法做了相应的对比(如图1所示)，最后得出了如下几点结论和展望：
1) 两种方法都显示了2005年至2015年云南省社会经济可持续发展竞争力平稳上升。
2) 五大支持系统的综合得分总体上逐年提升。其中生存支持系统的得分的变化在2009~2011年下降的幅度较大，这主要是因为2008年的雪灾等一系列自然灾害和2010~2011年的云南省大干旱所导致的水资源量减少和农田成灾率的增加所引起的。发展支持系统总体呈现一个较为平稳的上升趋势。这反映了云南在政府的一系列政策下，云南的经济快速平稳地增长，社会各项基础设施也充分地完善。环境支持系统经历了一个曲折上升的过程。在2006~2009年，环境支持系统呈现上升的趋势，这充分反映了“十一五”期间，云南省节能减排工作取得实质性的进展，有效地缓解和遏制了云南经济快速发展中所带来的生态环境冲击。社会支持系统也呈现了一个曲折的上升趋势。这显示了云南省社会迈向全面发展并进入到良性发展的轨道，也说明了政府的一系列政策使得云南人民的生活变得更加美好。智力支持系统呈现波动上升的趋势。证明政府对云南省科研经费的增加和2007年“两免一补”等政策的实施，使得云南省在教育、科技创新和政府管理方面取得较大的进步。
3) 由于数据来源的局限性，如相应的统计年鉴只出版至2016年的《云南统计年鉴》和《中国统计年鉴》等，故本文只选取了2005年~2015年的云南省社会经济可持续发展竞争力指标体系做综合评价，在今后的相关研究中也可对数据部分进行补充。
云南省哲学社会科学研究基地2015重点项目“云南省社会经济可持续发展竞争力指标体系研究”(JD2015ZD20)。
莫玉莲,费 宇. 基于全局主成分分析的云南省可持续发展竞争力评价 Evaluation of the Sustainable Development Competitiveness of Social Economy in Yunnan Province Based on Global Principal Component Analysis[J]. 可持续发展, 2018, 08(02): 85-93. https://doi.org/10.12677/SD.2018.82010