基于VAR模型的山東省產(chǎn)業(yè)發(fā)展系統(tǒng)評(píng)價(jià)

范夢(mèng)甜 王慧

摘要以山東省為例，以創(chuàng)新、協(xié)調(diào)、綠色、開(kāi)放、共享的發(fā)展理念為指導(dǎo)，構(gòu)建產(chǎn)業(yè)發(fā)展多指標(biāo)評(píng)價(jià)體系，測(cè)算改革開(kāi)放以來(lái)山東省產(chǎn)業(yè)發(fā)展演變情況，運(yùn)用基于VAR模型的脈沖響應(yīng)及方差分解分析產(chǎn)業(yè)發(fā)展的內(nèi)在機(jī)制，結(jié)果表明：①山東省產(chǎn)業(yè)發(fā)展水平不斷提高，經(jīng)歷了初期震蕩緩慢增長(zhǎng)、中期高速增長(zhǎng)和后期持續(xù)增長(zhǎng).②產(chǎn)業(yè)發(fā)展系統(tǒng)結(jié)構(gòu)逐漸優(yōu)化，呈現(xiàn)出產(chǎn)業(yè)規(guī)模、結(jié)構(gòu)、效益、創(chuàng)新、協(xié)調(diào)及綠色為一體的產(chǎn)業(yè)綜合發(fā)展態(tài)勢(shì)，目前主要面臨提升效益和發(fā)展外向型產(chǎn)業(yè)的雙重壓力.③從沖擊響應(yīng)結(jié)果來(lái)看，產(chǎn)業(yè)規(guī)模及外向型產(chǎn)業(yè)對(duì)產(chǎn)業(yè)綜合發(fā)展表現(xiàn)為正向效應(yīng)；產(chǎn)業(yè)和諧程度對(duì)產(chǎn)業(yè)綜合發(fā)展表現(xiàn)為負(fù)向效應(yīng)；產(chǎn)業(yè)效益及產(chǎn)業(yè)綠色化對(duì)產(chǎn)業(yè)綜合發(fā)展起初表現(xiàn)為正向效應(yīng)，之后作用不穩(wěn)定；產(chǎn)業(yè)結(jié)構(gòu)和產(chǎn)業(yè)科技創(chuàng)新能力對(duì)產(chǎn)業(yè)綜合發(fā)展的正向效應(yīng)存在一定滯后期.④預(yù)測(cè)方差結(jié)果表明，子系統(tǒng)對(duì)產(chǎn)業(yè)綜合發(fā)展的影響從高到低分別為：產(chǎn)業(yè)發(fā)展和諧程度，產(chǎn)業(yè)發(fā)展外向度，產(chǎn)業(yè)規(guī)模，產(chǎn)業(yè)科技創(chuàng)新能力，產(chǎn)業(yè)結(jié)構(gòu)，產(chǎn)業(yè)綠色發(fā)展，產(chǎn)業(yè)效益.

關(guān)鍵詞區(qū)域經(jīng)濟(jì)學(xué)；產(chǎn)業(yè)發(fā)展；VAR模型；山東省

中圖分類(lèi)號(hào)F127；F124文獻(xiàn)標(biāo)識(shí)碼A

Evaluation of industrial development system

in Shandong Province based on VAR model

Mengtian Fan， Hui Wang

（Department of Geography and Tourism， Qufu Normal University，Rizhao， Shandong276826， China）

AbstractThe study took innovation， coordination， green， open， sharing concept as an instruction and measured the evolution of industrial development by building multi-index evaluation system and analysed the internal mechanism of industrial development by using the impulse response and variance decomposition based on VAR model in Shandong Province since the reform and opening up and it found： （1） the industrial development level is enhancing gradually and it experienced the shock and slow growth in early period， high speed growth in middle period and sustained growth in end period； （2）the structure of industrial development system is optimizing gradually which presents the situation of industrial scale， structure， efficiency， innovation， harmonious and green as a whole and it currently mainly faces the double pressure of gaining efficiency and developing export-oriented industry；（3） the impulse response results show that the industrial scale and the export-oriented industry have positive effect on industrial comprehensive development， the industrial harmonious degree has negative effect on industrial comprehensive development， the industrial efficiency and the industrial green level have positive effect in early period and unstable effect in end period on industrial comprehensive development and the industrial structure and technological innovation ability have lagging positive effect on industrial comprehensive development； （4）variance results show that the influence of subsystem to industrial comprehensive development from high to low is： the harmonious degree of industrial development， the extroverted degree of industrial development， industrial scale， industrial technology innovation ability， industrial structure， industrial green development and industrial efficiency.

Key wordsregional economics；industrial development； VAR model；Shandong Province

1引言

在集聚效應(yīng)和規(guī)模經(jīng)濟(jì)推動(dòng)下，區(qū)域產(chǎn)業(yè)發(fā)展的規(guī)模和水平逐步提高，產(chǎn)業(yè)發(fā)展為區(qū)域經(jīng)濟(jì)社會(huì)全面發(fā)展奠定了物質(zhì)基礎(chǔ).目前，中國(guó)經(jīng)濟(jì)發(fā)展進(jìn)入新常態(tài)，產(chǎn)業(yè)發(fā)展不再局限于單一的數(shù)量擴(kuò)張，更注重產(chǎn)業(yè)的轉(zhuǎn)型升級(jí)，以創(chuàng)新、協(xié)調(diào)、綠色、開(kāi)放、共享的發(fā)展理念引領(lǐng)產(chǎn)業(yè)發(fā)展.在此背景下，合理評(píng)價(jià)產(chǎn)業(yè)發(fā)展系統(tǒng)，深入探討其內(nèi)在機(jī)制，客觀(guān)認(rèn)識(shí)區(qū)域產(chǎn)業(yè)發(fā)展態(tài)勢(shì)，對(duì)于促進(jìn)產(chǎn)業(yè)轉(zhuǎn)型發(fā)展具有重要的理論與實(shí)踐意義.

