基于向量自回歸模型北京城鎮(zhèn)化發(fā)展動力實證研究

徐海峰

（北京勞動保障職業(yè)學院，北京 100029）

一、引言

隨著經(jīng)濟形態(tài)向服務型經(jīng)濟轉(zhuǎn)型，服務業(yè)因其產(chǎn)業(yè)關聯(lián)性強加速推進產(chǎn)業(yè)集聚與人口集中，促進城鎮(zhèn)化發(fā)展。流通業(yè)和旅游業(yè)作為現(xiàn)代服務業(yè)的主導產(chǎn)業(yè)，其與城鎮(zhèn)化的互動作用日益顯著。流通業(yè)作為先導性產(chǎn)業(yè)和基礎性產(chǎn)業(yè)，因其產(chǎn)業(yè)關聯(lián)度大，聯(lián)動產(chǎn)業(yè)升級，為城鎮(zhèn)化提供強有力的產(chǎn)業(yè)環(huán)境，有效實現(xiàn)產(chǎn)城融合發(fā)展（鄭勇軍等，2014；顏雙波，2016），并因其就業(yè)吸納能力和帶動能力強，成為拉動城鎮(zhèn)就業(yè)的有效產(chǎn)業(yè)因素（王曉東、謝莉娟，2010）。旅游業(yè)作為現(xiàn)代服務業(yè)的重要組成部分，其與區(qū)域經(jīng)濟協(xié)調(diào)發(fā)展能夠促進城鎮(zhèn)化（鐘家雨等，2014），城鎮(zhèn)通過其公共產(chǎn)業(yè)及資源等旅游化促進區(qū)域經(jīng)濟發(fā)展，有助于提升城鎮(zhèn)化水平（林峰等，2013），旅游業(yè)與城鎮(zhèn)化之間存在共同演化關系（黃劍鋒等，2021）。本文基于產(chǎn)城融合發(fā)展的視角，以流通業(yè)和旅游業(yè)為例構(gòu)建產(chǎn)城多系統(tǒng)向量自回歸模型（VAR）實證研究城鎮(zhèn)化發(fā)展動力問題，探究產(chǎn)業(yè)發(fā)展促進城鎮(zhèn)化的演進規(guī)律。

二、實證研究設計

1.變量指標體系構(gòu)建

城鎮(zhèn)化水平通常用人口城鎮(zhèn)化率來表示，即城鎮(zhèn)常住人口占常住總?cè)丝诘谋戎兀洖閁RB。流通業(yè)發(fā)展水平表示為CIR，借鑒相關文獻（王曉東、謝莉娟，2010），本文分別從流通業(yè)產(chǎn)值、產(chǎn)業(yè)活動單位數(shù)和流通業(yè)就業(yè)量這4個方面構(gòu)建流通業(yè)發(fā)展水平指標體系。旅游業(yè)發(fā)展水平表示為TOU，結(jié)合相關文獻（周蕾、王沖，2017），本文分別從旅游收入與旅游人數(shù)兩方面構(gòu)建旅游業(yè)發(fā)展水平指標體系（見表1）。指標數(shù)據(jù)采用1999年～2016年北京常住人口數(shù)據(jù)和流通業(yè)、旅游業(yè)等產(chǎn)業(yè)數(shù)據(jù)，來源于《北京統(tǒng)計年鑒》。

表1 城鎮(zhèn)化與流通業(yè)、旅游業(yè)發(fā)展水平評價指標體系

2.產(chǎn)業(yè)發(fā)展水平測算

（1）產(chǎn)業(yè)數(shù)據(jù)標準化處理

采用目標值法對數(shù)據(jù)進行標準化處理。標準化公式為：

對于某產(chǎn)業(yè)，xij和Xij分別表示第i個樣本的第j個指標的標準化值和原始值。Mj為第j個指標的目標值，參考文獻（李江蘇等，2014），本文選擇第j個指標的最大值作為目標值。

（2）產(chǎn)業(yè)發(fā)展水平測算

產(chǎn)業(yè)發(fā)展水平采取線性加權(quán)的方法進行計算（張勇等，2013），如式（2）所示，其中，U表示該產(chǎn)業(yè)的發(fā)展水平，wj表示該產(chǎn)業(yè)第j個指標的權(quán)重。采用熵值賦權(quán)法確定指標權(quán)重。

3.向量自回歸模型設置

城鎮(zhèn)化與流通業(yè)、旅游業(yè)之間相互作用相互影響，存在著互動作用關系，采用向量自回歸模型（VAR）對這種關系進行實證分析更加具有科學性。建立城鎮(zhèn)化與流通業(yè)、旅游業(yè)的VAR模型為：

三、基于VAR模型實證分析

