山東省工業(yè)化、城鎮(zhèn)化與農(nóng)業(yè)現(xiàn)代化互動(dòng)關(guān)系研究

吳宗杰，劉廣亮，董會忠

（山東理工大學(xué)商學(xué)院，山東淄博255012）

●華東經(jīng)濟(jì)

山東省工業(yè)化、城鎮(zhèn)化與農(nóng)業(yè)現(xiàn)代化互動(dòng)關(guān)系研究

吳宗杰，劉廣亮，董會忠

（山東理工大學(xué)商學(xué)院，山東淄博255012）

文章選取山東省17地市2002-2014年面板數(shù)據(jù)，采用DEA－Malmquist指數(shù)法對農(nóng)業(yè)全要素生產(chǎn)率進(jìn)行測算并以其作為表示農(nóng)業(yè)現(xiàn)代化的指標(biāo)，運(yùn)用面板向量自回歸模型估計(jì)了工業(yè)化、城鎮(zhèn)化和農(nóng)業(yè)現(xiàn)代化（簡稱“三化”）的動(dòng)態(tài)關(guān)系。研究表明：城鎮(zhèn)化和工業(yè)化之間存在顯著的相互促進(jìn)作用，而農(nóng)業(yè)現(xiàn)代化與兩者（城鎮(zhèn)化和工業(yè)化）之間不存在顯著的關(guān)系；城鎮(zhèn)化與工業(yè)化互為格蘭杰因果關(guān)系，農(nóng)業(yè)現(xiàn)代化與兩者均不存在格蘭杰因果關(guān)系；工業(yè)化和城鎮(zhèn)化兩者與農(nóng)業(yè)現(xiàn)代化互動(dòng)關(guān)系較弱導(dǎo)致農(nóng)業(yè)現(xiàn)代化的發(fā)展滯后。在此基礎(chǔ)上提出了推進(jìn)山東省“三化”同步發(fā)展的政策建議。

工業(yè)化；城鎮(zhèn)化；DEA－Malmquist；面板向量自回歸

一、引言

工業(yè)化、城鎮(zhèn)化與農(nóng)業(yè)現(xiàn)代化發(fā)展水平是國家現(xiàn)代化進(jìn)程重要的考察方面，實(shí)現(xiàn)三者協(xié)調(diào)發(fā)展是促進(jìn)我國現(xiàn)代化進(jìn)程的必要前提。我國“十三五”規(guī)劃指出，農(nóng)業(yè)發(fā)展是我國全面建成小康社會的基礎(chǔ)，穩(wěn)步深化工業(yè)化、城鎮(zhèn)化進(jìn)程的同時(shí)必須加快轉(zhuǎn)變農(nóng)業(yè)現(xiàn)代化發(fā)展。“三化”均衡發(fā)展有利于實(shí)現(xiàn)國民經(jīng)濟(jì)穩(wěn)定健康發(fā)展，促進(jìn)我國社會經(jīng)濟(jì)產(chǎn)業(yè)結(jié)構(gòu)升級轉(zhuǎn)型，如果“三化”發(fā)展失衡必然會引起經(jīng)濟(jì)結(jié)構(gòu)性矛盾，導(dǎo)致經(jīng)濟(jì)發(fā)展停滯現(xiàn)象。當(dāng)前，我國依舊是農(nóng)業(yè)大國，農(nóng)業(yè)人口基數(shù)大，農(nóng)村的大量勞動(dòng)力和資金外流導(dǎo)致農(nóng)業(yè)現(xiàn)代化發(fā)展面臨嚴(yán)峻考驗(yàn)；工業(yè)發(fā)展處于中后期，產(chǎn)業(yè)結(jié)構(gòu)升級、調(diào)整等問題突出；城鎮(zhèn)化水平逐步深化，穩(wěn)步地推動(dòng)城鎮(zhèn)化發(fā)展能夠有效地解決城鄉(xiāng)二元結(jié)構(gòu)問題。因此，實(shí)現(xiàn)“三化”協(xié)調(diào)發(fā)展對于我國經(jīng)濟(jì)社會轉(zhuǎn)型具有重要的時(shí)代意義。