學(xué)術(shù)界基于地理學(xué)視角的產(chǎn)業(yè)發(fā)展研究主要集中在三次產(chǎn)業(yè)發(fā)展、產(chǎn)業(yè)結(jié)構(gòu)演進(jìn)、產(chǎn)業(yè)集群發(fā)展及產(chǎn)業(yè)綠色化等方面.李佳校（2016）[1]使用洛陽(yáng)市近年來(lái)的第一產(chǎn)業(yè)產(chǎn)值數(shù)據(jù)應(yīng)用灰色系統(tǒng)模型建模方法，對(duì)洛陽(yáng)市第一產(chǎn)業(yè)產(chǎn)值進(jìn)行了分析.章奇（2016）[2]對(duì)第二產(chǎn)業(yè)結(jié)構(gòu)調(diào)整與科技金融聯(lián)動(dòng)發(fā)展的關(guān)系進(jìn)行了研究.謝敏（2015）[3]等人以浙江省為例，采用空間基尼系數(shù)、區(qū)位熵對(duì)第三產(chǎn)業(yè)集聚程度進(jìn)行測(cè)算.對(duì)三次產(chǎn)業(yè)發(fā)展的研究，還集中在農(nóng)業(yè)[4]、工業(yè)[5]、戰(zhàn)略性新興產(chǎn)業(yè)[6]、高新技術(shù)產(chǎn)業(yè)[7]和文化創(chuàng)意產(chǎn)業(yè)[8]方面.對(duì)產(chǎn)業(yè)結(jié)構(gòu)演進(jìn)的研究主要集中于產(chǎn)業(yè)結(jié)構(gòu)特征、產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化及產(chǎn)業(yè)結(jié)構(gòu)調(diào)整方面.楊家偉（2013）[9]探究了河南省產(chǎn)業(yè)結(jié)構(gòu)演進(jìn)及其機(jī)理，宋周鶯（2013）[10]、楊智峰（2014）[11]、高遠(yuǎn)東（2015）[12]等人分別以河南省、西部地區(qū)和中國(guó)為例，研究了產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化問(wèn)題.何平（2014）[13]對(duì)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整進(jìn)行分析，發(fā)現(xiàn)目前中國(guó)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整并沒(méi)有取得“對(duì)經(jīng)濟(jì)發(fā)展方式轉(zhuǎn)變具有促進(jìn)作用”的優(yōu)化效果.對(duì)產(chǎn)業(yè)集群發(fā)展的研究主要集中于產(chǎn)業(yè)集群內(nèi)在機(jī)制和產(chǎn)業(yè)集群協(xié)同發(fā)展分析.劉媛媛（2014）[14]等人運(yùn)用回歸分析方法實(shí)證研究中國(guó)資源型產(chǎn)業(yè)集群形成動(dòng)因.萬(wàn)幼清（2015）[15]等人對(duì)產(chǎn)業(yè)集群協(xié)同創(chuàng)新風(fēng)險(xiǎn)進(jìn)行研究.趙春明（2016）[16]、李文（2016）[17]等人對(duì)產(chǎn)業(yè)綠色發(fā)展進(jìn)行了相關(guān)研究.以往的研究集中于產(chǎn)業(yè)發(fā)展的某一視角，將產(chǎn)業(yè)發(fā)展作為系統(tǒng)的研究相對(duì)較少.當(dāng)前，區(qū)域產(chǎn)業(yè)發(fā)展面臨經(jīng)濟(jì)新常態(tài)下產(chǎn)業(yè)轉(zhuǎn)型升級(jí)和復(fù)雜的國(guó)際背景，必須從產(chǎn)業(yè)發(fā)展各子系統(tǒng)出發(fā)，在創(chuàng)新、協(xié)調(diào)、綠色、開(kāi)放、共享的發(fā)展理念指導(dǎo)下，通盤(pán)謀劃區(qū)域產(chǎn)業(yè)發(fā)展.借鑒相關(guān)研究成果，以山東省為例，構(gòu)建產(chǎn)業(yè)發(fā)展及其子系統(tǒng)指標(biāo)評(píng)價(jià)體系，運(yùn)用計(jì)量方法測(cè)算改革開(kāi)放以來(lái)山東省產(chǎn)業(yè)發(fā)展系統(tǒng)演變態(tài)勢(shì)，從子系統(tǒng)視角探討產(chǎn)業(yè)發(fā)展的深層動(dòng)因，以期為山東省產(chǎn)業(yè)轉(zhuǎn)型升級(jí)提供決策參考.

2指標(biāo)體系構(gòu)建與數(shù)據(jù)來(lái)源

2.1指標(biāo)體系構(gòu)建

依據(jù)典型性、科學(xué)性、可獲取性的原則，以創(chuàng)新、協(xié)調(diào)、綠色、開(kāi)放、共享的發(fā)展理念為指導(dǎo)，借鑒已有研究成果對(duì)產(chǎn)業(yè)發(fā)展內(nèi)涵的界定[18-21]，從產(chǎn)業(yè)規(guī)模、結(jié)構(gòu)、效益、和諧程度、科技創(chuàng)新能力、綠色發(fā)展水平和外向度7個(gè)方面構(gòu)建產(chǎn)業(yè)綜合發(fā)展評(píng)價(jià)指標(biāo)體系.

多采用人均國(guó)民生產(chǎn)總值、地均二三產(chǎn)業(yè)增加值等指標(biāo)反映產(chǎn)業(yè)規(guī)模.產(chǎn)業(yè)結(jié)構(gòu)直接影響產(chǎn)業(yè)總體發(fā)展水平，其中包括制造業(yè)增加值占規(guī)模以上GDP比重以及二三產(chǎn)業(yè)就業(yè)比重等.采用公共財(cái)政預(yù)算收入占GDP比重和規(guī)模以上工業(yè)產(chǎn)值利稅率等度量區(qū)域產(chǎn)業(yè)效益.采用二三產(chǎn)業(yè)產(chǎn)值結(jié)構(gòu)與就業(yè)結(jié)構(gòu)偏離度和農(nóng)村人均純收入與城鎮(zhèn)居民人均可支配收入比率綜合反映區(qū)域產(chǎn)業(yè)間的協(xié)調(diào)態(tài)勢(shì).創(chuàng)新是產(chǎn)業(yè)發(fā)展和城鎮(zhèn)化的原始動(dòng)力，采用每萬(wàn)人在校大學(xué)生數(shù)和重要科技成果數(shù)量衡量區(qū)域產(chǎn)業(yè)創(chuàng)新能力.綠色化是全球可持續(xù)發(fā)展框架下對(duì)現(xiàn)代產(chǎn)業(yè)發(fā)展的要求，采用工業(yè)固體廢棄物綜合利用率和能源消費(fèi)彈性系數(shù)來(lái)反映產(chǎn)業(yè)發(fā)展的綠色化水平.改革開(kāi)放以來(lái)，中國(guó)依托外向型產(chǎn)業(yè)發(fā)展促進(jìn)區(qū)域產(chǎn)業(yè)發(fā)展，采用規(guī)模以上工業(yè)中外商和港澳臺(tái)投資比重和對(duì)外貿(mào)易依存度兩項(xiàng)指標(biāo)測(cè)算產(chǎn)業(yè)發(fā)展外向度.內(nèi)容層各指標(biāo)均包含兩項(xiàng)具體指標(biāo)，共14項(xiàng)指標(biāo)，構(gòu)建了產(chǎn)業(yè)發(fā)展的指標(biāo)評(píng)價(jià)體系（表1）.

