运用空间自相关分析我国入境旅游增长空间格局.docx

1、运用空间自相关分析我国入境旅游增长空间格局运用空间自相关分析中国入境旅游增长空间格局宋鸿/陈晓玲【专题名称】旅游治理【专 题 号】F9【复印期号】2006年07期【原文出处】世界地理研究(沪)2006年1期第99106页【作者简介】宋鸿，湖北大学商学院，武汉430062/宋鸿，湖北大学旅游进展研究院，武汉430062/宋鸿/陈晓玲，武汉大学测绘遥感信息工程国家重点实验室，武汉430079宋鸿，副教授，博士生，研究方向为区域经济研究中的空间分析方法与技术。【内容提要】空间自相关分析是空间统计学的一个重要组成部分，在设定显著性水平下研究邻近位置属性（或现象）之间的相关性，是认识空间格局的有效手段。

2、中国省级水平上的入境旅游区域增长空间自相关分析显示，1996年至2004年期间，除1997年外，我国在其它时期均不存在显著的全局空间自相关，表明邻近省区入境旅游区域增长之间的关系在整体上既不综合表现为趋同，也不综合表现为趋异，入境旅游区域增长整体空间格局为随机格局；各个时期均有一个或少数几个省区与邻近省区的局部空间自相关显著，在局部区域呈现出集聚或离散的空间格局。这类局部区域的个数随时刻变化，但2001年以后趋于稳定，在没有对入境旅游产生重大冲击事件发生的正常时期为2个，即以上海为核心的共同增长集聚格局区域和以广东省为中心的“中高周低”的离散格局区域（0.05，双侧显著性检验）。【关 键 词】

3、空间自相关/空间格局/Morans I/入境旅游1入境旅游增长空间格局研究的意义入境旅游作为我国旅游业的重要组成部分，一直是一个备受关注的领域，关注的核心是促进入境旅游增长。我国入境旅游增长来源于各省、自治区、直辖市（以下简称省区）等各区域单元的增长，各省区入境旅游增长空间格局是总增长形成的背景。因此，认识我国入境旅游区域增长空间特征及其演变过程，关于制定有效的入境旅游进展战略有着重要的意义。国内已有学者从空间的角度对我国入境旅游增长进行了研究。周云波等（1999）1 以区域经济非均衡增长理论为基础对中国国际旅游业的区域增长进行了分析，指出19901997年期间中国国际旅游业区域增长呈现明显的

4、非均衡态势，呈梯形进展格局，即：东部最发达，中部次之，西部最差；陆林等（2005）2 通过对入境旅游外汇收入省际差异的分析，指出我国入境旅游外汇收入的空间差异较大，19902002年期间我国高于当年平均水平的省区都集中在东部沿海，省区不超过6个；同时，由于进展基础和进展速度不一致，沿海地区的广东、上海、福建、江苏、浙江和北京等东部发达省区与全国其它省区之间，尽管相对差异不断缩小，然而绝对差异逐年扩大。这些研究深化了对我国入境旅游进展空间格局及其随时刻变化的认识，然而，对各省区入境旅游进展与邻近省区之间的相互关系缺乏探讨。空间自相关分析是研究邻近位置的属性（或现象）相关性的空间统计学方法，上个世

5、纪中期开始在美国、英国等西方国家引起关注，随后在世界范围内得以快速进展。目前，空间自相关分析正成为定量研究自然、经济和社会领域内涉及空间关系的各类问题的重要方法和分析空间格局的有效手段。因此，对我国入境旅游增长空间格局的空间自相关分析，将有利于加深对我国入境旅游区域增长空间特征的认识。从查阅的文献来看（维普资讯中文科技期刊数据库：1989年现在），尚未见此类研究成果。2差不多概念与理论2.1空间自相关与空间格局空间自相关（spatial autocorrelation）是指一个变量的观测值之间因观测点在空间上邻近而形成的相关性（auto即为“自”）3。依据分析空间范围的大小，空间自相关可分为全

6、局空间自相关（Global Spatial Autocorrelation）和局部空间自相关（Local Spatial Autocorrelation）。 全局空间自相关指研究范围内邻近位置同一属性相关性的综合水平，局部空间自相关指研究范围内各空间位置与各自周围邻近位置的同一属性相关性；依据空间自相关的性质，空间自相关可分为正空间自相关（Positive Spatial Autocorrelation）、 负空间自相关（Negative Spatial Autocorrelation）和无空间自相关。按Griffth3 的观点， 正空间自相关指邻近位置的同一变量观测值之间呈正相关，负空间自相

7、关指邻近位置的同一变量观测值之间呈负相关。与传统统计学的变量相关性一样，正空间自相关或负空间自相关也有强弱之分，可通过空间自相关指数度量。空间自相关分析是认识空间格局的有效手段。空间格局是观测属性及其在空间上的相互关系。以空间邻近位置属性的相似性为依据，空间格局能够分为集聚的（clustered）、离散的（dispersed）和随机的（random）等三种类型。空间集聚格局指相似属性在空间上邻近，空间离散格局指被关注的位置与其邻近位置在属性上迥异，其它情况归为随机空间格局。显然，空间自相关与空间格局存在着对应的关系，正空间自相关对应于集聚格局、负空间自相关对应于离散格局，当不存在空间自相关时，

