1. 引言
纵观当今世界,旅游业已经成长为世界上规模最大、种类最多的经济产业。旅游业在国民经济发展中的重要性已经得到了许多学者的讨论。在许多国家,它已经成为经济收入的主要来源,就业,企业成长和发展基础设施 [1] 。在近十年以来,云南省旅游总收入从最初2002~2003年7.7%的增长率到2016~2017年46%跳跃式增长,逐步成为拉动云南GDP增长的支柱产业。本文选取2000年至2017年云南省旅游业总收入和GDP作为抽样样本,再运用单位根以及协整检验等多种计量方法进行实证分析。结论为——发展旅游业的确能够拉动经济增长,且两者之间存在长期协整关系。我们计算了云南省旅游行业发展对云南省GDP拉动的贡献。结论表明,2000~2008年间云南旅游业对GDP的贡献呈现出稳定的势头。从2009年开始云南旅游业对GDP增长贡献率为13.14%,贡献率为20.8%;14年贡献率达41.74%。于此同时,在构建云南省旅游行业总收入和GDP的长期动态均衡模型后,发现云南省旅游产业增长率每增加或者降低1%,云南省GDP增加或降低0.2311755%。从长远来看,云南省的经济增长对旅游行业发展的依赖性较低。在短期内,云南省的经济增长对旅游业的发展高度依赖。
2. 文献综述
二十一世纪以来,特别是近年以来,随着旅游业在经济发展中重要性的日益增加,对旅游发展与经济增长关系的研究逐渐引起了全球学者的广泛关注。Gunduz和Hatemi-J通过应用杠杆辅助程序因果方法(Leveraged bootstrap causality)研究发现,旅游导致经济增长的假设在土耳其成立 [2] 。Balaguer和Cantavella Jorda通过构建实际GDP、实际国际旅游收入和实际有效利率等三个变量。西班牙旅游业与经济增长之间的关系是随着1975~1997年季度数据指标的讨论。结果表明,旅游业与经济增长之间存在协整关系,以及旅游业对西班牙经济增长的积极影响,从而支持旅游主导经济增长的研究假设 [3] 。1975~2001年韩国旅游业发展对经济增长贡献的研究发现,旅游业导致经济增长的假说尚未建立,旅游业发展与电子商务之间没有长期均衡关系 [4] 。Ditsakis利用GDP的协整向量、实际有效汇率和国际旅游收入用类似的方法对希腊的试验。通过修正误差模型的格兰杰因果检验证证明了经济增长和旅游业发展之间较强的因果关系 [5] 。同过去两个变量或者多变量模型不同的是,Durbarry构建了包含5个变量(资本存量、实际白糖出口量、人力资本、实际制造品出口量和实际旅游收入)的ECM模型来讨论1952~1999年旅游对毛里求斯共和国经济发展的影响,研究结果支持旅游促进经济增长的观点 [6] 。Kim,et al.构建两个变量的修正误差模型,基于1971~2003年的季度数据讨论了中国台湾经济发展与旅游业之间的联系,结论表明,经济发展和旅游业发展互为因果 [7] 。
近些年,我国出现了为数不多的有关中国旅游与经济增长的长期动态均衡计量研究。国内这项研究需要改进:首先,中国相较于国外研究样本容量较少。在以往研究中时间维度较短,大约为10年,受到一定程度测量结果的准确性的限制。其次,多使用国际旅游收入作为解释变量,从而降低了旅游总收入的价值,一定程度上使数据的估计出现偏差。在本质上,进行旅游业研究可以揭示旅游业对云南GDP增长的意义。本文基于经典的计量经济学模型检测方法,云南省2000~2017年的时间序列数据来检验云南省GDP和旅游业发展的长期均衡关系。
本文主要的研究思路如下:首先,对旅游业收入与的GDP序列数据的单位根检验确定其时间序列的稳定性。在此基础上,用E-G协整模型检验变量之间的协整关系,以确定两者之间的长期均衡关系,其次,利用格兰杰因果模型验证旅游业收入以及GDP之间的因果关系。最后利用ECM模型进行修正。
3. 云南省旅游业与经济增长关系的实证研究
3.1. 计量方法及模型
3.1.1. 单位根检验
包含如下三种回归模型:
模型一:
(3-1)
模型二:
(3-2)
模型三:
(3-3)
其中,δt代表线性趋势函数,t代表常数。原假设为H0: m = 0,则原序列存在一个单位根;备择假设为H1:m < 0,否则不存在单位根。
3.1.2. 协整检验
E-G两步协整检测法:
第一步,建立回归方程。
(3-4)
模型估计的残差为
(3-5)
第二步,是检验残差序列等是否稳定,如果它是平稳的,它有可能确定双方之间的协整关系,否则两者之间不存在协整关系 [8] 。
3.1.3. 格兰杰因果关系检验
在经济学研究中,有时需要判断经济变量的因果关系。格兰杰因果关系是指,若前变量可以在一定程度上说明或者预测后者,前者加上前者的滞后变量,则可以使这种说明或预测能力更加有效,则X是Y的格兰杰原因,它的检验方程如下所示 [8] :
(3-6)
(3-7)
3.1.4. 误差修正模型(ECM)
对于生活中采集的时间序列数据,若直接用于建立模型,可能会出现一些意想不到的错误。这是因为,在实际经济生活中,我们所采集的经济变量常受其滞后变量或者其他变量的影响 [9] 。此外,同时也会受到某些不相关因素的影响。所以,大多数情况下不能明确揭示经济变量之间的关系。同时,收集到的未经过处理的原始数据通常不表现为序列平稳。如果将它们直接作为抽样样本,则最终结果必然与实际情况相差较远,如果将原始数据先差分然后建模,会丢失一些有价值的信息,结论的准确性和可靠性将会降低。鉴于此,建模前需要进行协整检验。当变量之间存在协整关系时,建立ECM模型。
Engel-Granger两步法
1) 进行协整回归:
(3-8)
得到残差序列:
(3-9)
2) 若变量之间确存在协整关系,则用最小二乘方法进行检验。
3.2. 实证分析
3.2.1. 数据的选择与处理
本文以云南省旅游总业收入(记为X)和云南省国民生产总值GDP(记为Y)作为刻画旅游业和经济增长的指标,从云南省历年的统计年鉴中收集2000~2017年的时间序列数据作为样本。为了抑制异方差性,将所有数据取自然对数,变量的对数形式能够反映变量的长期弹性。经过数据处理以后,旅游业的收入(X)和云南省的GDP(Y)分别表示为lnX和lnY。实证检验数据见附录。
3.2.2. 数据相关性分析
为了直观的表示lnX和lnY之间的关系,我们绘制了两个变量的散点图。从图中可以了解到,lnX和lnY表现为正相关关系。
从图1可以看出,变量lnX和lnY有着明显的正相关关系,即云南省旅游业收入的变化与GDP增长的方向是相同。然而,基于散点图分析只是最初步的分析,其深层的关系则需要更加严格和科学的检测。
3.2.3. 单位根检验
为了确保分析的科学性以及后面实验的有序的进行,首先,我们测试最初时间序列的平稳性。如表1中所示,经济变量lnX、lnY和DlnX、DlnY的ADF检验统计量分别为1.186758、−1.328076和−0.836767、−2.776368。经济变量lnX、lnY和lnX,lnY的一阶差分变量DlnX和DlnY的ADF检验统计量如表1所示