近年來，政府和學(xué)術(shù)界對“三化”協(xié)調(diào)發(fā)展的相關(guān)問題進(jìn)行了大量研究與探索。劉易斯（1954）提出二元經(jīng)濟(jì)結(jié)構(gòu)理論，該理論假定工業(yè)代表著技術(shù)先進(jìn)的現(xiàn)代化部門，農(nóng)業(yè)代表著自給自足的傳統(tǒng)部門，農(nóng)業(yè)生產(chǎn)要素逐步向城市轉(zhuǎn)移，農(nóng)業(yè)部門長期為工業(yè)部門提供剩余勞動(dòng)力和農(nóng)業(yè)剩余價(jià)值［1］。城鎮(zhèn)化在農(nóng)業(yè)和工業(yè)互動(dòng)發(fā)展過程中的重要性逐步體現(xiàn)。錢納里（1989）的研究表明，隨著工業(yè)化進(jìn)程的推進(jìn)，人口向城市集聚促進(jìn)了城鎮(zhèn)化，城鎮(zhèn)化加快了工業(yè)和農(nóng)業(yè)互動(dòng)發(fā)展［2］。黃祖輝等（2013）基于我國“三化”發(fā)展現(xiàn)狀與特征，從經(jīng)濟(jì)學(xué)視角分析“三化”互動(dòng)關(guān)系，指出三者協(xié)調(diào)發(fā)展是區(qū)域經(jīng)濟(jì)協(xié)調(diào)發(fā)展的內(nèi)在要求［3］；曾福生等（2013）通過對城市與農(nóng)村區(qū)域工業(yè)水平對應(yīng)的勞動(dòng)工資、技術(shù)水平進(jìn)行對比研究，發(fā)現(xiàn)勞動(dòng)技能多元化、工資水平適中的區(qū)域“三化”發(fā)展的協(xié)調(diào)性較高［4-5］；王貝（2011）通過構(gòu)建VAR模型研究“三化”發(fā)展之間的動(dòng)態(tài)關(guān)聯(lián)關(guān)系，指出我國農(nóng)業(yè)現(xiàn)代化的發(fā)展與工業(yè)化和城鎮(zhèn)化的發(fā)展不同步，農(nóng)業(yè)現(xiàn)代化發(fā)展滯后問題會隨著工業(yè)化和城鎮(zhèn)化的深化而凸顯出來（姜會民等，2012；Weerdt J等，2013）［6-8］。發(fā)展滯后的農(nóng)業(yè)化對于工業(yè)化和城鎮(zhèn)化的進(jìn)程起到重要的影響作用［9］。因此，郭家虎等（2007）認(rèn)為我國需要進(jìn)行農(nóng)村要素市場化改革，實(shí)現(xiàn)城鄉(xiāng)之間要素的自由流動(dòng)，倡導(dǎo)工業(yè)反哺農(nóng)業(yè)，減少城鄉(xiāng)差異，實(shí)現(xiàn)城鄉(xiāng)協(xié)調(diào)發(fā)展［10］；可以說，在工業(yè)化進(jìn)程中，農(nóng)業(yè)的發(fā)展需要以先進(jìn)技術(shù)為依靠，以先進(jìn)生產(chǎn)力為依托，以政府扶持為保障（Hofmann A等，2013）［11］。

已有的文獻(xiàn)對“三化”研究在經(jīng)濟(jì)學(xué)解釋、理論分析和實(shí)證方面取得了一定成果。但是也存在一些不足：首先，農(nóng)業(yè)現(xiàn)代化發(fā)展水平僅選取農(nóng)村人均機(jī)械總動(dòng)力、土地利用率和農(nóng)業(yè)勞動(dòng)生產(chǎn)率這些指標(biāo)中某一指標(biāo)表示，不能全面地反映農(nóng)業(yè)現(xiàn)代化水平，使得整個(gè)文章的分析基礎(chǔ)缺乏科學(xué)性；其次，研究三者之間動(dòng)態(tài)關(guān)系的文獻(xiàn)多是基于VAR模型，VAR模型要求數(shù)據(jù)時(shí)序較長且研究結(jié)果會受到時(shí)間序列波動(dòng)大和自相關(guān)等問題干擾而造成偏差，而PVAR模型能夠運(yùn)用時(shí)序較短的面板數(shù)據(jù)并且可以有效地避免變量內(nèi)生性引起的偏差。因此本文與之前研究有所不同，通過DEA－Malmquist方法對農(nóng)業(yè)全要素生產(chǎn)率進(jìn)行測算，并以其表示農(nóng)業(yè)現(xiàn)代化，并構(gòu)建PVAR模型探究“三化”之間的互動(dòng)影響關(guān)系，為山東省“三化”協(xié)調(diào)發(fā)展提出政策建議。

二、“三化”互動(dòng)發(fā)展的作用機(jī)理

工業(yè)化、城鎮(zhèn)化和農(nóng)業(yè)現(xiàn)代化是一個(gè)多重關(guān)聯(lián)的系統(tǒng)，三者之間存在著促進(jìn)、配合、協(xié)調(diào)甚至抑制的多重關(guān)系，“三化”協(xié)調(diào)發(fā)展的實(shí)現(xiàn)需要三個(gè)子系統(tǒng)共同作用、相互配合（見圖1）。