運用MATLAB，根據(jù)式（2）計算得到歷年北京流通業(yè)發(fā)展水平指數(shù)（CIR）、旅游業(yè)發(fā)展水平指數(shù)（TOU），限于篇幅，計算結(jié)果略過。基于VAR模型采用Eviews 8.0實證分析。

1.格蘭杰因果檢驗

（1）單位根檢驗

采用ADF單位根檢驗對URB、CIR、TOU等變量進行平穩(wěn)性檢驗。表2顯示，3個變量是非平穩(wěn)序列，但是3個變量的一階差分序列都是平穩(wěn)的，且均為I（1）序列，即同階單整。

表2 ADF單位根平穩(wěn)性檢驗

（2）協(xié)整檢驗

采用Johansen檢驗法，用以確定變量URB、CIR、TOU之間有無協(xié)整關系。根據(jù)表3，VAR（p）模型的滯后階數(shù)p取值為3。

表3 VAR模型p值選擇表

基于VAR模型的協(xié)整檢驗的滯后階數(shù)為p-1，因此選擇滯后階數(shù)為2，并對協(xié)整方程的形式加以確定，通過聯(lián)合檢驗，選擇序列有線性確定性趨勢且協(xié)整方程僅有截距的檢驗模型。根據(jù)表4可知，3個變量之間存在協(xié)整關系。標準化的協(xié)整方程如式（6）所示。由此可知城鎮(zhèn)化與流通業(yè)、旅游業(yè)都是正相關關系，表明流通業(yè)和旅游業(yè)是城鎮(zhèn)化發(fā)展的長期原因。

表4 Johansen協(xié)整檢驗結(jié)果

（3）誤差修正模型

城鎮(zhèn)化與流通業(yè)、旅游業(yè)之間存在長期均衡關系，但是短期波動會偏離長期均衡，建立誤差修正模型（VEC）研究它們之間的短期動態(tài)關系。由式（6）得到協(xié)整方程殘差序列ut，令誤差修正項ECMt=ut，對模型進行系數(shù)估計后得到誤差修正模型如式（7）所示。根據(jù)估計結(jié)果，ΔCIR和ΔTOU的估計系數(shù)為正，說明流通業(yè)和旅游業(yè)發(fā)展都對城鎮(zhèn)化具有短期促進作用。ECMt-1項的系數(shù)為負，表明短期波動偏離長期均衡狀態(tài)時會受到反作用力以維持長期均衡。

（4）格蘭杰因果檢驗

根據(jù)表5可知，城鎮(zhèn)化與流通業(yè)之間互為格蘭杰原因；流通業(yè)是旅游業(yè)的格蘭杰原因，但是旅游業(yè)不是流通業(yè)的格蘭杰原因；旅游業(yè)是城鎮(zhèn)化的格蘭杰原因，但是城鎮(zhèn)化不是旅游業(yè)的格蘭杰原因。

表5 格蘭杰因果檢驗結(jié)果

2.脈沖響應分析

脈沖響應分析能夠反映城鎮(zhèn)化與流通業(yè)、旅游業(yè)VAR模型內(nèi)部變量間相互作用的結(jié)果，從而揭示三者之間的動態(tài)作用過程。AR根檢驗顯示該VAR模型是穩(wěn)定的。

在基期分別設置城鎮(zhèn)化、流通業(yè)和旅游業(yè)為1個單位的脈沖，考察沖擊作用的期限設置為10年，通過考察各變量對沖擊的動態(tài)反映，具體分析沖擊對其產(chǎn)生的影響。脈沖響應結(jié)果如圖1所示。當城鎮(zhèn)化發(fā)生正向沖擊后，流通業(yè)和旅游業(yè)發(fā)展會發(fā)生響應：流通業(yè)表現(xiàn)為正向波動，波動幅度先下降，到第2期下降到最低點，然后再上升，到第5期正向波動趨于平穩(wěn)，即城鎮(zhèn)化與流通業(yè)發(fā)展趨于平穩(wěn)發(fā)展態(tài)勢；旅游業(yè)也表現(xiàn)為正向波動，波動幅度先下降再緩慢上升，到第7期趨于平穩(wěn)，即城鎮(zhèn)化與旅游業(yè)發(fā)展趨于平穩(wěn)發(fā)展態(tài)勢。當流通業(yè)發(fā)生正向沖擊后，城鎮(zhèn)化與旅游業(yè)發(fā)展會發(fā)生響應：城鎮(zhèn)化表現(xiàn)為正向波動，波動幅度緩慢下降，到第4期趨于平穩(wěn)；旅游業(yè)表現(xiàn)為正向上升波動，到第3期達到最大值，然后趨于穩(wěn)定。當旅游業(yè)發(fā)展正向沖擊后，城鎮(zhèn)化表現(xiàn)為先正向上升波動，到第2期達到最大值，然后正向下降波動，到第3期變?yōu)樨撓虿▌樱缓缶徛仙▌樱诘?期波動基本消失；流通業(yè)表現(xiàn)為先負向下降波動，到第2期達到最小值，然后負向上升波動，到第4期波動基本消失。