Figure 1. Scatter plot based on the time series data of lnX and lnY in Yunnan from 2000 to 2017
图1. 基于云南2000~2017年时间序列数据lnX和lnY的散点图

Table 1. ADF test results of variables lnX, lnY, DlnX, DlnY, DDlnX, DDlnY
表1. 变量lnX、lnY、DlnX、DlnY、DDlnX、DDlnY的ADF检验结果
注:DlnX和DlnY都是原始序列的一阶差分序列。DDlnX和DDlnY都是原始序列的二阶差分序列。(C,T,N)表示截距项,时间趋势,滞后项。滞后项n以赤池信息准则AIC最小值为标准。
都大于其对应ADF检验临界值。结论显示,通过这种测试,上述四个经济变量有单位根,即lnX、lnY、DlnX和DlnY不是平稳的。lnX,lnY二阶差分变量DDlnX和DDlnY的ADF检验统计量的值都小于其对应的ADF检验的临界值,说明时间序列DDlnX和DDlnY是平稳的,则时间序列LnX和LnY是二阶平稳的。
关于时间趋势项和截距项的选择与否可以通过图2、图3和图4来进行判断:
如图2所示,序列lnX和lnY的平均值均不等于0,这意味两者应该都有截距项;此外,lnX和lnY的时序图可以看出两组序列都有时间趋势项。从图3可以看出,变量LnX和lnY的一阶差分序列DlnX和DlnY也应该有截距项。序列DlnX有时间趋势项,DlnY则没有。同理,如图4所示,很明显看出二阶差分序列DDlnX和DDlnY既不存在截距项同时不存在没有时间趋势项。
3.2.4. E-G两步协整检验
为了测试序列lnX和lnY之间的协整关系,需要进行协整检验 [10] 。

Figure 2. Sequence diagram of scatter plots based on the time series data lnX and lnY of Yunnan 2000-2017
图2. 基于云南2000~2017年时间序列数据lnX和lnY的散点图的时序图