農(nóng)業(yè)現(xiàn)代化是“三化”系統(tǒng)發(fā)展的基礎(chǔ)：農(nóng)業(yè)現(xiàn)代化的發(fā)展能夠提高勞動(dòng)生產(chǎn)率，單位勞動(dòng)力產(chǎn)出增加，促進(jìn)農(nóng)村剩余勞動(dòng)力向城市轉(zhuǎn)移。一方面農(nóng)副產(chǎn)品是工業(yè)化和城鎮(zhèn)化發(fā)展所必需的物質(zhì)資源；另一方面，農(nóng)村勞動(dòng)力向城市轉(zhuǎn)移能夠加快城鎮(zhèn)化進(jìn)程（蘇發(fā)金，2012）［12］。如果推進(jìn)工業(yè)化和城鎮(zhèn)化進(jìn)程是以犧牲農(nóng)村經(jīng)濟(jì)發(fā)展為代價(jià)，則難以實(shí)現(xiàn)三者可持續(xù)發(fā)展（姚旭兵，2016）［13］。工業(yè)化是“三化”系統(tǒng)發(fā)展的主導(dǎo)，是整個(gè)系統(tǒng)發(fā)展的驅(qū)動(dòng)力。工業(yè)化進(jìn)程穩(wěn)步推進(jìn)，能夠?yàn)檗r(nóng)村剩余勞動(dòng)力向城市轉(zhuǎn)移提供就業(yè)機(jī)會，同時(shí)為城鎮(zhèn)化和農(nóng)業(yè)現(xiàn)代化發(fā)展提供工業(yè)產(chǎn)品的供給（董梅生，2014）［14］；工業(yè)化進(jìn)程擠占了部分農(nóng)村資源，如果不能協(xié)調(diào)兩者發(fā)展關(guān)系則容易擴(kuò)大城鄉(xiāng)發(fā)展差距（夏顯力，2013）［15］。城鎮(zhèn)化是“三化”系統(tǒng)發(fā)展的載體與媒介，城鎮(zhèn)化發(fā)展加快了要素之間的流動(dòng)，勞動(dòng)人口轉(zhuǎn)移是工業(yè)發(fā)展勞動(dòng)要素的保障；城鎮(zhèn)化發(fā)展需要完善基礎(chǔ)設(shè)施，擴(kuò)大工業(yè)品消費(fèi)需求；城鎮(zhèn)化的發(fā)展有利于擴(kuò)大農(nóng)產(chǎn)品的消費(fèi)和提高農(nóng)業(yè)人均耕地面積（夏春萍，2012）［16］。

當(dāng)工業(yè)化和城鎮(zhèn)化穩(wěn)步發(fā)展時(shí)，可以有效地促進(jìn)農(nóng)村冗余勞動(dòng)力的轉(zhuǎn)移，擴(kuò)大農(nóng)副產(chǎn)品的需求，能夠?qū)崿F(xiàn)農(nóng)副產(chǎn)品和農(nóng)村冗余勞動(dòng)資源等要素優(yōu)化配置，對農(nóng)業(yè)現(xiàn)代化發(fā)展有正效應(yīng)。如果工業(yè)化和城鎮(zhèn)化的進(jìn)程過快，城市規(guī)模過分?jǐn)U張，加劇農(nóng)業(yè)部門的土地、資本等生產(chǎn)要素向非農(nóng)產(chǎn)業(yè)集聚，造成農(nóng)業(yè)現(xiàn)代化發(fā)展資源匱乏的困境，則對農(nóng)業(yè)現(xiàn)代化發(fā)展有負(fù)效應(yīng)。因此有必要研究三者的互動(dòng)關(guān)系，為山東省“三化”同步發(fā)展提出政策建議。

圖1 “三化”互動(dòng)發(fā)展作用機(jī)理

三、研究方法

（一）農(nóng)業(yè)現(xiàn)代化發(fā)展績效評價(jià)模型

包絡(luò)分析法（DEA）是一種非參數(shù)的效率評價(jià)方法，通過運(yùn)用包絡(luò)線來替代生產(chǎn)函數(shù)對全要素生產(chǎn)率進(jìn)行測算。Caves（1982）等人將這一理論與Malmquist生產(chǎn)率指數(shù)相結(jié)合，使其被廣泛應(yīng)用于生產(chǎn)率的測算［17］。DEA-Malmquist在測算全要素生產(chǎn)率的同時(shí)客觀地描述不同時(shí)期的全要素生產(chǎn)率變動(dòng)情況。因此，本文基于距離函數(shù)構(gòu)造Malmquist指數(shù)如下：

（二）“三化”互動(dòng)分析模型

面板數(shù)據(jù)向量自回歸模型（Panel-Vector-atuo-Regression）是Holtz等學(xué)者在VAR的基礎(chǔ)上提出的。PVAR主要用于預(yù)測和解釋相互關(guān)聯(lián)的序列系統(tǒng)中各個(gè)內(nèi)生變量的變動(dòng)造成其他內(nèi)生變量的動(dòng)態(tài)影響，模型估計(jì)運(yùn)用動(dòng)態(tài)面板GMM（Generalized method of moments）做參數(shù)估計(jì)。Holtz-Eakin認(rèn)為通常較長時(shí)間的數(shù)據(jù)是VAR模型取得較好效果的基礎(chǔ)，由于我國統(tǒng)計(jì)制度的原因造成許多區(qū)域數(shù)據(jù)時(shí)序較短不符合要求，甚至很多地市的宏觀經(jīng)濟(jì)數(shù)據(jù)公布的時(shí)間不足10年，PVAR模型能有效解決數(shù)據(jù)時(shí)序較短的問題［18-19］。相對于VAR模型，PVAR模型能夠充分利用面板數(shù)據(jù)信息量大、自由度高、共線性低的優(yōu)點(diǎn)，控制時(shí)間效應(yīng)和地區(qū)效應(yīng)而避免有偏估計(jì)，從而有效地避免變量內(nèi)生性引起的偏差（連玉君等，2007）［20］。其計(jì)算公式為：

