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

徐海峰

（北京勞動(dòng)保障職業(yè)學(xué)院，北京 100029）

一、引言

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

二、實(shí)證研究設(shè)計(jì)

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

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

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

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

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

采用目標(biāo)值法對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理。標(biāo)準(zhǔn)化公式為：

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

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

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

3.向量自回歸模型設(shè)置

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

三、基于VAR模型實(shí)證分析

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

1.格蘭杰因果檢驗(yàn)

（1）單位根檢驗(yàn)

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

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

（2）協(xié)整檢驗(yàn)

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

表3 VAR模型p值選擇表

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

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

（3）誤差修正模型

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

（4）格蘭杰因果檢驗(yàn)

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

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

2.脈沖響應(yīng)分析

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

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

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

3.方差分解

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

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

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

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

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

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

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