Figure 3. Sequence diagram of scatter plot based on the time series data DlnX and DlnY of Yunnan 2000-2017
图3. 基于云南2000~2017年时间序列数据DlnX和DlnY的散点图的时序图

Figure 4. Based on the time series of DlnX and DlnY time series data of Yunnan 2000-2017
图4. 基于云南2000~2017年时间序列数据DlnX和DlnY的散点图的时序图
1) 构建模型。假设上述序列模型为
(3-10)
2) 进行回归分析:
(3-11)
则残差为如下方程所示:
(3-12)
回归方程残差序列曲线图如图5所示。
对残差序列进行均值检验,其结果见表2,再结合图5和上述方法,能够判断出,此时间序列没有时间趋势项和截距项。
通过对残差序列做单位根检验,残差序列的ADF统计量为−4.223113 < −2.740613,故在1%显著水平下残差序列是平稳的(表3)。所以,序列lnX和lnY存在协整关系,即2000~2017年,云南省旅游业和经济发展的代理变量之间存在一种长期的动态均衡关系。
3.2.5. 格兰杰因果关系检验
根据上文内容中对云南省旅游总收入以及云南省国民经济总量的单位根检验结果,在进行下一步操作前,首先要对附录 [11] 中的数据进行差分处理,所得结果如表4所示。
经过处理后表5中二阶差分项的数据便是我们要进行格兰杰因果检验的数据项,根据上文中格兰杰因果检验方法的介绍,通过EViews9.0检验,得出如表5中的检验结果:

Figure 5. Regression equation residual timing diagram
图5. 回归方程残差时序图

Table 2. Mean test results of residual series
表2. 残差序列均值检验结果

Table 3. Unit root test results for residuals
表3. 残差的单位根检验结果
注:(C,T,N)表示截距项,时间趋势项和滞后阶数,滞后长度n是以AIC最小值为基准。

Table 4. Yunnan tourism total income X and Yunnan Province GDP difference and logarithmic processing
表4. 云南省旅游总收入X以及云南省GDP差分及对数处理
注:---表示缺省值。

Table 5. Granger causality test results
表5. 格兰杰因果检验结果
注:1. Null Hypothesis的意思是原假设,Lags是滞后值;2. Pro.意思是发生的概率;3. Pro.意思是发生的概率,F临界值是通过查F临界值分布表得到的在显著性为10%的情况下,F分布的临界值。
表5中在该计算过程中,我们可以看到如果选择5%显著水平,则假设“DDln(X) does not Granger Cause DDln(Y)”的p值大于0.05,且F分布临界值3.55047 < 4.1,假设被接受,旅游总收入增加不是经济增长的格兰杰原因;同时,假设“DDln(Y) does not Granger Cause DDln(X)”,的p值远远大于0.05,且F检验值0.55014小于临界值4.1,原假设被接受,即云南省经济增长不是旅游总收入增加的格兰杰原因,云南省旅游总收入与云南省经济增长之间不存在格兰杰因果关系。但是,在10%的显著水平下,原假设“DDln(X) does not Granger Cause DDln(Y)”的p值0.09220.1,且F捡验值3.55047大于F临界值2.92,原假设被拒绝,即云南省旅游总收入增长是经济增长的格兰杰原因;由原假设“DDln(Y) does not Granger Cause DDln(X)”的p值0.4772 > 0.1,且F检验值0.55014小于临界值2.92,原假设被接受,即云南省经济增长不是旅游总收入增长的格兰杰原因。
通过以上检验结果看出,在放宽显著水平的情况下,云南省旅游产业发展 对经济增长的贡献的影响是显著地,这也反应了云南省旅游产业发展与经济增长之间的现实情况。云南省旅游产业是一项综合性极强的产业,对其他产业的发展具有较强的关联效应,对云南省经济的提升具有积极促进作用。简单来说,云南省旅游产业的发展必然会引起与旅游产业相关的各个部门的投入,而这些投入又必然使得这些行业的规模不断扩大,取得行业上的发展和进步,从而使得云南省经济总量增长。
3.2.6. 误差修正模型
基于以上述检验的结果,先令
,然后对变量DlnX、DlnY以及ECMt-1进行回归,构建ECM模型 [12] :
(3-13)
通过回归可得数据如表6所示:
由表6可得误差修正模型:
(3-14)
由ECM模型中可以看出,在短期里,旅游业收入增大1%,GDP将会随之增大0.13%;在长期来看里,这个数字为0.66%。所以,对于云南省就其旅游业对其GDP增长效果而言,长期比短期更为显著。
3.3. 实证结果分析
1) 变量lnX与lnY之间的协整关系表明旅游业发展与云南省经济增长存在长期稳定的动态均衡关系。由(3-11)可知,前者增长1%,后者将会增长0.66%。由此可见,云南省旅游业发展能较为明显地带动云南省经济的发展。从(3-14)可以看出,造成lnY短期波动的因素由两个方面构成:第一是短期波动,即LnX短期波动所造成的;第二是长期动态均衡,即通过ECM长期调整达到均衡的调整过程。如果短期波动偶发偏离长期均衡状态,则会引起偏差但是这种偏差不会太久。跟随每一个周期(−0.351587)进行调整,并不断修正这个偏差直到它返回到平衡状态为止。
2) 从格兰杰因果检验可以看出,对于旅游业和经济增长,在5%的显著性水平下,旅游业发展是经济增长的格兰杰原因,但经济增长却不是云南旅游业发展的格兰杰原因,云南旅游业和经济发展存在单向因果关系。换句话说,在过去的18年里,旅游业能够较为明显地拉动云南经济增长,而云南经济增长未能明显地拉动旅游业收入的增长。
3) 从ECM模型中可以看出,在短期内,旅游业收入每增加1%,GDP将会提高0.13%,而对于(1)中我们通过计算得到这个数字是0.66%。因此,对于云南省而言,旅游业对其经济的增长长期相较于短期更为显著。与此同时,经济系统具有较强的自我调节功能,并能自发保持动态均衡。如果短期波动偶发偏离长期均衡状态,则会引起偏差但是这种偏差不会太久。跟随每一个周期(−0.351587)进行调整,并不断修正这个偏差直到它返回到平衡状态为止。
4) 经济变量对GDP的推动效率是衡量经济变量对于GDP推动效率的衡量指标。其计算公式为 [13] :
(4-15)
在(4-15)公式中,αi代表后者较之于前者的弹性系数,βi则表示这段时间中前者的总和相较于后者总和的比重,γi代表一段时间里前者对后者的拉动效率。
若γi < 1时,说明在这段时间内,这这个经济变量对于GDP增长推动贡献低于其所占GDP的比值,换句话说,它推动经济增长效率较低。