其中，yit為內(nèi)生變量向量；yit-n為yit的滯后項(xiàng)；yit中的i表示市區(qū)；t表示時(shí)間；A1,…,An表示待估系數(shù)矩陣；μit表示擾動(dòng)項(xiàng)。

四、實(shí)證研究

（一）研究區(qū)域概況

山東省經(jīng)濟(jì)基礎(chǔ)雄厚，重工業(yè)發(fā)展迅速，輕工業(yè)特別是紡織和食品工業(yè)相當(dāng)發(fā)達(dá)，山東省擁有規(guī)模以上工業(yè)企業(yè)3.8萬家，工業(yè)化水平高于全國平均水平。2014年山東省常住人口9 789.43萬人，城鎮(zhèn)人口5 385.17萬人，城鎮(zhèn)化率高于全國平均水平。雖然山東省是農(nóng)業(yè)大省，但是2014年第一產(chǎn)業(yè)就業(yè)比為30.6%，第一產(chǎn)業(yè)生產(chǎn)總值占總生產(chǎn)值的8%，說明近1/3的勞動(dòng)力只創(chuàng)造了8%的生產(chǎn)總值，即農(nóng)業(yè)現(xiàn)代化發(fā)展相對滯后。因此，測度山東省農(nóng)業(yè)現(xiàn)代化水平并研究山東省工業(yè)化、城鎮(zhèn)化、農(nóng)業(yè)現(xiàn)代化互動(dòng)發(fā)展關(guān)系，對于山東省統(tǒng)籌城鄉(xiāng)發(fā)展，促進(jìn)工業(yè)化、城鎮(zhèn)化、農(nóng)業(yè)現(xiàn)代化良性互動(dòng)具有一定的現(xiàn)實(shí)意義。

（二）變量選取與數(shù)據(jù)來源

許多學(xué)者對工業(yè)化、城鎮(zhèn)化與農(nóng)業(yè)現(xiàn)代化三者的評價(jià)指標(biāo)及體系進(jìn)行了研究，既有學(xué)者構(gòu)建綜合性指標(biāo)體系對其測算，也有學(xué)者以單個(gè)的代表性的指標(biāo)對其衡量。由于本文的研究圍繞山東省17地市展開，考慮到數(shù)據(jù)的可得性，本文只對農(nóng)業(yè)現(xiàn)代化進(jìn)行測算，工業(yè)化和城鎮(zhèn)化用單一指標(biāo)來表示。

農(nóng)業(yè)現(xiàn)代化測算變量。農(nóng)業(yè)產(chǎn)出變量：第一產(chǎn)業(yè)是指以利用自然力為主，生產(chǎn)不必經(jīng)過深度加工就可消費(fèi)的產(chǎn)品的部門。根據(jù)我國國家統(tǒng)計(jì)局對第一產(chǎn)業(yè)的劃分規(guī)定，第一產(chǎn)業(yè)是指農(nóng)業(yè)主要包括林業(yè)、牧業(yè)、漁業(yè)等，本文參照周端明（2009）的做法選取第一產(chǎn)業(yè)產(chǎn)值表示農(nóng)業(yè)產(chǎn)出（以山東2000年不變價(jià)格計(jì)算）［21］。參照多數(shù)學(xué)者的研究，本文選取勞動(dòng)投入、播種面積、機(jī)械化、化肥以及灌溉投入作為農(nóng)業(yè)投入變量。勞動(dòng)投入選取年底鄉(xiāng)村就業(yè)人員總數(shù)表示；播種面積以各市糧食作物合計(jì)播種面積表示；機(jī)械化選取農(nóng)村人均機(jī)械總動(dòng)力表示；化肥以實(shí)際農(nóng)用化肥施用量（折純計(jì)算）表示，包括氮、磷、鉀和復(fù)合肥等；灌溉投入以每年實(shí)際有效灌溉面積表示。

工業(yè)化通常被狹義地定義為第二產(chǎn)業(yè)產(chǎn)值（或工業(yè)就業(yè)人數(shù)）占國民生產(chǎn)總值（或就業(yè)總?cè)藬?shù)）的比重不斷上升的過程，但工業(yè)化進(jìn)程不是孤立的，是與農(nóng)業(yè)現(xiàn)代化和服務(wù)業(yè)發(fā)展相輔相成的。因此參照董梅生等（2014）的研究，本文選取第二、三產(chǎn)業(yè)增加值占GDP的比重表示工業(yè)化［14］。城鎮(zhèn)化是指隨著一個(gè)國家或區(qū)域工業(yè)化進(jìn)程中，科學(xué)技術(shù)的進(jìn)步、產(chǎn)業(yè)結(jié)構(gòu)調(diào)整和社會生產(chǎn)力發(fā)展促進(jìn)人口從農(nóng)業(yè)部門向非農(nóng)業(yè)部門轉(zhuǎn)移的歷史過程，本文借鑒王貝（2011）做法選取城鎮(zhèn)人口占總?cè)丝诘谋戎乇硎境擎?zhèn)化［6］。

原數(shù)據(jù)主要來源于《山東省統(tǒng)計(jì)年鑒》（2003-2015），部分?jǐn)?shù)據(jù)來源于山東省十七地市統(tǒng)計(jì)年鑒（2003-2015年）。

（三）農(nóng)業(yè)全要素生產(chǎn)率的測算結(jié)果