2.2數(shù)據(jù)來(lái)源

數(shù)據(jù)來(lái)源于《輝煌山東60年資料匯編》、《山東省統(tǒng)計(jì)年鑒》（1979～2016）及歷年山東省國(guó)民經(jīng)濟(jì)統(tǒng)計(jì)公報(bào).部分指標(biāo)個(gè)別年份數(shù)據(jù)缺失，采取插值處理.產(chǎn)業(yè)發(fā)展指標(biāo)中人均GDP、二三產(chǎn)業(yè)增加值指標(biāo)按可比價(jià)格進(jìn)行核算.為消除量綱影響，對(duì)統(tǒng)計(jì)指標(biāo)數(shù)據(jù)進(jìn)行極差標(biāo)準(zhǔn)化處理.

3產(chǎn)業(yè)綜合發(fā)展水平測(cè)度

3.1研究方法

為消除主觀(guān)因素影響，采用熵值法確定產(chǎn)業(yè)綜合發(fā)展各指標(biāo)權(quán)重.熵值法是根據(jù)各指標(biāo)信息載量的大小來(lái)確定指標(biāo)權(quán)重的客觀(guān)賦權(quán)法，其基本原理是對(duì)m個(gè)評(píng)價(jià)方案和n項(xiàng)評(píng)價(jià)指標(biāo)形成的數(shù)據(jù)矩陣X=xijn×m，某項(xiàng)指標(biāo)xj指標(biāo)值xij的離散程度越大，則該指標(biāo)在評(píng)價(jià)體系中發(fā)揮的作用越大，其權(quán)重越大；反之，該指標(biāo)的權(quán)重越小[22].熵值法計(jì)算公式為：

Wj=gj/∑mj=1gj（1≤j≤m），（1）

式中，wj為第j項(xiàng)指標(biāo)的權(quán)重，gj為第j項(xiàng)指標(biāo)的差異系數(shù)，gj=1-ej，ej為第j項(xiàng)指標(biāo)熵值，ej=-[∑ni=1pijln （pij）]/ln （n），pij為第j項(xiàng)指標(biāo)下第i年份占該指標(biāo)比重，pij=xij/∑ni=1xij，（i=1，2，…，n，j=1，2，…，m），xij為標(biāo)準(zhǔn)化之后的第j項(xiàng)指標(biāo)第i年份數(shù)值.在此基礎(chǔ)上計(jì)算出產(chǎn)業(yè)綜合發(fā)展指數(shù)得分，其計(jì)算公式為：

si=∑mj=1wjpij，（i=1，2，…，n）.（2）

3.2產(chǎn)業(yè)發(fā)展指標(biāo)權(quán)重

根據(jù)熵值法計(jì)算出山東省產(chǎn)業(yè)綜合發(fā)展各指標(biāo)權(quán)重（表2）.

3.3產(chǎn)業(yè)綜合發(fā)展水平評(píng)價(jià)

根據(jù)表2及公式（2），結(jié)合統(tǒng)計(jì)數(shù)據(jù)，計(jì)算出1978～2015年山東省產(chǎn)業(yè)發(fā)展水平綜合指數(shù)，并繪制折線(xiàn)圖（圖1）.

注：圖中虛線(xiàn)為協(xié)調(diào)度擬合曲線(xiàn)

由圖1，研究期內(nèi)山東省產(chǎn)業(yè)綜合發(fā)展水平不斷提高.產(chǎn)業(yè)發(fā)展綜合指數(shù)由研究期初的0.16上升至研究期末的0.66，年平均增長(zhǎng)率為1.36%.擬合曲線(xiàn)呈現(xiàn)加速上揚(yáng)趨勢(shì).產(chǎn)業(yè)綜合發(fā)展在時(shí)間上呈現(xiàn)出典型的階段性特點(diǎn).1978～1992年為震蕩緩慢增長(zhǎng)階段，產(chǎn)業(yè)發(fā)展波動(dòng)較大，經(jīng)歷兩次明顯起落，總體上產(chǎn)業(yè)增長(zhǎng)比較緩慢.這主要因?yàn)楦母镩_(kāi)放初期，山東省經(jīng)濟(jì)總量規(guī)模小，資本、技術(shù)等要素短缺，

產(chǎn)業(yè)結(jié)構(gòu)調(diào)整緩慢，改革開(kāi)放各方面的政策尚不協(xié)調(diào)，市場(chǎng)機(jī)制與政府調(diào)控交織影響.1992～2008年為高速增長(zhǎng)階段，1992年鄧小平同志發(fā)表了南巡講話(huà).與此同時(shí)，山東省市場(chǎng)經(jīng)濟(jì)體系逐漸完善，以及各方面改革政策的施行促進(jìn)了產(chǎn)業(yè)發(fā)展高速增長(zhǎng).1992～1994年產(chǎn)業(yè)發(fā)展年均增長(zhǎng)率高達(dá)4.28%，1994～1996年有所下降.1996～2008年開(kāi)始回升，受1998年亞洲“金融風(fēng)暴”影響，回升較慢，年均增長(zhǎng)率為2.40%，這說(shuō)明，隨著山東省產(chǎn)業(yè)發(fā)展市場(chǎng)化程度的加深，產(chǎn)業(yè)發(fā)展的國(guó)際化水平提升，由此帶來(lái)的產(chǎn)業(yè)發(fā)展對(duì)外依賴(lài)性明顯增強(qiáng).2008～2015年為持續(xù)增長(zhǎng)階段，產(chǎn)業(yè)發(fā)展年均增長(zhǎng)率為1.97%，2008年全球金融危機(jī)后，山東省產(chǎn)業(yè)發(fā)展進(jìn)入新的轉(zhuǎn)型期，產(chǎn)業(yè)增長(zhǎng)速度較前一階段明顯放緩，但依然維持較高的增長(zhǎng)速度.根據(jù)回歸模型，至2020年山東省產(chǎn)業(yè)綜合發(fā)展水平將繼續(xù)提高且速度將有所加快.說(shuō)明山東省產(chǎn)業(yè)發(fā)展基本擺脫了2008年金融危機(jī)的影響，產(chǎn)業(yè)發(fā)展趨于強(qiáng)勁.