Table 6. Regression analysis results
表6. 回归分析结果

Table 7. Pull efficiency calculation table
表7. 拉动效率计算表
若γi > 1时,说明在这段时间内,这个经济变量对于GDP增长推动贡献大于其所占GDP的比值,换句话说,它推动经济增长较为显著。
由此能够得到2000~2017年的拉动效率δi,如表7所示:
由表7可知,γi > 1,则说明在一定时间内,云南省旅游业积极地推动当地经济发展,且较为显著。
4. 政策建议
1) 政府部门在讨论制定当地经济发展和旅游业的计划和政策时,要将经济增长与旅游发展动态平衡和增强前瞻性的发展和长远发展考虑到其中。采取对应措施,以免短期经济波动对于当地旅游行业的冲击,保障当地旅游业和经济增长朝着更好更快的方向发展。
2) 由格兰杰因果检验可知虽然不存在经济增长到云南省旅游业发展的格兰杰原因 [14] ,但这并不能掩盖旅游业对于云南省经济发展的支柱地位,随着云南省旅游业不断地发展旅游发展,云南经济发展也随之不断的变化,其所占云南省GDP的份额也会逐步增大。从而旅游业对于推动云南省经济增长的假设是成立的。随着休闲时代的到来,休闲经济已经成为云南省经济发展的支柱产业。旅游经济的发展是云南省经济发展的战略性结构调整的“指南针”。因此,政府应当将大力发展休闲宜居型,环境友好型,资源节约型的旅游经济发展作为经济发展转型的方式之一。促进云南省接待国内和国际旅游游客的持续、稳定增长以此促进云南经济更好更健康的发展。
3) 旅游支付能力主要是消费者的可支配收入决定,这就要求企业对员工的工资进行提高,但是考虑到云南省的实际情况以及当前我国经济发展的现状,增加工资必须与企业发展与经济发展相适应。但是,在增加工资之外,企业可以考虑对内部员工开设旅游专户的奖励,该项奖励可以作为优秀员工考核绩效奖励,也可以作为年终奖励的部分内容,通过该项奖励不仅可以提高员工的积极性,同时也降低了企业的费用支出。
4) 加强国际旅游吸引力 [15] 。云南省地理资源、生物资源丰富,在国际市场上具有不可替代性,但是由于受到出行距离的影响,云南省内的国际旅游大多以东南亚等邻近国家游客为主,且重复旅游的几率较小。因此,云南省应当在此方面加大建设,提高国际旅游吸引力,吸引其他国家,尤其是欧美国家到云南进行旅游消费,增加旅游消费需求的有效性。
附录

Table A. Statistics on time series data of Tourism Total Income and GDP in Yunnan Province, 2000-2017
附表A. 云南省2000~2017年旅游总收入及其GDP相关时间数据