將農(nóng)業(yè)產(chǎn)出指標(biāo)和農(nóng)業(yè)投入指標(biāo)代入到DEAP2.1軟件中進(jìn)行求解，得到山東省17地市歷年的農(nóng)業(yè)全要素生產(chǎn)率（以其表示農(nóng)業(yè)現(xiàn)代化），見表1所列。

表1 山東省17地市2003-2014年的Malmquist指數(shù)

（四）“三化”的互動(dòng)關(guān)系研究

1.面板數(shù)據(jù)平穩(wěn)性檢驗(yàn)

考慮到取對數(shù)可以消除數(shù)據(jù)的異方差而不改變數(shù)據(jù)的穩(wěn)健性，分別對工業(yè)化、城鎮(zhèn)化和農(nóng)業(yè)現(xiàn)代化取對數(shù)表示為lnind、lnur和lnam。變量的平穩(wěn)性是進(jìn)行PVAR模型估計(jì)和脈沖響應(yīng)分析的基礎(chǔ)，分別采用ADF和LLC兩種檢驗(yàn)方法對面板數(shù)據(jù)進(jìn)行穩(wěn)健性檢驗(yàn)，見表2所列。由表2可以看出，在5%的顯著水平下，統(tǒng)計(jì)值均拒絕原假設(shè)，說明lnind、lnur和lnalp三個(gè)變量是平穩(wěn)的，符合PVAR模型分析要求。

表2 山東省“三化”指標(biāo)穩(wěn)健性檢驗(yàn)結(jié)果

2.PVAR滯后階數(shù)選擇

滯后階數(shù)選擇狀況見表3所列。

表3 滯后階數(shù)選擇

參照AIC、BIC和HQIC準(zhǔn)則和時(shí)序長度，確定PVAR模型滯后階數(shù)為1。

3.PVAR模型估計(jì)

估計(jì)PVAR模型不同于時(shí)序VAR模型，需要消除樣本的固定效應(yīng)，本文運(yùn)用向前均值差分法避免一般均值差分方法可能引起的偏誤。根據(jù)選取的滯后階數(shù)，基于GMM估計(jì)得到了工業(yè)化、城鎮(zhèn)化和農(nóng)業(yè)化的PVAR模型估計(jì)結(jié)果，見表4所列。

表4 山東省“三化”面板VAR（1）模型的參數(shù)估計(jì)結(jié)果

從第二列的回歸系數(shù)可以看出：滯后1階的城鎮(zhèn)化和工業(yè)化回歸系數(shù)顯著為正（在5%的水平上），滯后1階農(nóng)業(yè)現(xiàn)代化的回歸系數(shù)不顯著，表明滯后1階的城鎮(zhèn)化和工業(yè)化對工業(yè)化有顯著的正向作用，滯后1階農(nóng)業(yè)現(xiàn)代化對工業(yè)化的影響不顯著。從第三列的回歸系數(shù)可以看出：滯后1階的城鎮(zhèn)化和工業(yè)化回歸系數(shù)顯著為正（在5%的水平上），滯后1階農(nóng)業(yè)現(xiàn)代化的回歸系數(shù)不顯著，表明滯后1階的城鎮(zhèn)化和工業(yè)化對城鎮(zhèn)化有顯著的正向作用，滯后1階農(nóng)業(yè)現(xiàn)代化對城鎮(zhèn)化的影響不顯著。從第四列的回歸系數(shù)可以看出：滯后1階的城鎮(zhèn)化、工業(yè)化和農(nóng)業(yè)現(xiàn)代化的回歸系數(shù)均不顯著，表明滯后1階工業(yè)化、城鎮(zhèn)化和農(nóng)業(yè)現(xiàn)代化對農(nóng)業(yè)現(xiàn)代化的影響均不顯著。

從三者的回歸系數(shù)可以看出：城鎮(zhèn)化和工業(yè)化之間存在顯著的相互促進(jìn)的作用，而農(nóng)業(yè)現(xiàn)代化與兩者（城鎮(zhèn)化和工業(yè)化）之間不存在顯著的關(guān)系，表明城鎮(zhèn)化和工業(yè)化之間存在良好的互動(dòng)關(guān)系，而農(nóng)業(yè)現(xiàn)代化與兩者的互動(dòng)關(guān)系較弱。工業(yè)的快速發(fā)展擠占了部分農(nóng)業(yè)資源，工業(yè)發(fā)展能夠加快城鎮(zhèn)化進(jìn)程；城鎮(zhèn)化水平的提高為工業(yè)化提供消費(fèi)需求和生產(chǎn)要素的保障，兩者共同發(fā)展。工業(yè)發(fā)展對農(nóng)業(yè)發(fā)展具有兩種效應(yīng)：促進(jìn)效應(yīng)，工業(yè)發(fā)展為農(nóng)業(yè)發(fā)展提供化肥、農(nóng)藥等方面的支持，促進(jìn)農(nóng)業(yè)現(xiàn)代化實(shí)現(xiàn)機(jī)械化和高效化；擠占效應(yīng)，工業(yè)化發(fā)展引起農(nóng)村資源向城市轉(zhuǎn)移，工業(yè)發(fā)展?fàn)帄Z了部分農(nóng)業(yè)資源，體現(xiàn)在土地要素、勞動(dòng)要素、資本要素等方面。當(dāng)擠占效應(yīng)大于促進(jìn)效應(yīng)或者兩者中和時(shí)，城鎮(zhèn)化和工業(yè)化與農(nóng)業(yè)現(xiàn)代化之間的互動(dòng)較弱，引起農(nóng)業(yè)現(xiàn)代化發(fā)展滯后。