山東省產(chǎn)業(yè)綜合發(fā)展水平不斷提升的同時(shí)，產(chǎn)業(yè)發(fā)展系統(tǒng)結(jié)構(gòu)也不斷演化.根據(jù)表2及公式（2），結(jié)合統(tǒng)計(jì)數(shù)據(jù)，計(jì)算出1978～2015年產(chǎn)業(yè)發(fā)展子系統(tǒng)綜合指數(shù)，并據(jù)此繪制1985，1995，2005及2015年產(chǎn)業(yè)發(fā)展系統(tǒng)結(jié)構(gòu)演化雷達(dá)圖（圖2）.

1985年，子系統(tǒng)差異較大，除產(chǎn)業(yè)結(jié)構(gòu)、產(chǎn)業(yè)效益及產(chǎn)業(yè)和諧程度外，其余子系統(tǒng)均處于起步階段.改革開(kāi)放初期，以聯(lián)產(chǎn)承包責(zé)任制為核心的農(nóng)村經(jīng)濟(jì)體制改革率先實(shí)施，促進(jìn)了農(nóng)村經(jīng)濟(jì)繁榮發(fā)展和農(nóng)民收入收入水平提升，產(chǎn)業(yè)和諧程度綜合指數(shù)超過(guò)0.6，達(dá)到研究期內(nèi)最高水平.

1995年，產(chǎn)業(yè)發(fā)展子系統(tǒng)差異減小，產(chǎn)業(yè)外向度、產(chǎn)業(yè)綠色化水平及產(chǎn)業(yè)結(jié)構(gòu)提升較快，產(chǎn)業(yè)效益及和諧程度呈現(xiàn)衰減態(tài)勢(shì).深化企業(yè)改革，進(jìn)一步擴(kuò)大對(duì)外開(kāi)放，積極推動(dòng)外向型經(jīng)濟(jì)等政策的施行，促進(jìn)了山東省產(chǎn)業(yè)結(jié)構(gòu)升級(jí)和外向型發(fā)展，大量新技術(shù)的應(yīng)用提升了產(chǎn)業(yè)發(fā)展的綠色化水平.

2005年，產(chǎn)業(yè)發(fā)展子系統(tǒng)差異進(jìn)一步減小，產(chǎn)業(yè)結(jié)構(gòu)和產(chǎn)業(yè)科技創(chuàng)新能力訊速提高，產(chǎn)業(yè)外向度則呈現(xiàn)衰減態(tài)勢(shì).這個(gè)時(shí)期社會(huì)主義市場(chǎng)經(jīng)濟(jì)體制初步建立，山東省以公有制為主體多種所有制經(jīng)濟(jì)共同發(fā)展的產(chǎn)業(yè)發(fā)展格局基本形成，一大批優(yōu)勢(shì)產(chǎn)業(yè)在市場(chǎng)競(jìng)爭(zhēng)中發(fā)展壯大，原有的基礎(chǔ)設(shè)施“瓶頸”基本消除，“科教興魯”戰(zhàn)略實(shí)施效果明顯.

2015年，產(chǎn)業(yè)外向度進(jìn)一步衰減，綜合指數(shù)低于0.2，其余子系統(tǒng)綜合指數(shù)不斷提高.2008年“金融危機(jī)”帶來(lái)的外向型產(chǎn)業(yè)發(fā)展受阻，產(chǎn)業(yè)外向度降低.產(chǎn)業(yè)長(zhǎng)期發(fā)展積累下的矛盾和問(wèn)題，促使產(chǎn)業(yè)在科技創(chuàng)新、綠色發(fā)展、和諧發(fā)展、轉(zhuǎn)型升級(jí)等方面不斷提高.

研究期內(nèi)，產(chǎn)業(yè)發(fā)展系統(tǒng)結(jié)構(gòu)逐漸優(yōu)化，除產(chǎn)業(yè)效益及和諧程度整體下降，產(chǎn)業(yè)發(fā)展其余子系統(tǒng)均得到不同程度的提升，但提升幅度存在較大差異.產(chǎn)業(yè)規(guī)模、結(jié)構(gòu)、創(chuàng)新能力提升幅度較高，產(chǎn)業(yè)綠色化水平次之，產(chǎn)業(yè)外向度先增加后降低.可見(jiàn)山東省產(chǎn)業(yè)發(fā)展正逐漸轉(zhuǎn)變?yōu)榧a(chǎn)業(yè)創(chuàng)新發(fā)展、協(xié)調(diào)發(fā)展、綠色發(fā)展為一體的產(chǎn)業(yè)綜合發(fā)展.目前主要面臨質(zhì)量提升和外向型經(jīng)濟(jì)發(fā)展的雙重壓力.

4產(chǎn)業(yè)發(fā)展機(jī)制分析

為探討各子系統(tǒng)對(duì)產(chǎn)業(yè)綜合發(fā)展的影響，運(yùn)用eviews6.0構(gòu)建產(chǎn)業(yè)綜合發(fā)展及其子系統(tǒng)VAR模型，運(yùn)用脈沖響應(yīng)及方差分解分析產(chǎn)業(yè)綜合發(fā)展指數(shù)對(duì)各子系正向沖擊產(chǎn)生的響應(yīng)及子系統(tǒng)的平均貢獻(xiàn)率.

4.1模型與方法

VAR模型是把系統(tǒng)中每個(gè)內(nèi)生變量作為系統(tǒng)所有內(nèi)生變量滯后值的函數(shù)來(lái)構(gòu)造模型，從而將單變量自回歸模型推廣到由多元時(shí)間序列變量組成的“向量”自回歸模型[23-24].

其數(shù)學(xué)表達(dá)式為：

yt=A1yt-1+…+Apyt-p+Bxt+εt，t=1，2，…，T，（3）

式中，yt是h維的內(nèi)生變量向量，xt是k維的外生變量向量，p是滯后階數(shù)，T是樣本個(gè)數(shù)，h×h維矩陣A1，…，AP和h×h維矩陣B是被估計(jì)的系數(shù)矩陣，εt是k維的擾動(dòng)向量，其可以同期相關(guān)但不與自己的滯后期及等式右邊的變量相關(guān)[25-26].