8、属性观测值呈随机分布。2.2空间自相关指数与空间格局空间自相关指数是定量测度空间自相关的统计量。由前述分析可知，空间自相关与空间格局存在对应的关系，因此，空间格局能够通过空间自相关指数进行定量分析。常用的空间自相关指数有Morans I、Geary Ratio C、General G-Statistics、LISA、Local General G-Statistics等等。这些指数的差不多思想类似，然而计算方法不同。下面要紧对本文将用于实际计算的Morans I和LISA予以介绍。欲了解其它指数，请参考相应文献4,5。2.2.1全局Morans I与空间格局全局Morans I（Global

9、Morans I）是最常用的全局空间自相关指数，计算公式如下6：式中：n为区域单元个数x,i为区域单元i的属性观测值，i1，2，nW,ij为空间权重系数矩阵，表示各区域单元空间邻近关系区域单元间邻近关系通常按邻接标准或距离标准进行定义，空间邻近关系的判定与空间权重系数矩阵的表达之间的关系见表1。表1区域单元空间邻近关系与空间权重系数矩阵 权重系数邻接标准距离标准区域i与j邻接区域i与j不邻接d,ijd,0d,ij d,0W,ij1 0 1 0注：d,ij指区域i与j之间的距离，d,0为设定距离（能够是空间距离、时刻距离等）。用邻接标准定义区域单元间的邻近关系时，邻接又分为几种不同的具体情况。依

10、照邻接方式，可分为Rook（共边为邻接）、Bishop（共点为邻接）和Queen（共边或共点均为邻接）；依照是否是直接邻接，可分为一阶邻接（first order spatial contiguity，即直接邻接）、二阶邻接（second order spatial contiguity，通过一阶邻接区域单元与其他区域单元形成的邻接）、高阶邻接（higher order spatial contiguity，二阶邻接的推广）。在实际的运用中，多采纳Rook方式、一阶邻接定义邻近关系，构建空间权重系数矩阵。由Global Morans I的大小可对空间格局进行定量分析。当Global Moran

11、s I接近1时，研究区域在整体上存在强的正空间自相关，观测属性呈集聚空间格局；反之，Global Morans I接近1时，研究区域在整体上存在强的负空间自相关，观测属性呈离散空间格局。当观测属性不存在空间自相关时，Morans I1（n1），观测属性在空间上呈随机分布。由于Global Morans I的计算使用的是样本数据，因此，当Global Morans I1（n1）时，Global Morans I与1（n1）之间的差异存在两种可能：一是显著不同；二是不存在显著不同，Global Morans I与1（n1）之间的差异是抽样的随机性所致。因此，基于Global Morans I分析观

12、测变量的空间格局时，需要进行显著性检验，以便在一定概率下保证推断结论的正确性。对Global Morans I的显著性检验通常采纳双侧检验。检验的零假设为H,0：Global Morans I1（n1），即观测变量呈随机空间分布；备择假设H,1：Global Morans I1（n1），观测变量呈非随机空间分布；检验统计量用Z统计量：式中：I为待检验的由样本数据计算得到的Global Morans IE（I）为Global Morans I的期望值，其值为1（n1）（n为样本容量）Var（I）为Global Morans I的方差，计算方法与抽样假设有关。抽样假设有两种：正态抽样假设和随机抽样

13、假设。更详细的内容及计算方法请参考相应文献6。最后依照Z值大小，在设定显著性水平下做出同意或拒绝零假设的推断。取0.05，则当Z1.96或Z1.96时，拒绝零假设，观测变量的空间自相关显著，观测属性在空间上呈离散格局（Z1.96）或集聚格局（Z1.96）；反之，则同意零假设，观测变量在空间上呈随机分布。2.2.2LISA与空间格局Global Morans I是对研究区域空间自相关的综合度量，是各区域单元与相邻单元间空间自相关性质与水平的综合反映。显然，从研究区域内部来看，各局部区域的空间自相关完全一致的情况是专门少见的，常常是存在着不同水平与性质的空间自相关，这种现象称为空间异质性（spat

14、ial heterogeneity）。揭示空间自相关的空间异质性可用LISA（Local Indicators of Spatial Association，简称LISA）。LISA是一组指数的总称，如Local Morans I、Local Geary Ratio C、Local Gamma等等。下面对本文研究中使用的Local Morans I进行简单介绍。区域单元i的Local Morans I计算该单元与邻近单元间的空间自相关，计算公式如下6：式中：为区域单元i属性观测值的标准化值。W,ijW，一般为行标准化空间权重系数矩阵，即。现在，Global Morans I当I,i为正值且较大