圖1 城鎮(zhèn)化與流通業(yè)、旅游業(yè)的脈沖響應

3.方差分解

利用Cholesky分解法進行方差分解，進一步解釋城鎮(zhèn)化與流通業(yè)、旅游業(yè)的相互作用程度。限于篇幅，此處略過變量預測誤差的方差分解結(jié)果。在城鎮(zhèn)化過程中，自身作用對其發(fā)展的貢獻率在第1期達到100%，由于流通業(yè)和旅游業(yè)的發(fā)展，城鎮(zhèn)化自身作用呈現(xiàn)先大幅下降后小幅上升的趨勢，到第15期貢獻率變?yōu)?3.16%，與此同時，流通業(yè)發(fā)展、旅游業(yè)發(fā)展對城鎮(zhèn)化作用的貢獻率都呈現(xiàn)先大幅上升后小幅下降的趨勢，由第1期的0%分別變?yōu)榈?5期的13.15%和3.69%，三者貢獻率的變化幅度越到后期越小，表明流通業(yè)發(fā)展和旅游業(yè)發(fā)展對城鎮(zhèn)化的作用逐漸趨于穩(wěn)定的平衡狀態(tài)，互動作用進入良性循環(huán)，城鎮(zhèn)化發(fā)展到新階段。

在流通業(yè)發(fā)展過程中，城鎮(zhèn)化對流通業(yè)發(fā)展起到了極大的推動作用，貢獻率由第1期的25.16%增長到第15期的74.50%。旅游業(yè)發(fā)展對流通業(yè)發(fā)展的貢獻率在第15期穩(wěn)定在3.31%。流通業(yè)發(fā)展對自身作用的貢獻率由第1期的74.84%下降到第15期的22.19%。在后期變化的幅度越來越小，即逐漸趨于平穩(wěn)，從第4期開始流通業(yè)發(fā)展的方差分解出現(xiàn)拐點，即城鎮(zhèn)化對流通業(yè)發(fā)展的貢獻率超過了流通業(yè)發(fā)展對于自身的貢獻。

在旅游業(yè)發(fā)展過程中，城鎮(zhèn)化對流通業(yè)發(fā)展起到了極大的推動作用，貢獻率由第1期的24.95%增長到第15期的72.72%。流通業(yè)發(fā)展對旅游業(yè)發(fā)展的貢獻率由1.39%提高到20.19%。旅游業(yè)發(fā)展對自身作用的貢獻率由第1期的73.66%下降到第15期的7.09%。在后期變化的幅度越來越小，即逐漸趨于平穩(wěn)，從第5期開始旅游業(yè)發(fā)展的方差分解出現(xiàn)拐點，即城鎮(zhèn)化、流通業(yè)發(fā)展對旅游業(yè)發(fā)展的貢獻率都超過了旅游業(yè)對于自身的貢獻。

四、研究結(jié)論與政策啟示

現(xiàn)代服務業(yè)促進城鎮(zhèn)化發(fā)展的作用日益凸顯，流通業(yè)和旅游業(yè)作為現(xiàn)代服務業(yè)的主導產(chǎn)業(yè)，其與城鎮(zhèn)化的動態(tài)互動作用顯著。本文以流通業(yè)和旅游業(yè)為例，建立產(chǎn)城多系統(tǒng)VAR模型分析產(chǎn)業(yè)發(fā)展對城鎮(zhèn)化發(fā)展的推動作用。實證分析結(jié)果表明：

（1）城鎮(zhèn)化與流通業(yè)、旅游業(yè)之間存在著長期穩(wěn)定的均衡關系，流通業(yè)和旅游業(yè)是城鎮(zhèn)化發(fā)展的長期原因；（2）三者之間也存在短期波動，短期波動時系統(tǒng)會受到反作用力被調(diào)整回長期均衡狀態(tài)。

政策啟示如下。第一，在城鎮(zhèn)化發(fā)展規(guī)劃中，要更加重視產(chǎn)業(yè)發(fā)展對城鎮(zhèn)化的支撐作用，促進流通業(yè)、旅游業(yè)等現(xiàn)代服務業(yè)高質(zhì)量發(fā)展，為高質(zhì)量城鎮(zhèn)化提供產(chǎn)業(yè)支撐。第二，充分發(fā)揮政府在產(chǎn)城融合發(fā)展中的引導作用。政府在推動城鎮(zhèn)化發(fā)展過程中，要積極推進產(chǎn)業(yè)結(jié)構(gòu)轉(zhuǎn)型升級，促進城鎮(zhèn)化高質(zhì)量發(fā)展。