假設(shè)模型外生變量為常數(shù)項(xiàng)，內(nèi)生變量為產(chǎn)業(yè)發(fā)展及其子系統(tǒng)綜合指數(shù).為消除可能存在的異方差獲得平穩(wěn)性的時(shí)間序列，對(duì)所有數(shù)據(jù)進(jìn)行對(duì)數(shù)化處理（表3）.

4.2模型驗(yàn)證

VAR模型的檢驗(yàn)重點(diǎn)是模型整體的穩(wěn)定性水平，只有在模型穩(wěn)定的基礎(chǔ)上，才能利用脈沖響應(yīng)和方差分解來(lái)研究隨機(jī)擾動(dòng)對(duì)變量系統(tǒng)的動(dòng)態(tài)沖擊[27].

以時(shí)間序列數(shù)據(jù)建立VAR模型，首先要保證變量是平穩(wěn)的，否則將產(chǎn)生偽回歸現(xiàn)象[28].運(yùn)用ADF檢驗(yàn)方法對(duì)LnY、LnX1、LnX2、LnX3、LnX4、LnX5、LnX6以及LnX7進(jìn)行單位根檢驗(yàn)（表4）.

可以看出，時(shí)間序列LnY、LnX1、LnX3、LnX4及LnX6為非平穩(wěn)序列，對(duì)一階差分后的序列進(jìn)行ADF檢驗(yàn)，均顯示為平穩(wěn)序列，說(shuō)明時(shí)間序列一階單整.因此，建立DLnY、DLnX1、DLnX2、DLnX3、DLnX4、DLnX5、DLnX6以及DLnX7的VAR模型來(lái)分析產(chǎn)業(yè)綜合發(fā)展與其子系統(tǒng)之間的動(dòng)態(tài)關(guān)系.

對(duì)時(shí)間序列數(shù)據(jù)差分后，通過(guò)協(xié)整檢驗(yàn)確定各變量之間是否存在長(zhǎng)期穩(wěn)定的均衡關(guān)系.首先確定協(xié)整檢驗(yàn)的滯后階數(shù)，表5為VAR模型直至最大滯后階數(shù)的各種信息標(biāo)準(zhǔn)，可以看出，VAR模型的最優(yōu)滯后階數(shù)為2，則協(xié)整檢驗(yàn)的最優(yōu)滯后階數(shù)為1.

在5%的顯著水平下，DLnY、DLnX1、DLnX2、DLnX3、DLnX4、DLnX5、DLnX6以及DLnX7之間存在協(xié)整關(guān)系，且至少有7個(gè)協(xié)整方程（表6）.

隨后可建立VAR模型，為達(dá)到分析各變量之間關(guān)系的目的，還須對(duì)整個(gè)模型進(jìn)行平穩(wěn)性檢驗(yàn).

VAR模型的全部特征根全部落在單位圓內(nèi)（圖3），說(shuō)明該VAR模型系統(tǒng)穩(wěn)定且后續(xù)分析結(jié)果合理，可以在此基礎(chǔ)上，進(jìn)行產(chǎn)業(yè)發(fā)展脈沖響應(yīng)及方差分解分析.

4.3產(chǎn)業(yè)發(fā)展影響機(jī)制分析

脈沖響應(yīng)能夠刻畫(huà)出某個(gè)變量的變動(dòng)或沖擊對(duì)其它變量產(chǎn)生影響的軌跡，并沖過(guò)脈沖響應(yīng)圖來(lái)展現(xiàn)每個(gè)影響因素的影響過(guò)程[26，27].圖4為產(chǎn)業(yè)發(fā)展對(duì)其子系統(tǒng)的脈沖響應(yīng)路徑圖.其中，實(shí)線(xiàn)表示脈沖響應(yīng)函數(shù)曲線(xiàn)，虛線(xiàn)表示正負(fù)兩倍標(biāo)準(zhǔn)差偏離帶曲線(xiàn).

圖4（a）為產(chǎn)業(yè)綜合發(fā)展對(duì)其自身的脈沖響應(yīng)路徑.給產(chǎn)業(yè)綜合發(fā)展一個(gè)正單位的沖擊，其自身初期基本保持正效應(yīng)但逐步衰減，第5.5期之后保持穩(wěn)定的負(fù)效應(yīng).這表明，產(chǎn)業(yè)綜合發(fā)展水平的提高對(duì)其自身初期表現(xiàn)為帶動(dòng)作用，但作用強(qiáng)度逐漸降低，之后表現(xiàn)為抑制作用.改革開(kāi)放初期，產(chǎn)業(yè)不斷發(fā)展集聚了大量人力、物力、財(cái)力，促使產(chǎn)業(yè)進(jìn)一步向高水平躍遷，發(fā)展到一定規(guī)模后，市場(chǎng)需求趨于飽和，勞動(dòng)力、自然資源、技術(shù)等要素的供給約束性增強(qiáng)，產(chǎn)業(yè)發(fā)展趨于緩慢，對(duì)自身發(fā)展表現(xiàn)為抑制作用.

圖4（b）為產(chǎn)業(yè)綜合發(fā)展對(duì)產(chǎn)業(yè)規(guī)模的脈沖響應(yīng)路徑.給產(chǎn)業(yè)規(guī)模一個(gè)正單位的沖擊，產(chǎn)業(yè)綜合發(fā)展一直保持正效應(yīng)，初期響應(yīng)較強(qiáng)且波動(dòng)性較大，第8期之后變化平穩(wěn).這表明，產(chǎn)業(yè)規(guī)模的擴(kuò)大對(duì)山東省產(chǎn)業(yè)綜合發(fā)展有明顯的帶動(dòng)作用，且早期影響顯著，之后影響程度明顯降低.因此，在產(chǎn)業(yè)發(fā)展已經(jīng)達(dá)到較高水平的情況下，不能夠再單一依靠數(shù)量擴(kuò)張促進(jìn)產(chǎn)業(yè)發(fā)展.

圖4（c）為產(chǎn)業(yè)綜合發(fā)展對(duì)產(chǎn)業(yè)結(jié)構(gòu)的脈沖響應(yīng)路徑.給產(chǎn)業(yè)結(jié)構(gòu)一個(gè)正單位的沖擊，產(chǎn)業(yè)綜合發(fā)展的正效應(yīng)存在一定滯后期，前6.5期響應(yīng)不穩(wěn)定且初期表現(xiàn)為負(fù)效應(yīng)，之后一直保持穩(wěn)定的正效應(yīng).這說(shuō)明，山東省應(yīng)將產(chǎn)業(yè)結(jié)構(gòu)轉(zhuǎn)型升級(jí)作為區(qū)域經(jīng)濟(jì)發(fā)展的長(zhǎng)期戰(zhàn)略，不能急于求成.