4.因果關(guān)系檢驗(yàn)

上述的PVAR（1）模型的參數(shù)估計(jì)結(jié)果可以確定山東省“三化”之間的相關(guān)性，利用Granger因果檢驗(yàn)進(jìn)一步明確這三個(gè)變量間的因果關(guān)系，見表5所列。

表5 山東省“三化”指標(biāo)因果關(guān)系檢驗(yàn)結(jié)果

表5的結(jié)果顯示：城鎮(zhèn)化是工業(yè)化的格蘭杰原因。一方面，城鎮(zhèn)化提高需要完善城市基礎(chǔ)設(shè)施且日常生活需要消費(fèi)大量的工業(yè)產(chǎn)品，擴(kuò)大了工業(yè)需求；另一方面，城鎮(zhèn)化提高為工業(yè)生產(chǎn)提供了大量的資源、勞動(dòng)要素。農(nóng)業(yè)現(xiàn)代化不是工業(yè)化的格蘭杰原因，是由于山東工業(yè)化早期發(fā)展是建立在政府主導(dǎo)機(jī)制下引起的。工業(yè)化是城鎮(zhèn)化的格蘭杰原因，工業(yè)化發(fā)展為勞動(dòng)力從農(nóng)村區(qū)域向城市區(qū)域轉(zhuǎn)移提供了就業(yè)機(jī)會，促進(jìn)了城鎮(zhèn)化水平的提高。農(nóng)業(yè)現(xiàn)代化不是城鎮(zhèn)化的格蘭杰原因，說明山東省城鎮(zhèn)化與農(nóng)業(yè)現(xiàn)代化的協(xié)調(diào)發(fā)展需要加強(qiáng)。城鎮(zhèn)化和工業(yè)化均不是農(nóng)業(yè)現(xiàn)代化的格蘭杰原因，說明山東省工業(yè)化和城鎮(zhèn)化發(fā)展擠占了農(nóng)業(yè)資源，山東省人多地少，農(nóng)業(yè)科技人才相對其他行業(yè)十分匱乏，導(dǎo)致城鄉(xiāng)發(fā)展差距擴(kuò)大，造成山東省部分區(qū)域農(nóng)業(yè)現(xiàn)代化長期落后于城鎮(zhèn)化、工業(yè)化。格蘭杰因果檢驗(yàn)進(jìn)一步驗(yàn)證了PVAR（1）模型的參數(shù)估計(jì)結(jié)果。

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

脈沖響應(yīng)函數(shù)表示在PVAR模型中，其他變量保持不變的情況下，某一變量t期的擾動(dòng)項(xiàng)發(fā)生變動(dòng)，會引起模型中其他變量之間發(fā)生一系列交互作用。脈沖響應(yīng)函數(shù)用于描述整個(gè)系統(tǒng)對于沖擊擾動(dòng)在不同滯后期的反映關(guān)系，即用于預(yù)測一個(gè)變量在不同滯后期對另一個(gè)變量的作用時(shí)滯。本文在PVAR（1）模型的基礎(chǔ)上，使用Monte Carlo方法模擬200次，得到每個(gè)變量的沖擊對各個(gè)變量滯后0-6期的作用（見圖2）。橫軸表示沖擊反映期數(shù)，最大值為6期，縱軸表示沖擊反映程度，中間的線表示各期反映程度，上下線為95%的置信區(qū)間。

圖2 lnalp、lnind和lnur的脈沖響應(yīng)曲線

由圖2第1列城鎮(zhèn)化對其他變量沖擊反映曲線可以看出：城鎮(zhèn)化對自身的新息沖擊反映呈現(xiàn)降低趨勢，從第一期0.1降到第6期的0.001；城鎮(zhèn)化對農(nóng)業(yè)現(xiàn)代化的新沖擊的反映置信區(qū)間均包括零，表明在95%的置信水平下不能證明農(nóng)業(yè)現(xiàn)代化對城鎮(zhèn)化信息沖擊能夠引起城鎮(zhèn)化反映；城鎮(zhèn)化對工業(yè)化的信息沖擊反映較弱。由圖2第2列農(nóng)業(yè)現(xiàn)代化對其他變量沖擊反映曲線可以看出：農(nóng)業(yè)現(xiàn)代化對于城鎮(zhèn)化和工業(yè)化的信息沖擊均為負(fù)數(shù)，表明需要加強(qiáng)城鎮(zhèn)化、工業(yè)化與農(nóng)業(yè)現(xiàn)代化的聯(lián)系；農(nóng)業(yè)現(xiàn)代化對自身的沖擊反映在第一期存在顯著正向作用。由圖2第3列工業(yè)化對其他變量沖擊反映曲線可以看出：城鎮(zhèn)化和農(nóng)業(yè)現(xiàn)代化對工業(yè)化沖擊反映曲線的置信區(qū)間均包括零，說明在95%的置信區(qū)間內(nèi)不能證明農(nóng)業(yè)現(xiàn)代化和城鎮(zhèn)化對工業(yè)化的信息沖擊會引起工業(yè)化反映；工業(yè)化對自身沖擊反映在1-6期均大于0，說明工業(yè)化對自身發(fā)展有顯著促進(jìn)作用。