15、时，区域单元i与相邻单元的观测属性存在较强的正空间自相关，呈局部空间聚拢（高值聚拢或低值聚拢）；反之，存在较强的负空间自相关，讲明区域单元i的属性值相关于相邻单元为高值或者低值离群点（即所谓的“热点”或者“冷点”）。与Global Morans I一样，用基于样本数据的Local Morans I分析局部空间格局时，需对Local Morans I进行显著性检验，检验方法同Global Morans I，详见文献6。通过显著性检验得到每个Ii的检验统计量值，并据此做出同意或拒绝零假设的结论。3指标选取与数据来源以省、自治区、直辖市为研究区域单元，以入境旅游外汇收入逐期增加值为反映各省区入境旅游

16、区域增长的指标。数据来源于中国统计年鉴（2005），时刻期限为1995年至2004年共10年。重庆直辖市成立于1997年，中国统计年鉴上只有1997年及以后的数据。通过计算发觉，重庆直辖后与四川省入境旅游外汇收入之比相对稳定，因此，出于延长时刻序列的考虑，依据1997年至2004年期间重庆直辖市与四川省入境旅游外汇收入之比的平均值，将四川省1995年和1996年的入境旅游外汇收入一分为二，对重庆直辖市成立之前的1995年和1996年，在研究中也按31个省区进行分析。4研究结果与讨论研究中采纳Rook方式、一阶邻接定义空间邻近关系，使用行标准化空间权重系数矩阵，显著性水平0.05，显著性检验类型

17、为双侧检验，计算任务用文献6 提供的程序完成。4.1我国入境旅游区域增长的全局空间自相关分析1996年至2004年期间，我国省级水平上的入境旅游外汇收入逐期增加值的Global Morans I及显著性检验统计量值（随机抽样假设）计算结果见表2。表中结果显示，尽管各年Global Morans I均不等于E（I），或大或小，然而，除1997年外，其它各年与E（I）均不存在显著差异，都同意零假设（HO：我国31个大陆省区入境旅游外汇收入逐期增加值在空间上呈随机分布），不存在显著的全局空间自相关，表明我国邻近省区入境旅游区域增长既不综合表现为趋同，也不综合表现为趋异，入境旅游区域增长整体格局为随机

18、空间格局。表2入境旅游外汇收入逐期增加值的Global Morans I及Z值 时期(年)Morans IE(I) VAR,R(I)Z,R(I)1996 0.035 -0.0333 0.0110.641997 0.202 -0.0333 0.0132.081998 -0.146-0.0333 0.014-0.961999 0.023 -0.0333 0.0090.582000 0.002 -0.0333 0.0070.432001 -0.019-0.0333 0.0120.132002 0.018 -0.0333 0.0121.302003 -0.071-0.0333 0.009-0.3820

19、04 0.043 -0.0333 0.0130.674.2我国入境旅游区域增长的局部空间自相关分析由Global Morans I可知，在1996年至2004年期间，除1997年外，我国大陆各省区入境旅游外汇收入逐期增加值整体上表现出随机分布的特征。但这只是揭示了全局综合特征，还不能表明各省区与周边邻近省区之间的局部空间关系在各个年份差不多上如此。4.2.1Local Morans I的时空变化1996年至2004年期间，各省区入境旅游外汇收入逐期增加值的Local Morans I箱线图见图1。图中显示，Local Morans I随时刻发生变化，变化幅度存在显著省区差异。上海、江苏、海南、

20、北京、广东、福建、天津等省区明显较大，特不是上海市尤为突出；江西及其右边的18个省区则专门小。Local Morans I的变化幅度反映空间关系的稳定性，变化幅度小，表明空间关系随时刻变化小。图11996年至2004年入境旅游外汇收入逐期增加值Local Morans I箱线图从图中能够看出，各省区Local Morans I的中位数也存在较大差异。上海市、海南省、广东省、福建省等省区Local Morans I的中位数绝对值较大，其余27个省区都专门小。这27个省区Local Morans I的中位数绝对值均小于0.5、四分位距IQR均小于0.7（除江苏外的26个省区，IQR均小于0.4），

21、各省区的Local Morans I与期望值E（I,i）之间的差值专门小且稳定（当使用行标准化空间权重系数矩阵时，所有区域单元i的Local Morans I的期望值相同，E（I,i）E （I）1（n1）0.033）。表明在1996年至2004年期间的多数时期内，大部分省区入境旅游外汇收入逐期增加值与相邻省区构成的空间格局趋向于随机格局；同时，Local Morans I中位数绝对值的大小与Local Morans I年际变化幅度有着明显的联系。绝对值大的省区，Local Morans I的年际变化大，反之则小，讲明在我国内部各局部区域形成的随机空间格局具有专门强的稳定性。4.2.2局部区域的空间格局分析结合各省区Local Morans I显著性检验结果进一步分析表明，1996年至2004年期间，我国入境旅游外汇收入逐期增加值每年都在少数局部区域存在显著的局部空间自相关，呈现出或聚拢或离散的空间格局，同时，这些局部区域及其空间格局的类型存在着年际变化（表3）。表3入境旅游外汇收入逐期增加值局部空间自相关性显著的省区及其Local Morans I