圖4（d）為產(chǎn)業(yè)綜合發(fā)展對(duì)產(chǎn)業(yè)效益的脈沖響應(yīng)路徑.給產(chǎn)業(yè)效益一個(gè)正單位的沖擊，產(chǎn)業(yè)發(fā)展在前2期保持逐漸增強(qiáng)的正效應(yīng)，之后不斷衰減，第4～5期保持負(fù)效應(yīng)，之后一直保持微弱的正效應(yīng).這表明，產(chǎn)業(yè)效益對(duì)產(chǎn)業(yè)綜合發(fā)展的帶動(dòng)難以長(zhǎng)期維持.

圖4（e）為產(chǎn)業(yè)綜合發(fā)展對(duì)產(chǎn)業(yè)發(fā)展和諧程度的脈沖響應(yīng)路徑.給產(chǎn)業(yè)發(fā)展協(xié)調(diào)程度一個(gè)正單位的沖擊，產(chǎn)業(yè)綜合發(fā)展一直保持負(fù)效應(yīng)，且初期強(qiáng)度較高.這表明，農(nóng)村經(jīng)濟(jì)體制改革，盡管增加了農(nóng)民收入，縮小了城鄉(xiāng)收入差距，但并沒(méi)有從根本上成為二三產(chǎn)業(yè)發(fā)展的主要驅(qū)動(dòng)力，且目前山東省城鄉(xiāng)“二元化”及其帶來(lái)的二三產(chǎn)業(yè)產(chǎn)值結(jié)構(gòu)與就業(yè)結(jié)構(gòu)偏離依然制約著產(chǎn)業(yè)的進(jìn)一步發(fā)展.

圖4（f）為產(chǎn)業(yè)綜合發(fā)展對(duì)產(chǎn)業(yè)科技創(chuàng)新能力的脈沖響應(yīng)路徑.給產(chǎn)業(yè)科技創(chuàng)新能力一個(gè)正單位的沖擊，產(chǎn)業(yè)綜合發(fā)展的正效應(yīng)存在一定滯后期，前2.5期為負(fù)效應(yīng)，之后一直保持穩(wěn)定的正效應(yīng).這表明，產(chǎn)業(yè)科技創(chuàng)新能力的提高對(duì)產(chǎn)業(yè)綜合發(fā)展的帶動(dòng)作用存在一定滯后期.山東省應(yīng)全面實(shí)施創(chuàng)新驅(qū)動(dòng)戰(zhàn)略，通過(guò)創(chuàng)新驅(qū)動(dòng)推進(jìn)產(chǎn)業(yè)轉(zhuǎn)型升級(jí)，提升產(chǎn)業(yè)效益，促進(jìn)產(chǎn)業(yè)綜合發(fā)展.

圖（g）為產(chǎn)業(yè)綜合發(fā)展對(duì)產(chǎn)業(yè)綠色發(fā)展水平的脈沖響應(yīng)路徑.給產(chǎn)業(yè)綠色發(fā)展一個(gè)正單位的沖擊，產(chǎn)業(yè)綜合發(fā)展在前5期保持正效應(yīng)，第5～8期為負(fù)效應(yīng)，之后一直保持正效應(yīng).這表明，產(chǎn)業(yè)綠色發(fā)展水平的提高對(duì)產(chǎn)業(yè)綜合發(fā)展作用不穩(wěn)定，初期表現(xiàn)為帶動(dòng)作用，之后表現(xiàn)為短暫的抑制作用，之后又表現(xiàn)為帶動(dòng)作用.說(shuō)明山東省產(chǎn)業(yè)綠色化發(fā)展尚未對(duì)產(chǎn)業(yè)整體發(fā)展帶來(lái)顯著且穩(wěn)定的影響.但從長(zhǎng)期來(lái)看，“節(jié)能減排”等綠色發(fā)展對(duì)產(chǎn)業(yè)綜合發(fā)展會(huì)產(chǎn)生積極作用.

圖（h）為產(chǎn)業(yè)綜合發(fā)展對(duì)產(chǎn)業(yè)發(fā)展外向度的脈沖響應(yīng)路徑.給產(chǎn)業(yè)發(fā)展外向度一個(gè)正單位的沖擊，產(chǎn)業(yè)綜合發(fā)展一直保持穩(wěn)定的正效應(yīng).這表明，外向型產(chǎn)業(yè)的發(fā)展對(duì)產(chǎn)業(yè)綜合發(fā)展有明顯的帶動(dòng)作用且作用穩(wěn)定.山東省應(yīng)積極融入“一帶一路”戰(zhàn)略，積極參與國(guó)際競(jìng)爭(zhēng)與合作，促進(jìn)外向型經(jīng)濟(jì)發(fā)展.

綜上，產(chǎn)業(yè)發(fā)展各子系統(tǒng)對(duì)產(chǎn)業(yè)綜合發(fā)展的影響機(jī)制不同，從沖擊響應(yīng)結(jié)果來(lái)看，產(chǎn)業(yè)規(guī)模及外向型產(chǎn)業(yè)對(duì)產(chǎn)業(yè)綜合發(fā)展表現(xiàn)為正向效應(yīng)，產(chǎn)業(yè)和諧程度對(duì)產(chǎn)業(yè)綜合發(fā)展表現(xiàn)為負(fù)向效應(yīng)，產(chǎn)業(yè)效益及產(chǎn)業(yè)綠色化水平對(duì)產(chǎn)業(yè)綜合發(fā)展起初表現(xiàn)為正向效應(yīng)，之后作用不穩(wěn)定，產(chǎn)業(yè)結(jié)構(gòu)和產(chǎn)業(yè)科技創(chuàng)新能力對(duì)產(chǎn)業(yè)綜合發(fā)展的正向效應(yīng)存在一定滯后期.

為進(jìn)一步分析長(zhǎng)期產(chǎn)業(yè)發(fā)展過(guò)程中各子系統(tǒng)貢獻(xiàn)程度大小，對(duì)產(chǎn)業(yè)綜合發(fā)展進(jìn)行方差分解（圖5）.方差分解是通過(guò)分析每個(gè)結(jié)構(gòu)沖擊對(duì)內(nèi)生變量變化的貢獻(xiàn)程度，進(jìn)而評(píng)價(jià)不同結(jié)構(gòu)沖擊的重要性[29].