五、結(jié)論與建議

本文基于山東省17地市2002-2014年面板數(shù)據(jù)，采用DEA－Malmquist指數(shù)法對農(nóng)業(yè)全要素生產(chǎn)率進(jìn)行測算并以其作為農(nóng)業(yè)現(xiàn)代化的指標(biāo)，運(yùn)用面板向量自回歸模型估計(jì)了工業(yè)化、城鎮(zhèn)化和農(nóng)業(yè)現(xiàn)代化的動(dòng)態(tài)關(guān)系。研究表明工業(yè)化與城鎮(zhèn)化之間存在顯著的互動(dòng)關(guān)系，農(nóng)業(yè)現(xiàn)代化與兩者的關(guān)系不顯著；工業(yè)化與城鎮(zhèn)化存在雙向格蘭杰原因，農(nóng)業(yè)現(xiàn)代化與兩者均不存在格蘭杰原因。各變量的脈沖響應(yīng)分析表明，對未來的預(yù)測各變量對自身沖擊反映占了很大部分。城鎮(zhèn)化和工業(yè)化與農(nóng)業(yè)現(xiàn)代化互動(dòng)較少是造成農(nóng)業(yè)現(xiàn)代化發(fā)展滯后的主要原因。

根據(jù)以上研究提出以下建議：

充分發(fā)揮山東科技優(yōu)勢，引進(jìn)國內(nèi)外先進(jìn)的農(nóng)業(yè)科技成果，依托現(xiàn)代科技促進(jìn)傳統(tǒng)農(nóng)業(yè)向現(xiàn)代化農(nóng)業(yè)轉(zhuǎn)變，實(shí)現(xiàn)生物技術(shù)和體制機(jī)制方面的創(chuàng)新，提高農(nóng)業(yè)全要素生產(chǎn)率。充分利用山東工業(yè)化和城鎮(zhèn)化發(fā)展積累的現(xiàn)實(shí)條件，用現(xiàn)代物質(zhì)裝備農(nóng)業(yè)，提高農(nóng)業(yè)技術(shù)裝備現(xiàn)代化水平，進(jìn)一步提高農(nóng)業(yè)現(xiàn)代化發(fā)展水平。優(yōu)化山東省產(chǎn)業(yè)結(jié)構(gòu)，農(nóng)業(yè)現(xiàn)代化發(fā)展滯后在一定程度上限制了山東省經(jīng)濟(jì)綜合發(fā)展水平，因此需要加大農(nóng)業(yè)政策扶持力度，完善農(nóng)村的衛(wèi)生、醫(yī)療、教育等社會保障體系。加強(qiáng)農(nóng)業(yè)現(xiàn)代化與工業(yè)化、城鎮(zhèn)化的互動(dòng)聯(lián)系，實(shí)施工業(yè)反哺農(nóng)業(yè)戰(zhàn)略部署，建立扶持農(nóng)業(yè)的長效機(jī)制，以期實(shí)現(xiàn)三化協(xié)調(diào)發(fā)展。

［1］Lewis W A.Economic Development with Unlimited Supply of Labor［J］.The Manchester School of Economic and So?cial Studies，1954，22（2）:131-191.

［2］H·錢納里，M·塞爾昆.發(fā)展的格局［M］.北京：中國財(cái)政經(jīng)濟(jì)出版社，1989.

［3］黃祖輝，邵峰，朋文歡.推進(jìn)工業(yè)化、城鎮(zhèn)化和農(nóng)業(yè)現(xiàn)代化協(xié)調(diào)發(fā)展［J］.中國農(nóng)村經(jīng)濟(jì)，2013（1）：8-14.

［4］曾福生，高鳴.中國農(nóng)業(yè)現(xiàn)代化、工業(yè)化和城鎮(zhèn)化協(xié)調(diào)發(fā)展及其影響因素分析——基于現(xiàn)代農(nóng)業(yè)視角［J］.中國農(nóng)村經(jīng)濟(jì)，2013（1）：24-39.

［5］Young A.Inequality，the Urban-Rural Gap，and Migration［J］.The Quarterly Journal of Economics，2013，128（4）：1727-1785.

［6］王貝.中國工業(yè)化、城鎮(zhèn)化和農(nóng)業(yè)現(xiàn)代化關(guān)系實(shí)證研究［J］.城市問題，2011，194（9）：21-25.

［7］姜會民，王振華.吉林省工業(yè)化、城鎮(zhèn)化和農(nóng)業(yè)現(xiàn)代化關(guān)系實(shí)證分析［J］.地理科學(xué)，2012，32（5）：591-595.

［8］Christiaensen L，Weerdt J，Todo Y.Urbanization and Pover?ty Reduction：The Role of Rural Diversification and Second?ary Towns［J］.Agricultural Economics，2013，44（4/5）：435-447.

