基于ESDA的山东省四维城镇化协调发展空间差异
吴云青，王多多，郭 佳
（天津工业大学经济与管理学院，天津300387）
摘要： 在构建评价指标体系的基础上，运用多指标综合评价法、协调发展度模型、探索性空间分析等方法探讨山东省人口、土地、经济与社会城镇化协调发展的空间差异。结果表明，大多数地市的城镇化都处于中度及以上协调发展状态；人口、土地、经济和社会城镇化协调发展度存在显著的正空间自相关特点，城镇化协调发展的均质性区域明显多于异质性区域，城镇化协调发展的空间依赖性和差异性同时存在；H-H聚集类型区主要集中在东部沿海和鲁中地区，H-L聚集类型区内只有济南市；L-L聚集类型区、L-H聚集类型区主要分布在鲁西和鲁南地区。
关键词： 城镇化；协调发展；空间差异；探索性空间分析
长期以来，我国的城镇化采用传统模式，质量不高，存在诸多限制性因素，尤其表现为城镇化内部各要素不协调、不匹配［1］，经济城镇化不稳，土地城镇化冒进，人口城镇化滞后，社会城镇化脆弱［2］。党的十八届三中全会提出要坚持走中国特色新型城镇化道路，而新型城镇化是内涵式发展的城镇化，离不开其内部要素相互协调、良性互动，城镇化内部各要素能否协调发展及其协调发展程度决定了新型城镇化进程的可持续性及其效益。国内相关研究大多将城镇化划分为人口、土地、经济和社会城镇化4个要素；国外并没有我国这种“要素式”城镇化内涵界定与分类，更没有城镇化内部各要素协调发展的研究，但他们对人口迁移［3］、城市空间扩张［4］、城市人口增长［5］的研究成果可为本研究提供有益借鉴。近年来，着眼于“健康、有序、绿色、可持续”城镇化，田时中等对人口、土地、经济和社会城镇化协调发展状况开展了系列研究，但这些研究多停留在协调度测量层面，聚焦于要素协调发展空间差异的研究并不多见［6-8］。另外，已有为数不多的城镇化协调发展空间差异的研究也大多采用传统分析方法，忽视了空间因素的作用，缺乏对空间信息的有效体现，对城镇化协调发展在地理位置上的关联度、聚集性的研究较少，难以揭示人口、土地、经济与社会城镇化协调发展的空间特征。因此，本研究以山东省城镇化发展的实际状况为例，构建人口、土地、经济和社会城镇化各要素评价指标体系，利用协调发展度模型、探索性空间分析（ESDA）研究山东省人口、土地、经济和社会城镇化协调发展的空间差异，这对于有针对性地制定创新相关政策，促进新型城镇化健康、可持续发展，提升新型城镇化质量具有重要的现实意义。
1 评价指标体系构建及研究方法
1.1 指标体系构建
基于研究目的，遵循评价指标构建的系统性、科学性、代表性、可比性等原则，利用理论分析法、频度统计法、专家咨询法甄选人口、土地、经济、社会城镇化评价指标体系。为了能够获得客观的原始数据，尽可能选择统计年鉴、统计公报中有统计数据的指标或依据统计年鉴、统计公报经数学计算能够得到数据的指标。为便于比较不同区域城镇化发展状况，指标尽量少采用总量指标，较多采用人均和相对指标，构建评价指标体系（表1）。
1.2 研究方法
1.2.1 人口、土地、经济与社会城镇化发展评价方法
1.2.1.1 基于组合赋权法的评价指标权重计算 本研究将主观赋权的三角模糊法与客观赋权的熵值法结合计算指标组合权重。按照模糊数学原理，三角模糊法计算的标权重取决于专家的模糊判断。首先，由专家给出最悲观估计、最可能估计和最乐观估计，形成三角模糊合成矩阵；其次，运用算子M（⊗，⊕）进行模糊合成；最后，根据各指标的模糊得分计算出最终权重，所得权重向量可表示为：
在传统金融系统中，账户名是可以公开的，但账本内容是必须保密的。为储户保密指账户中的记录，不是账户号。区块链的设计反过来了，账户号是匿名的，但账本内容是公开的。区块链匿名性说的是账户名匿名，不是账户中的记录。
式中：u1j表示利用三角模糊法计算得到的第j个指标的权重。
利用熵值法计算指标权重，首先要计算出指标的熵值；其次，计算出指标各自对应的变异系数；最后，确定各个指标的权重，所得权重向量记为：
假设各个指标的组合权重为：
2.2.3 城镇化协调发展度计算 利用公式（8）至公式（10）可计算出山东省2015年综合城镇化指数、城镇化耦合系数和城镇化协调发展度（表2）。本研究将山东省人口、土地、经济与社会城镇化协调发展度（E）划分为4个等级：0＜E≤0.04为低度协调；0.4＜E≤0.6为中度协调，0.6＜E≤0.8为良好协调，0.8＜E≤1为优度协调。
面对公摊引发的诸多问题，一些专家认为，多年来我国法律界就有声音认为“公摊争议严重影响《物权法》有效实施”，广州、甘肃、北京、重庆都曾出台征求意见稿，尝试推行套内面积计价。艾振强认为，采取“套内面积”计价的主要好处在于“让房产交易更加透明，最大程度减少后续纠纷，保护购房者合法权益”，同时最大程度避免公摊面积给房产税征收埋下“不公平”隐患。
本研究在对运用相对熵求解贴近度算法进行改进的基础上计算组合系数θk，运用求最优解的方法获取1个与单个赋权方法所得权重结果最贴近的集结权重，记为X＝（x1，x2，…，xN）。
表1 人口、土地、经济、社会城镇化评价指标体系及权重
假设U＝（u1，u2，…，uQ）为各个单一赋权方法所得权重的集合，ukj为基于第 k种方法的第j项指标的权重（k＝1，2，…，Q；j＝1，2，…，N），运用相对熵思想构建优化模型。
这个故事对我触动很大，从小父母对我进行的安全教育中，说的最多的就是“不要和陌生人说话”。自己当了母亲后，我对孩子的安全教育很重视，但是我一直在反思，作为安全教育的这句话到底对不对呢？
求得公式（4）的最优解为：
1.2.1.3 人口、土地、经济与社会城镇化指数计算 根据各指标标准化值及其对应的组合权重，运用多指标综合评价法分别计算人口、土地、经济与社会城镇化指数值。
计算各个单一赋权所得结果uk与集结权重X＝（x1，x2，…，xN）的贴近度 h［uk，x］（k＝1，2，…，Q），再依据贴近度计算可信度即组合系数［9］。
1.2.1.2 指标数据的标准化处理 采用极差标准化方法对各个评价指标数据进行规范化处理，从而得到标准化值［10］。
2.2.1 评价指标权重的计算 首先，分别运用三角模糊法、熵值法计算指标权重，得到权重向量u1j、u2j；其次，利用公式（4）至公式（6）可以得到人口、土地、经济与社会城镇化评价中三角模糊法和熵值法的组合系数θk分别为0.364、0.636，0.198、0.802，0.244、0.756，0.194、0.806；最后，依据公式（3）可以计算得到各个指标的组合权重（表1）。
陈诚的警卫排长夏国忠早已在门外等候，听见总司令叫他，赶紧几步来到门口，朗声报告。陈诚招招手，示意让他进去。夏国忠走进接待室，给陈诚敬了礼，又向胡琏敬了一个礼。说“胡长官好。”
纹枯病是水稻生长过程中较为常见的一种疾病，主要是在稻秆与稻叶中，早期症状是稻秆及稻叶上有圆形斑点，中后期圆形斑点不断增大并伴有白色菌丝，随着会病情加重会逐渐变为褐色菌斑。发生纹枯病的主要原因是在进行水稻种植时没有根据栽种标准进行栽培，导致养分不足或者过量，为此纹枯病需要通过以下几个方面进行处理：1.应严格进行田间管理，科学提供养分与灌溉，这在较大程度上能够避免真菌的侵袭；2.适量农药的使用能够有效杀灭有害菌，可以使用井冈美素A可湿粉剂进行喷洒治疗，对治疗纹枯病有较好的治疗效果。
Moran'sI指数取值介于［-1，1］，当 Moran'sI＜0且显著时，表示存在负空间自相关；当Moran'sI＞0且显著时，表示存在正空间自相关；当Moran'sI＝0且显著时，表示不存在空间自相关。
式中：P表示综合城镇化指数；D表示城镇化耦合系数；E表示城镇化协调发展度；F1、F2、F3、F4分别表示人口、土地、经济、社会城镇化指数；a、b、c、d分别表示人口、土地、经济、社会城镇化对城镇化总系统的贡献度，考虑到各子系统对城镇化发展同等重要，故取 a＝b＝c＝d＝1／4。
其次，我国应当继续保持对黑社会性质组织犯罪的高压态势。在持续的高压打击下，当前我国黑社会性质组织犯罪的发展势头受到了一定的遏制，但我们不能因此而掉以轻心。研究表明，我国大陆过半数的黑社会性质组织的存在时间都不长，一般是1-5年，但具有升级演变快、再生能力强的特点。近几年在一些地区，一年半载就能形成较成熟的黑恶势力，还存在一定区域或行业打掉一个组织不久就形成一个新的犯罪组织替代的现象，也就是所谓的“割韭菜”现象明显。［18］因此，在打击黑恶犯罪时，一定要对其保持高压态势，坚持“打早打小、露头就打”的方针，防止其成型坐大，泛滥成灾。
1.2.3 探索性空间数据分析（ESDA）
1.2.3.1 全局空间自相关分析 全局空间自相关体现了观测变量在研究区域内空间相关性的总体特征，可以定量表征某种地理现象在空间上的总体关联与差异程度。本研究运用Moran提出的GlobalMoran'sI指数衡量全局自相关情况［12］，计算公式为：
式中：n表示区域数量；xi、xj分别表示观测变量x在第i、第j个研究区域的观测值；
表示观测变量x的均值；wij表示研究区域之间的空间权重，如果研究区域i与j相邻，则wij＝1，否则 wij＝0。
运用标准化Z对Moran'sI进行显著性检验，检验公式为：
式中：Z（I）表示空间自相关的显著性水平；E（I）表示 Moran's I的数学期望；VAR（I）表示Moran'sI的方差。
1.2.2 人口、土地、经济与社会城镇化协调发展度计算 引入协调发展度模型对山东省人口、土地、经济与社会城镇化协调发展情况进行定量评价［11］。
1.2.3.2 局部空间自相关分析 局部空间自相关可以分析观测变量的局部空间变化，能具体测量每个区域与周边区域的局部空间关联和空间差异程度，用Moran散点图和Local Moran'sI（LISA）指数来衡量。Moran散点图的横轴是变量在不同位置的观测值向量，纵轴是该向量空间滞后向量的所有取值，每个空间滞后是该区域周边观测值的加权平均值，Moran散点图分为4个象限。然而，Moran散点图无法定量揭示各研究区域与其相邻区域的自相关程度及其显著性，Local Moran'sI（LISA）指数恰恰可弥补这一缺陷，可用于定量分析研究区域i与j之间的空间自相关性，计算公式为：
80例患者经超声心动图检查结果显示EKG与典型ST段抬56例，EKG表现正常2例，完全性左束支传导阻滞3例，非ST段抬高的患者且有轻微的ST-T变化7例，金标准数据显示为61例、4例、5例、10例。与金标准相比，80例患者经超声心动图检查结果显示，灵敏度为85.00%，特异度为0.00%，准确率为85.00%。见表1所示。
2 实证结果与分析
2.1 研究区域与数据来源
作为全国的农业大省、工业大省和经济大省，山东省近年来城镇化进程不断推进。但山东省城镇化也具有“速度快、协调性不够、区域差异明显”的特点，其城镇化进程具有较强的现实背景和代表性，所得研究结论和启示对于其他地区的城镇化协调发展具有重要的参考意义。
本试验是基于2015年山东省地级市尺度的研究，共有17个研究单元，数据主要来源于2016年《中国城市统计年鉴》、2016年《山东省统计年鉴》、山东省2015年国民经济和社会发展统计公报等。
2.2 城镇化协调发展度计算
吵也好，打也好，真真假假地闹过离婚也好，齐海峰和李红磕磕绊绊地走过了几年，跌跌撞撞地又走了几年。这一走就走了大半辈子。
2.2.2 城镇化指数计算 在对各评价指标数据进行标准化处理的基础上，可计算出2015年山东省人口、土地、经济与社会城镇化指数值（表2）。
利用采样点集SAB、SBC、SDE、SEF,对公式(4)应用最小二乘法,可获得k、bAB、bBC、bDE、bEF的估计值,得到以上4条直线在世界坐标系中的直线方程。
表2 山东省人口、土地、经济与社会城镇化指数与协调发展度评价结果
式中：wj表示组合权重；θk表示组合系数；ukj表示第k种赋权方法计算的第j个指标的权重值。
选取2016年7月～2018年5月在我院接急诊PCI术治疗的急性心肌梗死伴泵衰竭的高危患者28例作为研究对象。入选标准：所有患者均符合ST段抬高型心肌梗死（STEMI）诊断标准、泵功能衰竭Killip分级Ⅱ～Ⅳ级者；行急诊PCI术者；排除标准：肿瘤，肝、肾和血管等严重器质性疾病者，有心脏手术史者。依据护理方式差异将其分为试验组和对照组。其中，试验组男14例，女6例，年龄46～81岁，平均年龄（62.1±3.6）岁；对照组男21例，女7例，年龄45～80岁，平均年龄（65.1±2.7）岁。
式中：
分别表示观测变量 x在第 i、j研究区域的观测值；wij表示研究区域之间的空间权重矩阵，如果研究区域i与j相邻，则wij＝1，否则wij＝0
表示观测变量x的均值。统计检验也采用标准化Z检验。如果Moran's Ii＞0且显著，表示某研究区域与周围区域之间差异小，存在高-高值或低-低值空间聚集；如果Moran'sIi＜0且显著，表示研究区域与周围区域之间差异大，存在高（H）-低（L）值或低-高值空间聚集。将这些通过显著性检验的区域可视化到地图上即可形成LISA集聚图，可测量局部空间的异质性，识别局部空间聚集的“热点”与“冷点”。
2.3 城镇化协调发展度计算结果分析
以山东17个市作为横轴标签，以表2中人口城镇化指数、土地城镇化指数、经济城镇化指数、社会城镇化指数、综合城镇化指数为主纵坐标轴，以城镇化耦合系数、城镇化协调发展度为次纵坐标轴，可得到2015年山东省17个市城镇化协调发展度变量的二维坐标图（图1）。
由图1可知，山东17个市之间城镇化耦合度、协调发展度差异明显，但耦合系数和协调发展度2条折线的折拐方向基本一致，说明山东省城镇化耦合度、协调发展度的空间布局特点较相似。济南市城镇化协调发展水平最高，其次依次是青岛市、东营市和威海市，其协调发展度均大于0.8，属于优度协调发展。这4个市不仅各系统的城镇化水平在全省位居前列，且发展较均衡，无明显短板。淄博、烟台、潍坊、莱芜等市城镇化协调发展度均大于0.6但小于0.8，属于良好协调发展。其中，淄博市土地城镇化水平较高，位居全省第2位，人口、经济和社会城镇化水平均居全省上游；烟台市土地、经济、社会城镇化水平居全省中上游且发展较均衡，但人口城镇化发展相对滞后，列全省第11位。潍坊市、莱芜市均属于在城镇化某单一系统表现突出，但在其他方面发展较差。滨州、日照、泰安、枣庄、济宁等市处于中度协调发展状态，其各系统城镇化水平大多位列全省中下游水平，但发展较均衡。临沂、德州、聊城、菏泽等市处于低度协调发展状态，其各维度城镇化水平都处于中下游位置，其中菏泽市更是各方面都居末位。
2.4 城镇化协调发展的空间差异分析
2.4.1 全局空间自相关分析 基于上述协调发展度计算结果，依据公式（10）并借助OpenGeoDa软件可以计算出2015年山东省人口、土地、经济与社会城镇化协调发展度的Global Moran'sI指数为0.191，运用公式（11）的Z值法对空间自相关进行显著性检验，结果表明GlobalMoran'sI指数的正态统计量Z值为2.00，大于0.05置信水平的临界值1.96，显著性检验通过。说明山东省人口、土地、经济和社会城镇化协调发展度存在显著的正空间自相关特点。
2.4.2 局部空间自相关分析 上述全局空间自相关分析无法有效识别空间聚类的内部状态，也可能会隐藏偏离整体布局特征的空间关系。本研究利用Moran散点图和Local Moran'sI（LISA）指数进一步分析山东省人口、土地、经济和社会城镇化协调发展的局部空间特征。以城镇化协调发展度的标准化向量
为横坐标，以城镇化协调发展度滞后向量W＿为纵坐标，以散点横纵坐标的均值为中心坐标，可以得到包括4个基本象限的Moran散点图（图2），每个象限对应不同的城镇化协调发展度局部空间聚集类型。由图2可知，山东省人口、土地、经济和社会城镇化协调发展度散点主要分布在H-H聚集类型区、L-L聚集类型区，约占全部散点的70.59%，而H-L聚集类型区、L-H聚集类型区的散点较少，约占全部散点的29.41%，山东省城镇化协调发展的均质性区域明显多于异质性区域，这也验证了上述全局空间自相关的分析结果。
利用OpenGeoDa软件计算出2015年山东省城镇化协调发展度的LISA值，在Z检验的基础上（P＜0.05）绘制出LISA聚集图（图3）。由图3可知，山东省大部分地市的城镇化协调发展度与相邻市之间的关联性并不显著，只有青岛市、烟台市、济南市和济宁市的LISA值通过5%显著性水平的Z统计检验，相对于其他地市，这4个市对山东省城镇化协调发展度空间自相关程度的贡献更大。
（1）H-H聚集类型。属于该类型的地市共有7个，分别是青岛市、东营市、威海市、淄博市、烟台市、潍坊市、莱芜市，以面状形式连片分布在山东半岛区域和鲁中地区。总体来说，该区域内各市自身的城镇化协调发展水平较高，周围相邻区域的城镇化协调发展度也较高，成为高值集簇区。但值得注意的是，尽管从全省层面看，该区域内各市城镇化协调发展度较高，但这7个市也存在较大差异，青岛、东营、威海等市各维度城镇化指数均位居全省前列，协调发展度均列全省前4位，但潍坊市和莱芜市城镇化各系统发展却差异较大，优点和缺点均较明显，这2个市城镇化协调度比较靠近全省协调度均值线，随时都有可能落入L-H聚集类型区。LISA指数分析结果表明，该聚集类型区的7个市中，只有青岛市和烟台市的LISA值通过了显著性检验。一方面，说明青岛市和烟台市本身城镇化协调发展能力较强，向东、西2个方向辐射到了威海、东营、潍坊、淄博等市；另一方面，也说明青岛市和烟台市与周边各市同处山东半岛沿海和胶济铁路沿线，空间联系密切、文化风俗相似，在城镇化进程中相互协作，良性竞争，共同推动青岛市和烟台市成为城镇化协调发展的“热点”区。潍坊、东营、淄博、莱芜等市的LISA值没有通过显著性检验，说明这4个地市城镇化协调发展能力还较弱，还不能有效辐射带动周边的滨州、日照、临沂、泰安等市区域城镇化的协调发展。
（2）L-H聚集类型。属于该类型的市有4个，分别是滨州市、泰安市、德州市和日照市，以点状形态零散分布在鲁西和鲁南地区。从历史基础条件来看，这4个区域一直都是山东省的农业大市，工业化、城镇化基础较差，本身的城镇化协调发展度较低；而其周边的东营、淄博、济南、莱芜、潍坊等市尽管城镇化协调发展度较高，但其本身的辐射力还不够强大，对于周边区域的带动能力有限，因此这4个市成为城镇化协调发展的低高空心区。
（3）L-L聚集类型。属于该类型的市有5个，分别是临沂市、枣庄市、济宁市、聊城市、菏泽市，以面状形态分布在鲁西和鲁南地区。该区域内各市自身的城镇化协调发展水平较低，而相邻区域的水平也较低，相对落后的聚集特征明显，成为低值萧索区。一直以来，这些市都是山东省的欠发达地区，工业化基础薄弱，城镇化起步晚，产业支撑不足；这些市在地理位置上均是三面环绕山东省的欠发达地区，而南部又与河北省、河南省、江苏省的欠发达市相邻，缺乏发达城市体的辐射带动。LISA指数分析结果表明，该类型区的5个市中，只有济宁市的LISA值通过了显著性检验，成为“冷点”区域，主要是因为相对于其他处于L-L聚集类型区的地市，济宁市的周边为泰安市、聊城市、菏泽市、枣庄市、临沂市，而这些地市均为城镇化协调度的低值区域，低水平聚集的特征更加明显。
（4）H-L聚集类型。属于该类型的市只有济南市，其LISA值也通过了显著性检验，说明济南市自身城镇化协调发展水平较高，位列全省第1位，但其周边的德州、聊城、泰安、滨州等市协调发展度却较低，成为高值孤立区。这种空间相关格局表明济南市城镇化发展的极化效应显著，长期以来济南市都是山东省的政治、文化中心，但其城镇化协调发展能力并没有强大到可以有效带动周边区域，与相邻市没有形成“近域扩散”效应；相反，由于自身的优势地位，济南市吸引了周边区域的人才、资本、技术向其聚集，“马太效应”明显，影响了相邻区域的发展。
3 结论与建议
本研究基于2015年山东省的相关数据，运用ESDA分析探讨人口、土地、经济与社会城镇化协调发展的空间差异，结果表明，山东省大多数地市的城镇化都处于中度及以上协调发展状态；山东省人口、土地、经济和社会城镇化协调发展度存在显著的正空间自相关特点；山东省城镇化协调发展的均质性区域明显多于异质性区域，城镇化协调发展的空间依赖性和差异性同时存在，东部沿海和西部内陆不均衡现象突出。H-H聚集类型区内有7个市，主要集中在山东东部沿海和鲁中地区，是城镇化协调发展的“热点”地区；H-L聚集类型区域内只有济南市；L-L聚集类型区内有5个市，主要分布在鲁西和鲁南地区，是城镇化协调发展的“冷点”地区；L-H聚集类型区内只有4个市，且其显著性不明显，零星散布在鲁西和鲁南地区。
综上，提出如下政策建议：第一，践行科学的城镇化理念，包括协调发展理念、共享理念、顺应市场规律理念。第二，实施差异化城镇化协调发展战略，提升山东半岛城市群各市的城镇化协调发展能力；加强宏观调控，扶持鲁西、鲁南各市城镇化协调发展。第三，创新保障城镇化协调发展的体制机制，如户籍制度、城乡一体化的社会保障体系、进城农民就业培训体系、集体建设用地入市制度等。
大学生综合素质测评是一个重要且现实意义很强的工作。本文在全面回顾现有综合测评理论和方法的基础上，引入关键控制点理论，构建了基于关键控制点的综合素质测评指标体系，并运用标准差法和理想点法进行综合评价。论文的研究成果可以为高校合理的奖惩、促进和决策提供相应的帮助。论文今后将依据评价结果对部分大学生的综合素质结果进行深度挖掘，从而提高现实应用价值。
参考文献：
［1］林爱文，樊 星.湖北省人口城镇化与土地城镇化协调发展分析［J］.地域研究与开发，2015，34（6）：14-18.
［2］姚士谋，陆大道，陈振光，等.顺应我国国情条件的城镇化问题的严峻思考［J］.经济地理，2012，32（5）：1-6.
［3］Diebold W，Despard L E.The evolution of the international economic order［J］.International Economic Policies and Their Theoretical Foundations，1992，9（4）：3-25.
［4］The price effects of international airline alliances［J］.Journal of Law and Economics，2000，43（2）：503-546.
［5］Henderson V.The urbanization process and economic growth：the sowhat question［J］.Journal of Economic Growth，2003，8（1）：47-71.
［6］田时中，方眉玉，李光龙.人口-土地-经济-社会城镇化水平测度及耦合协调分析——基于31个省份面板数据的实证［J］.河北地质大学学报，2017，40（4）：49-57.
［7］多淑杰.城镇化四维分析模型的构建与协调发展评价［J］.企业经济，2014（8）：126-130.
［8］刘 勇，田 蕾，金 浩，等.中国四维城镇化协调性空间格局及演化研究［J］.天津大学学报（社会科学版），2015，17（6）：513-517.
［9］张 目，黄春燕，李 岩.基于相对熵和可变模糊集理论的战略性新兴产业企业信用评价［J］.数学的实践与认识，2014，44（13）：18-26.
［10］迟国泰，齐 菲，张 楠.基于最优组合赋权的城市生态评价模型及应用［J］.运筹与管理，2012，21（2）：183-191.
［11］刘耀彬，李仁东，宋学锋.中国城市化与生态环境耦合度分析［J］.自然资源学报，2005，20（1）：105-112.
［12］Moran P A P.Notes on continuous stochastic phenomena［J］.Biometrika，1950，37（1／2）：17-23.
中图分类号： F291.3
文献标志码： A
文章编号： 1002-1302（2019）13-0337-05
吴云青，王多多，郭 佳.基于ESDA的山东省四维城镇化协调发展空间差异［J］.江苏农业科学，2019，47（13）：337-341.
doi： 10.15889／j.issn.1002-1302.2019.13.078
收稿日期： 2018-03-13
基金项目： 国家社会科学基金（编号：18BZZ077）。
作者简介： 吴云青（1978—），男，山东临沂人，博士，教授，主要从事乡村发展与新型城镇化研究。E-mail：28640067＠qq.com。
通信作者： 王多多，硕士研究生，主要从事土地资源管理研究。E-mail：2316549156＠qq.com。
标签：城镇化论文; 协调发展论文; 空间差异论文; 探索性空间分析论文; 天津工业大学经济与管理学院论文;