產(chǎn)業(yè)綜合發(fā)展受自身沖擊的影響呈逐步遞減的趨勢(shì)，從初始的100%遞減到14.50%.各子系統(tǒng)對(duì)產(chǎn)業(yè)綜合發(fā)展的影響整體上都呈增長(zhǎng)趨勢(shì).產(chǎn)業(yè)協(xié)調(diào)程度的貢獻(xiàn)率最高，且影響越來(lái)越顯著.產(chǎn)業(yè)規(guī)模初期影響較強(qiáng)，中后期逐漸減弱.產(chǎn)業(yè)結(jié)構(gòu)和產(chǎn)業(yè)效益的影響先上升后下降，除初期個(gè)別時(shí)段較強(qiáng)外，整體上穩(wěn)定在5.5%左右.產(chǎn)業(yè)科技創(chuàng)新影響逐漸增強(qiáng)，后期貢獻(xiàn)率約為初期的7倍.此外，產(chǎn)業(yè)綠色化及產(chǎn)業(yè)外向度的影響也逐漸增強(qiáng).

從目前來(lái)看，山東省產(chǎn)業(yè)綜合發(fā)展除14.50%受本身決定之外，受子系統(tǒng)影響程度從高到低分別為：產(chǎn)業(yè)發(fā)展協(xié)調(diào)程度，產(chǎn)業(yè)發(fā)展外向度，產(chǎn)業(yè)規(guī)模，產(chǎn)業(yè)科技創(chuàng)新能力，產(chǎn)業(yè)結(jié)構(gòu)，產(chǎn)業(yè)綠色發(fā)展，產(chǎn)業(yè)效益.

5結(jié)論

改革開(kāi)放以來(lái)，山東省產(chǎn)業(yè)發(fā)展水平不斷提高且呈現(xiàn)出典型的階段性特點(diǎn).改革開(kāi)放初期為震蕩緩慢增長(zhǎng)階段，1992～2008年為高速增長(zhǎng)階段，2008～2015年為持續(xù)增長(zhǎng)階段，未來(lái)五年山東省產(chǎn)業(yè)綜合發(fā)展水平將繼續(xù)提高且速度將有所加快.

從產(chǎn)業(yè)綜合發(fā)展系統(tǒng)結(jié)構(gòu)演化來(lái)看，山東省產(chǎn)業(yè)發(fā)展系統(tǒng)結(jié)構(gòu)逐漸優(yōu)化.改革開(kāi)放初期產(chǎn)業(yè)系統(tǒng)結(jié)構(gòu)不平衡，進(jìn)入新世紀(jì)后，產(chǎn)業(yè)發(fā)展逐漸呈現(xiàn)出產(chǎn)業(yè)規(guī)模、結(jié)構(gòu)、效益、創(chuàng)新、協(xié)調(diào)、綠色為一體的產(chǎn)業(yè)綜合發(fā)展態(tài)勢(shì)，目前主要面臨質(zhì)量提升和外向型經(jīng)濟(jì)發(fā)展的雙重壓力.

從沖擊響應(yīng)結(jié)果來(lái)看，產(chǎn)業(yè)規(guī)模及外向型產(chǎn)業(yè)對(duì)產(chǎn)業(yè)綜合發(fā)展表現(xiàn)為正向效應(yīng)，產(chǎn)業(yè)協(xié)調(diào)程度對(duì)產(chǎn)業(yè)綜合發(fā)展表現(xiàn)為負(fù)向效應(yīng)，產(chǎn)業(yè)效益的產(chǎn)業(yè)綠色化對(duì)產(chǎn)業(yè)綜合發(fā)展起初表現(xiàn)為正向效應(yīng)，之后作用不穩(wěn)定，產(chǎn)業(yè)結(jié)構(gòu)和產(chǎn)業(yè)科技創(chuàng)新能力對(duì)產(chǎn)業(yè)綜合發(fā)展的正向效應(yīng)存在一定滯后期.

預(yù)測(cè)方差結(jié)果表明，產(chǎn)業(yè)發(fā)展本身及各子系統(tǒng)對(duì)產(chǎn)業(yè)綜合發(fā)展都有一定的貢獻(xiàn)度，產(chǎn)業(yè)綜合發(fā)展除14.50%受本身決定之外，受子系統(tǒng)影響程度從高到低分別為：產(chǎn)業(yè)發(fā)展和諧程度，產(chǎn)業(yè)發(fā)展外向度，產(chǎn)業(yè)規(guī)模，產(chǎn)業(yè)科技創(chuàng)新能力，產(chǎn)業(yè)結(jié)構(gòu)，產(chǎn)業(yè)綠色發(fā)展，產(chǎn)業(yè)效益.

產(chǎn)業(yè)發(fā)展系統(tǒng)復(fù)雜，包括產(chǎn)業(yè)規(guī)模、結(jié)構(gòu)、效益等多方面，其對(duì)產(chǎn)業(yè)綜合發(fā)展的貢獻(xiàn)率和影響機(jī)制在產(chǎn)業(yè)發(fā)展的不同階段差異明顯.山東省應(yīng)在適度合理控制產(chǎn)業(yè)規(guī)模的基礎(chǔ)上，把對(duì)外開(kāi)放與滿(mǎn)足內(nèi)需有效結(jié)合起來(lái)，積極融入國(guó)家“一帶一路”戰(zhàn)略，推動(dòng)外向型產(chǎn)業(yè)新一輪發(fā)展.并采取有效措施，處理好產(chǎn)業(yè)和諧發(fā)展問(wèn)題，盡量減少其對(duì)產(chǎn)業(yè)綜合發(fā)展的負(fù)向效應(yīng)，將產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化升級(jí)和提高產(chǎn)業(yè)科技創(chuàng)新能力納入長(zhǎng)遠(yuǎn)目標(biāo)，通過(guò)創(chuàng)新驅(qū)動(dòng)推進(jìn)產(chǎn)業(yè)升級(jí).

在指標(biāo)體系的構(gòu)建和數(shù)學(xué)模型的選擇上，存在一定的局限性，僅僅從全省整體視角出發(fā)，探討山東省產(chǎn)業(yè)發(fā)展演變及其影響機(jī)制，基于地市或縣域空間尺度的地域分異視角研究有待進(jìn)一步深入.

參考文獻(xiàn)

[1]李佳校， 武新乾. “十三五”時(shí)期洛陽(yáng)第一產(chǎn)業(yè)產(chǎn)值灰色系統(tǒng)預(yù)測(cè)[J]. 農(nóng)村經(jīng)濟(jì)與科技， 2016， 27（1）：119-120.