［9］蔡昉.城鄉(xiāng)收入差距與制度變革的臨界點(diǎn)［J］.中國社會科學(xué)，2003（5）：16-26.

［10］郭家虎.工業(yè)反哺農(nóng)業(yè)長效機(jī)制：收入反哺向要素反哺的過渡［J］.改革，2007，116（12）：13-17.

［11］Hofmann A，Wan G.Determinants of Urbanization［R］. ADB Economics Working Paper Series，2013.

［12］蘇發(fā)金.工業(yè)化、城鎮(zhèn)化與農(nóng)業(yè)現(xiàn)代化：基于VAR模型［J］.統(tǒng)計(jì)與決策，2012（11）：147-150.

［13］姚旭兵，羅光強(qiáng)，寧瑞芳.城鎮(zhèn)化與農(nóng)業(yè)經(jīng)濟(jì)增長的區(qū)域效應(yīng)研究［J］.西南大學(xué)學(xué)報(bào)，2016，42（3）：60-68.

［14］董梅生，楊德才.工業(yè)化、信息化、城鎮(zhèn)化和農(nóng)業(yè)現(xiàn)代化互動(dòng)關(guān)系研究——基于VAR模型［J］.農(nóng)業(yè)技術(shù)經(jīng)濟(jì)2014（4）：14-24.

［15］夏顯力，郝晶輝.陜西省工業(yè)化、城鎮(zhèn)化與農(nóng)業(yè)現(xiàn)代化互動(dòng)關(guān)系的實(shí)證研究［J］.華中農(nóng)業(yè)大學(xué)學(xué)報(bào)，2013（1）：19-24.

［16］夏春萍，劉文清.農(nóng)業(yè)現(xiàn)代化與城鎮(zhèn)化、工業(yè)化協(xié)調(diào)發(fā)展關(guān)系的實(shí)證研究基于VAR模型的計(jì)量分析［J］.農(nóng)業(yè)技術(shù)經(jīng)濟(jì)，2012（5）：79-85.

［17］Caves D，Christensen L，Diewert W E.The Economic Theo?ry of Index Numbers and the Measurement of Input，Out?put and Productivity［J］.Econometrica，1982，50（6）：1393-1414.

［18］Holtz-Eakin D，Newey W，Rosen H.EstimatingVector Au?toregressions with Panel Data［J］.Econometrica，1998，56（6）：1371-1395.

［19］李小勝，張煥明.中國經(jīng)濟(jì)增長、污染排放與能源消費(fèi)間動(dòng)態(tài)關(guān)系研究［J］.山東財(cái)經(jīng)大學(xué)學(xué)報(bào)，2013，35（13）：25-30.

［20］連玉君，程建.投資-現(xiàn)金流敏感性：融資約束還是代理成本？［J］.財(cái)經(jīng)研究，2007，33（2）：37-46.

［21］周端明.技術(shù)進(jìn)步、技術(shù)效率與中國農(nóng)業(yè)生產(chǎn)率增長［J］.農(nóng)業(yè)技術(shù)經(jīng)濟(jì)，2009（12）：70-82.

A Study of the Dynamic Relationship between Industrialization, Urbanization and Agricultural Modernization in Shandong Province

WU Zong-jie,LIU Guang-liang，DONG Hui-zhong
（School of Business,Shandong University of Technology,Zibo 255012,China)

Based on the panel data of 17 prefecture-level cities in Shandong province from 2002 to 2014,this paper estimates agricultural total factor productivity by applying the DEA-Malmquist index method and takes it as the indicator of agricultural modernization.The pa?per also analyzes the dynamic relationship of industrialization,urbanization and agricultural modernization by employing the panel vector auto regression model.The results show that:There is a significant mutual promoting effect between urbanization and industrialization,but there is no significant mutual relationship between agricultural modernization and the two(urbanization and industrialization);Industrial?ization and urbanization present mutual Granger causality,whereas there is no Granger relationship between agricultural modernization and the two.The weaker interactive relationship between industrialization,urbanization and agricultural modernization results in the devel?opment of agricultural modernization lagging behind.Accordingly,this paper puts forward policy suggestions on how to push forward the coordinated development of industrialization,urbanization and agricultural in Shandong province.

industrialization;urbanization;DEA-Malmquist;panel vector auto regression

F127

A

1007-5097（2017）08-0012-06

［責(zé)任編輯：余志虎］

10.3969/j.issn.1007-5097.2017.08.002

2016-10-25

國家軟科學(xué)研究計(jì)劃項(xiàng)目（2012GXS4D090）；國家自然科學(xué)基金項(xiàng)目（71371112）；山東省自然科學(xué)基金項(xiàng)目（ZR2012GM020）

吳宗杰（1962-），男，山東淄博人，教授，經(jīng)濟(jì)學(xué)博士，研究方向：低碳經(jīng)濟(jì)；劉廣亮（1989-），男，山東臨沂人，碩士研究生，研究方向：區(qū)域經(jīng)濟(jì)分析與規(guī)劃；董會忠（1968-），男，山東煙臺人，副教授，博士，研究方向：能源系統(tǒng)分析與工業(yè)生態(tài)化。