[2]章奇. 第二產(chǎn)業(yè)結(jié)構(gòu)調(diào)整與科技金融聯(lián)動(dòng)發(fā)展的關(guān)系研究[J]. 科學(xué)管理研究， 2016（3）：109-112.

[3]謝敏， 趙紅巖， 朱娜娜，等. 浙江省第三產(chǎn)業(yè)空間集聚特征與成因[J]. 經(jīng)濟(jì)地理， 2015， 35（9）：96-102.

[4]龍冬平， 李同昇， 苗園園，等. 中國(guó)農(nóng)業(yè)現(xiàn)代化發(fā)展水平空間分異及類(lèi)型[J]. 地理學(xué)報(bào)， 2014， 69（2）：213-226.

[5]黃群慧.“新常態(tài)”、工業(yè)化后期與工業(yè)增長(zhǎng)新動(dòng)力[J]. 中國(guó)工業(yè)經(jīng)濟(jì)， 2014（10）：5-19.

[6]吳艷， 張訓(xùn)， 翟歡歡. 戰(zhàn)略性新興產(chǎn)業(yè)發(fā)展的金融支持研究[J]. 經(jīng)濟(jì)數(shù)學(xué)， 2013， 30（1）：105-110.

[7]曹虹劍， 李康. 中國(guó)高技術(shù)產(chǎn)業(yè)自主創(chuàng)新能力的影響因素[J]. 經(jīng)濟(jì)數(shù)學(xué)， 2016， 33（3）：77-82.

[8]洪涓， 劉甦， 孫黛琳，等. 北京與倫敦文化創(chuàng)意產(chǎn)業(yè)發(fā)展比較研究[J]. 城市問(wèn)題， 2013（6）：38-41.

[9]楊家偉， 喬家君. 河南省產(chǎn)業(yè)結(jié)構(gòu)演進(jìn)與機(jī)理探究[J]. 經(jīng)濟(jì)地理， 2013， 33（9）：93-100.

[10]宋周鶯， 劉衛(wèi)東. 西部地區(qū)產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化路徑分析[J]. 中國(guó)人口·資源與環(huán)境， 2013， 23（10）：31-37.

[11]楊智峰， 陳霜華， 汪偉. 中國(guó)產(chǎn)業(yè)結(jié)構(gòu)變化的動(dòng)因分析——基于投入產(chǎn)出模型的實(shí)證研究[J]. 財(cái)經(jīng)研究， 2014， 40（9）.

[12]高遠(yuǎn)東， 張衛(wèi)國(guó)， 陽(yáng)琴. 中國(guó)產(chǎn)業(yè)結(jié)構(gòu)高級(jí)化的影響因素研究[J]. 經(jīng)濟(jì)地理， 2015， 35（6）：96-101.

[13]何平， 陳丹丹， 賈喜越. 產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化研究[J]. 統(tǒng)計(jì)研究， 2014， 31（7）：31-37.

[14]劉媛媛， 孫慧. 資源型產(chǎn)業(yè)集群形成機(jī)理分析與實(shí)證[J]. 中國(guó)人口·資源與環(huán)境， 2014， 24（11）：103-111.

[15]萬(wàn)幼清， 張妮， 魯平俊. 產(chǎn)業(yè)集群協(xié)同創(chuàng)新風(fēng)險(xiǎn)及其形成機(jī)理研究[J]. 管理世界， 2015（2）：182-183.

[16]趙春明. “一帶一路”戰(zhàn)略與中國(guó)綠色產(chǎn)業(yè)發(fā)展[J]. 學(xué)海， 2016， 2016（1）：137-142.

[17]李文， 馬潤(rùn)平. 金融支持“一帶一路”綠色產(chǎn)業(yè)發(fā)展效度檢驗(yàn)-基于17個(gè)省份的面板數(shù)據(jù)模型[J]. 甘肅金融， 2016（11）：34-35.

[18]朱江麗，李子聯(lián).長(zhǎng)三角城市群產(chǎn)業(yè)—人口—空間耦合協(xié)調(diào)發(fā)展研究[J].中國(guó)人口·資源與環(huán)境，2015，25（2）： 75-82.

[19]李存貴.中國(guó)省域城鄉(xiāng)產(chǎn)業(yè)協(xié)調(diào)發(fā)展綜合評(píng)價(jià)[J]. 統(tǒng)計(jì)與決策，2016（9）： 61-64.

[20]梁威，劉滿(mǎn)鳳.戰(zhàn)略性新興產(chǎn)業(yè)與區(qū)域經(jīng)濟(jì)耦合協(xié)調(diào)發(fā)展研究——以江西省為例[J].華東經(jīng)濟(jì)管理，2016，30（5）： 14-19.

[21]李道芳，陳琳，趙恒志.產(chǎn)業(yè)發(fā)展現(xiàn)狀評(píng)價(jià)的指標(biāo)體系及實(shí)證研究[J].中國(guó)行政管理，2007（09）：70-72.

[22]喬家君.改進(jìn)的熵值法在河南省可持續(xù)發(fā)展能力評(píng)估中的應(yīng)用[J].資源科學(xué). 2004，26（1）： 113-119.

[23]ALLEN P G， MORZUCH B J. Twenty-five years of progress， problems， and conflicting evidence in econometric forecasting. What about the next 25 years？[J]. International journal of forecasting， 2006， 22（3）：475-492.

[24]高鐵梅.計(jì)量經(jīng)濟(jì)分析方法與建模[M].北京：清華大學(xué)出版社， 2009.

[25]樊重俊.貝葉斯向量自回歸分析方法及其應(yīng)用[J].數(shù)理統(tǒng)計(jì)與管理， 2010， 29（6）：1060-1066.

[26]朱慧明，韓玉啟，鄭進(jìn)城.基于正態(tài)-Gamma共軛先驗(yàn)分布的貝葉斯AR（p）預(yù)測(cè)模型[J].統(tǒng)計(jì)與決策， 2005（2）：8-9.

[27]朱慧明，劉智偉.時(shí)間序列向量自回歸模型的貝葉斯推斷理論[J]. 統(tǒng)計(jì)與決策， 2004（1）：11-12.

[28]王瑞鵬，馮曉華.基于VAR模型的新疆城市化動(dòng)力機(jī)制研究[J]. 企業(yè)經(jīng)濟(jì)， 2013（2）：149-153.

[29]高鐵梅.計(jì)量經(jīng)濟(jì)分析方法與建模：EVIEWS應(yīng)用及實(shí)例[M].北京：清華大學(xué)出版社，2005.