孫學(xué)濤

(1.山東社會科學(xué)院 農(nóng)村發(fā)展研究所，山東 濟(jì)南 250002；2.中國海洋大學(xué) 經(jīng)濟(jì)學(xué)院，山東 青島 266100)

一、問題提出

城鄉(xiāng)收入差距是全面建成小康社會的難點。改革開放以來農(nóng)民收入顯著提升，2019年中國農(nóng)村居民人均可支配收入為16 021元，是1978年的119倍。但農(nóng)民收入水平快速增長的同時城鄉(xiāng)收入差距也不斷擴(kuò)大，2019年中國城鄉(xiāng)居民收入比達(dá)到2.64。與此同時，黨和政府也高度重視城鄉(xiāng)收入差距問題。2017年中國提出了到2050年實現(xiàn)鄉(xiāng)村全面振興和農(nóng)民富的目標(biāo)，2020年中央一號文件提出了促進(jìn)農(nóng)民持續(xù)增收的具體措施。促進(jìn)農(nóng)民增收、縮小城鄉(xiāng)收入差距成為實施鄉(xiāng)村振興戰(zhàn)略的關(guān)鍵所在。因此，在鄉(xiāng)村振興的背景下，研究農(nóng)民收入差距的來源以及促進(jìn)農(nóng)民增收的路徑變得尤為迫切。

在脫貧攻堅大背景下，針對城鄉(xiāng)居民收入差距逐漸擴(kuò)大的現(xiàn)象，學(xué)者從經(jīng)濟(jì)發(fā)展的工業(yè)化偏好理論[1-2]和城市偏向理論[3-4]等角度解釋了城鄉(xiāng)收入差距產(chǎn)生的原因。同時，國內(nèi)學(xué)者也從戶籍制度[5]和政治制度[6]等角度解釋了城鄉(xiāng)收入差距產(chǎn)生的原因?，F(xiàn)有文獻(xiàn)主要從制度角度分析城鄉(xiāng)收入差距，但城鄉(xiāng)收入差距問題產(chǎn)生的根源在于城市居民收入相對較高且增速快，而農(nóng)村居民收入水平相對較低且增速慢。因此不僅需要從城鄉(xiāng)收入差距的“城”角度解釋城鄉(xiāng)收入差距，而且還需要從城鄉(xiāng)收入差距的“鄉(xiāng)”這一角度解釋城鄉(xiāng)收入差距產(chǎn)生的原因[7]。

由農(nóng)業(yè)經(jīng)濟(jì)學(xué)理論可知，農(nóng)業(yè)機(jī)械化不僅是提高農(nóng)業(yè)生產(chǎn)率和增加農(nóng)民收入最為重要的動力[8]，也是縮小城鄉(xiāng)收入差距的主要途徑[2]。農(nóng)業(yè)機(jī)械化水平提升可以實現(xiàn)農(nóng)業(yè)機(jī)械對勞動要素的替代[9]，提高勞動生產(chǎn)率[10]，進(jìn)而增加農(nóng)民的農(nóng)業(yè)經(jīng)營性收入。與此同時，農(nóng)業(yè)機(jī)械對勞動要素的替代還為農(nóng)民的非農(nóng)就業(yè)創(chuàng)造了機(jī)會，間接地促進(jìn)了城市化[11]，增加了農(nóng)民的工資性收入，進(jìn)而縮小了城鄉(xiāng)收入差距，最終實現(xiàn)鄉(xiāng)村振興。但這些文獻(xiàn)僅局限于直觀的判斷，缺乏實證分析，特別是農(nóng)業(yè)機(jī)械化對城鄉(xiāng)居民收入差距的具體作用機(jī)理和影響路徑等方面缺乏深入研究。

雖然農(nóng)業(yè)機(jī)械化促進(jìn)農(nóng)民增收的觀點已經(jīng)成為學(xué)界的普遍共識，但在脫貧攻堅背景下農(nóng)業(yè)機(jī)械化對城鄉(xiāng)居民收入差距作用機(jī)理的研究成果還相對缺乏，特別是現(xiàn)有文獻(xiàn)還不能很好地解決農(nóng)業(yè)機(jī)械化與農(nóng)民收入之間的內(nèi)生性問題。本文將利用2002—2017年中國縣域面板數(shù)據(jù)分析農(nóng)村勞動力轉(zhuǎn)移背景下農(nóng)業(yè)機(jī)械化對農(nóng)民收入和城鄉(xiāng)收入差距的影響。

二、文獻(xiàn)回顧

隨著中等收入陷阱理論的發(fā)展，學(xué)者開始逐漸關(guān)注城鄉(xiāng)收入差距的相關(guān)問題，并就形成原因進(jìn)行了多視角研究。戶籍制度[12]、農(nóng)村金融發(fā)展[13]和教育[14]是引起城鄉(xiāng)收入差距的重要原因。除從二元經(jīng)濟(jì)結(jié)構(gòu)視角研究城鄉(xiāng)收入差距的成因外，部分學(xué)者發(fā)現(xiàn)農(nóng)業(yè)自身發(fā)展水平是影響城鄉(xiāng)居民收入差距的主要因素[10]。改革開放以來中國基尼系數(shù)不斷上升，由1978年的0.18上升到2018年的0.47。學(xué)者研究發(fā)現(xiàn)中國最主要的不平等表現(xiàn)在城鄉(xiāng)收入差距[15]，其中縮小城鄉(xiāng)收入差距的重點在于“鄉(xiāng)”[2]，即通過提高農(nóng)民收入的方式縮小城鄉(xiāng)收入差距。農(nóng)業(yè)機(jī)械化不僅能夠提高農(nóng)民收入[10]，而且還是縮小城鄉(xiāng)收入差距的重要方式[2]。因此，農(nóng)業(yè)機(jī)械化或是農(nóng)民增收和緩解城鄉(xiāng)收入差距最具解釋力的變量之一。

理論上，農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響主要表現(xiàn)在兩個方面。一方面，提高農(nóng)業(yè)機(jī)械化水平意味著農(nóng)業(yè)要素投入增加，進(jìn)而提高農(nóng)業(yè)生產(chǎn)率，從而緩解城鄉(xiāng)收入差距；另一方面，雖然農(nóng)業(yè)發(fā)展方式逐步由粗放型增長轉(zhuǎn)向高質(zhì)量型發(fā)展，但農(nóng)業(yè)部門內(nèi)土地非農(nóng)化現(xiàn)象經(jīng)常發(fā)生。因此，通過增加農(nóng)業(yè)部門要素投入的方式促進(jìn)農(nóng)業(yè)高質(zhì)量發(fā)展很難實現(xiàn)。而在農(nóng)業(yè)經(jīng)濟(jì)高質(zhì)量發(fā)展的背景下，通過優(yōu)化農(nóng)業(yè)部門的要素投入結(jié)構(gòu)，促進(jìn)農(nóng)民增收，緩解城鄉(xiāng)收入差距的路徑是切實可行的?；谝陨侠碚摲治?，部分學(xué)者采用家庭數(shù)據(jù)[16]、省級區(qū)域(不包括香港、澳門和臺灣，本文研究同)的樣本數(shù)據(jù)[10]、縣域數(shù)據(jù)[17]等證實了農(nóng)業(yè)機(jī)械化水平提升促進(jìn)農(nóng)民增收，緩解了城鄉(xiāng)收入差距。

但現(xiàn)有文獻(xiàn)存在不足。一是現(xiàn)有文獻(xiàn)忽視了農(nóng)業(yè)機(jī)械化與城鄉(xiāng)收入差距之間的內(nèi)生性。農(nóng)業(yè)部門內(nèi)部通過投入農(nóng)業(yè)機(jī)械促進(jìn)農(nóng)民增收，進(jìn)而縮小了城鄉(xiāng)收入差距；同時，隨著農(nóng)民收入水平的不斷提高，農(nóng)民會繼續(xù)增加對農(nóng)業(yè)機(jī)械要素的投入，進(jìn)而影響農(nóng)業(yè)機(jī)械化水平。因此只探討農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響必然會面臨著內(nèi)生性問題，進(jìn)而可能對實證結(jié)果的有效性產(chǎn)生損失。二是學(xué)者在探討農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響過程中忽視了農(nóng)業(yè)機(jī)械化的空間相關(guān)性，只考慮到本地區(qū)農(nóng)業(yè)機(jī)械化對本地區(qū)城鄉(xiāng)收入差距的影響，忽視了周邊地區(qū)對其城鄉(xiāng)收入差距的影響。

根據(jù)現(xiàn)有文獻(xiàn)存在的不足，本文從以下幾個方面進(jìn)一步拓展：第一，在探討農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距影響的同時還采用縣域所在地的平均坡度作為農(nóng)業(yè)機(jī)械化的工具變量，利用空間雙重差分模型(SDID)進(jìn)一步研究農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響；第二，由于農(nóng)業(yè)機(jī)械化不僅會對農(nóng)民收入和城鄉(xiāng)收入差距產(chǎn)生直接影響，而且還可能存在著農(nóng)村勞動力轉(zhuǎn)移這一中介效應(yīng)，因此在農(nóng)村勞動力流動背景下討論了農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的作用機(jī)制；第三，選取縣域數(shù)據(jù)探討農(nóng)業(yè)機(jī)械化與城鄉(xiāng)收入差距之間的關(guān)系，空間相關(guān)性要求不僅要分析農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響，而且還需要考慮其空間溢出效應(yīng)[18]。

三、計量模型

如果樣本數(shù)據(jù)不僅依賴于自身，而且還依賴于其他地區(qū)的數(shù)據(jù)，即樣本數(shù)據(jù)存在著一定的空間依賴關(guān)系，需要采用空間計量模型研究樣本數(shù)據(jù)之間的關(guān)系。樣本數(shù)據(jù)的空間依賴性在空間計量模型中表現(xiàn)為變量的自回歸項和滯后項。以往的空間計量模型在處理自回歸項和滯后項時，分別運用空間自回歸模型和空間誤差模型處理，但空間自回歸項和空間誤差項可能同時存在[18]?；诖?，本部分嘗試采用帶空間自回歸誤差項的空間自回歸模型(SARAR模型)討論農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響，因為該模型能夠同時處理空間自回歸項和空間誤差項。

SARAR模型主要用于研究相鄰地區(qū)之間的相互關(guān)系，即分析自變量對因變量的影響及其空間溢出效應(yīng)，樣本的空間溢出效應(yīng)體現(xiàn)在空間自回歸項上，空間依賴關(guān)系體現(xiàn)在空間誤差項上。SARAR模型的數(shù)學(xué)表達(dá)式為：

Y=ρWY+Xβ+μ

(1)

μ=λWμ+ε

(2)

其中，ρ表示城鄉(xiāng)收入差距的反應(yīng)系數(shù)，衡量的是相鄰縣域的城鄉(xiāng)收入差距對本縣域城鄉(xiāng)收入差距的影響程度；W表示縣域的空間權(quán)重矩陣；Y表示城鄉(xiāng)收入差距；WY表示城鄉(xiāng)收入差距的空間滯后因變量，衡量的是相鄰縣域的變化對本縣域城鄉(xiāng)收入差距的影響程度；X表示自變量；β為解釋變量的回歸系數(shù)；λ表示空間誤差系數(shù)，衡量的是相鄰縣域的城鄉(xiāng)收入差距誤差對本縣域城鄉(xiāng)收入差距的影響程度；Wμ表示空間滯后誤差項，衡量的是其他相鄰縣域的誤差；ε表示殘差擾動項。

為了分析農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響，式(1)和式(2)中的W為地理距離權(quán)重矩陣，其衡量方式為根據(jù)縣域的經(jīng)緯度并運用歐式距離公式計算縣域i與縣域j之間的距離。

四、數(shù)據(jù)來源

本文的研究對象是縣域。包含全國縣域指標(biāo)的統(tǒng)計年鑒只有《中國區(qū)域經(jīng)濟(jì)統(tǒng)計年鑒》(2003—2015年)和《中國縣(市)社會經(jīng)濟(jì)統(tǒng)計年鑒》(2003—2018年)。同時，城鄉(xiāng)收入差距在一定程度上受到價格的影響，因此，從《中國城市統(tǒng)計年鑒》(2003—2018年)收集各地級市的居民消費價格指數(shù)，對涉及價格的指標(biāo)進(jìn)行平減。由于部分省份的相關(guān)指標(biāo)缺失，本文相關(guān)數(shù)據(jù)指標(biāo)還來源于相關(guān)地市的統(tǒng)計公報(2003—2018年)。以上統(tǒng)計年鑒內(nèi)部分指標(biāo)缺失，本文利用省級統(tǒng)計年鑒對相關(guān)指標(biāo)進(jìn)行補充?？h域所在地的平均坡度數(shù)據(jù)主要由ArcGIS 12.0軟件在地圖上抓取，縣域所在地的經(jīng)緯度數(shù)據(jù)和地理距離權(quán)重矩陣的生成過程中所需要的縣域經(jīng)緯度數(shù)據(jù)均來源于經(jīng)緯度網(wǎng)(1)由于基于不同地圖數(shù)據(jù)的經(jīng)緯度值存在一定的差異，本文所采用的經(jīng)緯度為常用地圖網(wǎng)站經(jīng)緯度的平均值。。

本文的被解釋變量為城鄉(xiāng)收入差距，采用農(nóng)村居民人均純收入占城鎮(zhèn)職工平均工資的比重(%)來衡量。由于本文從城鄉(xiāng)收入差距的“鄉(xiāng)”角度討論城鄉(xiāng)收入差距產(chǎn)生的原因，因此被解釋變量還包含了農(nóng)民收入。農(nóng)民收入采用農(nóng)村居民人均純收入(元)取自然對數(shù)來衡量。

本文的解釋變量為農(nóng)業(yè)機(jī)械化水平。農(nóng)業(yè)機(jī)械化水平采用單位耕地機(jī)械總動力(KW/hm2)來衡量。同時為了解決農(nóng)業(yè)機(jī)械化與城鄉(xiāng)收入差距之間的內(nèi)生性問題，本文還嘗試采用縣域所在地的平均坡度作為農(nóng)業(yè)機(jī)械化的工具變量?？h域所在地的平均坡度數(shù)據(jù)主要由ArcGIS 12.0軟件在地圖上抓取。

控制變量包括縣域互聯(lián)網(wǎng)水平和人力資本等。

變量的統(tǒng)計描述如表1所示。

表1 縣域數(shù)據(jù)統(tǒng)計描述

五、實證分析

(一)空間相關(guān)性檢驗

進(jìn)行空間計量分析之前，首先要分析農(nóng)業(yè)機(jī)械化與城鄉(xiāng)收入差距的空間依賴關(guān)系，如果農(nóng)業(yè)機(jī)械化不僅會對本縣域的城鄉(xiāng)收入差距產(chǎn)生影響，而且還會對相鄰縣域的城鄉(xiāng)收入差距產(chǎn)生影響，才可使用空間計量模型。本文運用莫蘭指數(shù)(Moran’sI)考察農(nóng)業(yè)機(jī)械化的空間相關(guān)關(guān)系?？h域農(nóng)業(yè)機(jī)械化的莫蘭指數(shù)檢驗結(jié)果如表2所示。

表2 莫蘭指數(shù)檢驗結(jié)果

由表2可知，中國縣城農(nóng)業(yè)機(jī)械化的莫蘭指數(shù)值為負(fù)，除2004年和2017年外其余年份均顯著。這表明縣域農(nóng)業(yè)機(jī)械化存在顯著的負(fù)向相關(guān)效應(yīng)，即農(nóng)業(yè)機(jī)械化水平較高的地區(qū)與農(nóng)業(yè)機(jī)械化水平較低的地區(qū)相鄰，說明本縣域農(nóng)業(yè)機(jī)械化水平提升會抑制相鄰縣域農(nóng)業(yè)機(jī)械化水平提升。

(二)空間計量分析

農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響機(jī)制是農(nóng)業(yè)機(jī)械化通過促進(jìn)農(nóng)民增收，進(jìn)而影響城鄉(xiāng)收入差距，因此在討論農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距影響的同時還需要討論農(nóng)業(yè)機(jī)械化對農(nóng)民收入的影響。如果農(nóng)業(yè)機(jī)械化對農(nóng)民收入不存在影響，那么討論農(nóng)業(yè)機(jī)械化與城鄉(xiāng)收入差距就沒有意義。如果農(nóng)業(yè)機(jī)械化對農(nóng)民收入存在影響，那么農(nóng)業(yè)機(jī)械化才有可能影響城鄉(xiāng)收入差距。因此，在分析農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響之前，需要分析農(nóng)業(yè)機(jī)械化對農(nóng)民收入的影響。由于本文采用的數(shù)據(jù)是面板數(shù)據(jù)，實證分析農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響之前需要討論混合效應(yīng)、隨機(jī)效應(yīng)和固定效應(yīng)對于分析農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距影響的適用性。BP檢驗結(jié)果顯著，說明農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響更適合采用隨機(jī)效應(yīng)。豪斯曼(Hausman)檢驗結(jié)果顯著，說明農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響更適合采用固定效應(yīng)。本文采用面板數(shù)據(jù)討論農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響，不僅縣域樣本變動會影響城鄉(xiāng)收入差距，而且時間變動也會對其產(chǎn)生影響。因此，本文嘗試建立的模型為時點個體固定效應(yīng)模型。利用SARAR模型分析農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響，模型的具體估計結(jié)果如表3所示。

表3 農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響

表3給出了農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響，從表3可知。計量模型內(nèi)的自回歸項的系數(shù)顯著為正，說明不僅本地區(qū)的機(jī)械化會影響本地區(qū)的城鄉(xiāng)收入差距，周邊地區(qū)的機(jī)械化也會影響本地區(qū)的城鄉(xiāng)收入差距。由于不同模型之間的結(jié)果基本一致，本部分以縣域樣本為例分析農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響，具體如下：

第一，分析農(nóng)業(yè)機(jī)械化對農(nóng)民收入的影響。由表3可知，農(nóng)業(yè)機(jī)械化對農(nóng)民收入的影響為正，并且通過了顯著性水平檢驗。這說明隨著農(nóng)業(yè)機(jī)械化水平的不斷提高，農(nóng)民收入水平也在不斷提高[10]?？赡艿慕忉屖?，一方面，隨著農(nóng)業(yè)機(jī)械化水平的不斷提高，勞動要素從農(nóng)業(yè)內(nèi)釋放出來，這部分釋放出來的勞動要素流向非農(nóng)產(chǎn)業(yè)，農(nóng)民獲得更多的工資性收入，最終提高了農(nóng)民收入。另一方面，農(nóng)業(yè)機(jī)械化水平提高的過程是資本要素逐漸代替勞動要素的過程。根據(jù)邊際技術(shù)替代率，農(nóng)業(yè)部門內(nèi)農(nóng)業(yè)機(jī)械與農(nóng)業(yè)勞動力之間存在著替代關(guān)系，農(nóng)業(yè)部門增加農(nóng)業(yè)機(jī)械運用，則會減少勞動投入，即會提高勞動產(chǎn)出效率，農(nóng)民進(jìn)而獲得更多的農(nóng)業(yè)經(jīng)營性收入，最終提高農(nóng)民收入。

第二，分析農(nóng)業(yè)機(jī)械化在城鄉(xiāng)收入差距縮小過程中的作用。由表3可知，農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響為正，并且通過了顯著性水平檢驗。這說明隨著農(nóng)業(yè)機(jī)械化水平的不斷提高，會促進(jìn)農(nóng)業(yè)勞動生產(chǎn)率的提高，進(jìn)而縮小收入差距[19]?？赡艿慕忉屖?，一方面，農(nóng)業(yè)機(jī)械化水平提升，意味著農(nóng)業(yè)部門內(nèi)資本要素投入不斷增加，隨著農(nóng)業(yè)部門內(nèi)資本要素和勞動要素投入結(jié)構(gòu)的不斷優(yōu)化，勞動要素生產(chǎn)效率不斷提高，進(jìn)而增加農(nóng)民收入。另一方面，農(nóng)業(yè)機(jī)械補貼等惠農(nóng)政策的實施，意味著非農(nóng)部門向農(nóng)業(yè)部門轉(zhuǎn)移資本要素，農(nóng)業(yè)部門與非農(nóng)部門之間要素投入結(jié)構(gòu)的不斷優(yōu)化也會提高農(nóng)民的收入水平。無論是通過優(yōu)化農(nóng)業(yè)部門內(nèi)部的要素配置結(jié)構(gòu)還是通過優(yōu)化農(nóng)業(yè)部門與非農(nóng)部門之間的要素配置結(jié)構(gòu)，都能促進(jìn)農(nóng)民增收，從而實現(xiàn)縮小收入差距的目標(biāo)。

第三，分析農(nóng)業(yè)機(jī)械化的區(qū)域差異。農(nóng)業(yè)機(jī)械化對東部縣域農(nóng)民收入和城鄉(xiāng)收入差距的影響均為正，并且通過了顯著性水平檢驗；對中部縣域農(nóng)民收入和城鄉(xiāng)收入差距的影響為負(fù)，但不顯著；對西部農(nóng)民收入和城鄉(xiāng)收入差距的影響顯著為負(fù)。實證結(jié)果說明農(nóng)業(yè)機(jī)械化會以提高農(nóng)民收入的方式縮小東部縣域的城鄉(xiāng)收入差距；但農(nóng)業(yè)機(jī)械化不利于西部縣域農(nóng)民收入水平的提升。可能的解釋是，東部地區(qū)平原相對較多，同時地勢起伏也較小，適合農(nóng)業(yè)機(jī)械大規(guī)模作業(yè)，而中、西部縣域的地形地貌使得部分地區(qū)農(nóng)業(yè)機(jī)械無法進(jìn)入。因此，中、西部地區(qū)農(nóng)業(yè)機(jī)械化可能會成為農(nóng)民增收的負(fù)擔(dān)。

第四，分析控制變量對收入的影響。為了準(zhǔn)確地分析農(nóng)業(yè)機(jī)械化對收入的影響，在實證分析過程中加入了控制變量。實證分析過程中信息化水平和醫(yī)療水平等對收入的影響沒有通過顯著水平檢驗；居民消費水平、人力資本水平、金融發(fā)展水平、政府干預(yù)水平和規(guī)模經(jīng)濟(jì)水平等對收入的影響均通過了顯著水平檢驗。限于篇幅，控制變量對收入的影響結(jié)果以及可能的解釋不再一一分析。

六、內(nèi)生性討論

農(nóng)業(yè)機(jī)械化在一定程度上會受到農(nóng)民收入和地形等因素的影響，即農(nóng)業(yè)機(jī)械化不是外生隨機(jī)變量，而是在一定程度上可能是內(nèi)生變量。農(nóng)民收入水平相對較高的地區(qū)農(nóng)業(yè)機(jī)械化可能相對發(fā)達(dá)，而農(nóng)民收入水平相對較低的地區(qū)農(nóng)業(yè)機(jī)械化可能正在發(fā)展中。如果不考慮農(nóng)業(yè)機(jī)械化的內(nèi)生性問題，那么估計出來的農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距影響的系數(shù)只能是農(nóng)業(yè)機(jī)械化與城鄉(xiāng)收入差距之間的相關(guān)關(guān)系，而不是因果關(guān)系。本文嘗試尋找農(nóng)業(yè)機(jī)械化的工具變量以分析其對城鄉(xiāng)收入差距的影響。好的工具變量必須是相關(guān)的和外生的，即工具變量只與農(nóng)業(yè)機(jī)械化相關(guān)，同時工具變量與模型內(nèi)的其他控制變量無關(guān)。基于此，本文借鑒李欣澤等(2017)[20]的研究成果，使用縣域所在地的平均坡度作為農(nóng)業(yè)機(jī)械化的工具變量。因為縣域所在地的平均坡度只通過農(nóng)業(yè)機(jī)械化來影響城鄉(xiāng)收入差距，同時縣域所在地的平均坡度不會對城鄉(xiāng)收入差距產(chǎn)生直接影響。首先，分析縣域所在地的平均坡度的相關(guān)性。農(nóng)業(yè)機(jī)械化與縣域所在地的平均坡度相關(guān)是因為縣域內(nèi)能否大規(guī)模地推行農(nóng)業(yè)機(jī)械會受到地勢和地形的影響。地勢相對平坦的地區(qū)農(nóng)業(yè)機(jī)械大規(guī)模作業(yè)相對容易，而地勢相對復(fù)雜的地區(qū)農(nóng)業(yè)機(jī)械作業(yè)相對困難。其次，分析縣域所在地平均坡度的外生性?？h域所在地的平均坡度屬于自然地理變量。至今還沒有看到相關(guān)文獻(xiàn)證實地區(qū)經(jīng)濟(jì)發(fā)展和農(nóng)民收入水平能夠改變當(dāng)?shù)氐牡貏??；谝陨戏治?，本文運用縣域所在地的平均坡度作為農(nóng)業(yè)機(jī)械化的工具變量來分析農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距影響的結(jié)果，結(jié)果見表4。

表4內(nèi)生性檢驗部分報告了采用縣域所在地的平均坡度作為農(nóng)業(yè)機(jī)械化工具變量的估計結(jié)果。限于篇幅，表4工具變量部分沒有匯報農(nóng)業(yè)機(jī)械化的工具變量對收入差距的影響。由表4的工具變量部分可知，農(nóng)業(yè)機(jī)械化的工具變量對全國縣域和東部縣域農(nóng)民收入的影響為正，并且通過了顯著性水平檢驗；農(nóng)業(yè)機(jī)械化的工具變量對中部縣域農(nóng)民收入的影響為正，但沒有通過顯著性水平檢驗；農(nóng)業(yè)機(jī)械化的工具變量對西部縣域農(nóng)民收入的影響顯著為負(fù)。以上表明縣域坡度可以看成農(nóng)業(yè)機(jī)械化的工具變量。

接下來，分析農(nóng)業(yè)機(jī)械化的工具變量與農(nóng)業(yè)機(jī)械化的差異。將表4工具變量部分與表3農(nóng)民收入部分對比，發(fā)現(xiàn)縣域坡度的系數(shù)大于農(nóng)業(yè)機(jī)械化的系數(shù)。這說明在分析農(nóng)業(yè)機(jī)械化對農(nóng)民收入的影響過程中，可能存在著一定的內(nèi)生性，內(nèi)生性問題會導(dǎo)致實證分析過程中低估了農(nóng)業(yè)機(jī)械化對農(nóng)民收入的影響。

表4 內(nèi)生性檢驗

為了避免樣本的自選擇問題，需要運用農(nóng)業(yè)機(jī)械化水平相對較高的縣域與農(nóng)業(yè)機(jī)械化水平相對較低的縣域進(jìn)行對比，研究農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響。因此本文借鑒范巧(2018)[21]的空間雙重差分法(SDID)將實驗組與控制組相匹配，將控制組作為對照組，研究實驗組與對照組之間的區(qū)別。為了推動農(nóng)業(yè)機(jī)械化，政府于2004年試點實施并于2009年在全國縣域推廣農(nóng)機(jī)購置補貼政策。因此，本部分以農(nóng)機(jī)購置補貼政策作為依據(jù)，基于空間自回歸的雙重差分法研究農(nóng)業(yè)機(jī)械化對農(nóng)民收入影響的模型具體為：

(3)

式(3)內(nèi)的T表示農(nóng)機(jī)購置補貼政策試點實施的年份，2004年之前T=0，2004年以后T=1。D表示縣域是否實施農(nóng)機(jī)購置補貼政策，實施農(nóng)機(jī)購置補貼D=1，否則D=0。其余變量設(shè)置與式(1)和式(2)的變量設(shè)置一致。基于SDID分析的農(nóng)業(yè)機(jī)械化對農(nóng)民收入影響的估計結(jié)果如表4所示。

表4報告了空間雙重差分法估計的農(nóng)機(jī)購置補貼政策對農(nóng)民收入影響的估計結(jié)果，限于表格，本部分沒有匯報農(nóng)機(jī)購置補貼政策對居民收入差距影響的估計結(jié)果。從表4中SDID部分的估計結(jié)果可以看出，農(nóng)機(jī)購置補貼政策對全國縣域和東部縣域農(nóng)民收入的影響顯著為正；對中部縣域農(nóng)民收入的影響為正，但不顯著；對西部縣域農(nóng)民收入的影響顯著為負(fù)。這說明農(nóng)機(jī)購置補貼政策能夠影響農(nóng)民收入水平，且與前文結(jié)果一致。

最后，分析基于不同方法估計的農(nóng)業(yè)機(jī)械化對農(nóng)民收入影響系數(shù)的差異。本部分嘗試討論基于SDID模型的估計系數(shù)與基于SARAR模型的估計系數(shù)的差異。將表4的SDID模型部分與表3的SARAR模型對比發(fā)現(xiàn)，無論是全國縣域還是東部縣域、中部縣域和西部縣域，SARAR模型同樣會低估農(nóng)業(yè)機(jī)械化對農(nóng)民收入的影響?？赡艿慕忉屖荢ARAR模型在分析農(nóng)業(yè)機(jī)械化對農(nóng)民收入的影響過程中將所有的縣域樣本看作一個整體，沒有區(qū)分實施農(nóng)機(jī)購置補貼的縣域與未實施農(nóng)機(jī)購置補貼的縣域。而SDID模型在研究農(nóng)業(yè)機(jī)械化對農(nóng)民收入影響的過程中依據(jù)農(nóng)機(jī)購置補貼政策將實施農(nóng)機(jī)購置補貼政策的縣域與未實施農(nóng)機(jī)購置補貼政策的縣域進(jìn)行匹配，并作為樣本的反事實參照組，去除隨機(jī)因素對農(nóng)民收入的影響。因此基于SARAR模型估計的系數(shù)在一定程度上可能會低估農(nóng)業(yè)機(jī)械化對農(nóng)民收入的影響。

七、作用機(jī)制分析

根據(jù)文獻(xiàn)綜述和前文的分析，本部分將從農(nóng)村勞動力轉(zhuǎn)移視角來研究農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的傳導(dǎo)機(jī)制。在作用機(jī)制分析過程中，本部分借鑒陳詩一和陳登科(2018)[22]的研究成果，將農(nóng)業(yè)機(jī)械化和農(nóng)村勞動力轉(zhuǎn)移滯后一期，分析農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響?？紤]到數(shù)據(jù)的一致性和準(zhǔn)確性，并借鑒伍山林(2016)[23]衡量農(nóng)村勞動力流動的方法，本部分農(nóng)村勞動力流動的衡量方式為縣域從業(yè)人員數(shù)減去第二、三產(chǎn)業(yè)從業(yè)人員數(shù)進(jìn)而得到第一產(chǎn)業(yè)從業(yè)人員數(shù)，并從鄉(xiāng)村從業(yè)人員數(shù)中扣除第一產(chǎn)業(yè)從業(yè)人員數(shù)，進(jìn)而得到農(nóng)村勞動力流動數(shù)量(2)按照該方法計算的2 851個縣級行政單元的農(nóng)村勞動力流動數(shù)量與2009年以后國家統(tǒng)計局公布的《中國農(nóng)民工監(jiān)測調(diào)查報告》中的農(nóng)民工數(shù)量接近。。農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距作用機(jī)制的估計結(jié)果如表5所示。

表5 作用機(jī)制分析

第一，分析農(nóng)村勞動力轉(zhuǎn)移對城鄉(xiāng)收入差距的影響。由表5報告結(jié)果可知，農(nóng)村勞動力轉(zhuǎn)移對城鄉(xiāng)收入差距的影響為正，并且通過了顯著性水平檢驗；農(nóng)村勞動力轉(zhuǎn)移滯后項對城鄉(xiāng)收入差距的影響為正，但沒有通過顯著性水平檢驗。說明農(nóng)村勞動力轉(zhuǎn)移能夠提高農(nóng)民收入進(jìn)而縮小城鄉(xiāng)收入差距[24]，但農(nóng)村勞動力轉(zhuǎn)移對城鄉(xiāng)收入差距的影響不存在滯后效應(yīng)?？赡艿慕忉屖?，一方面，由于農(nóng)業(yè)部門內(nèi)勞動力過剩，農(nóng)村勞動力轉(zhuǎn)移會優(yōu)化農(nóng)業(yè)部門的要素配置結(jié)構(gòu)、釋放農(nóng)業(yè)部門的結(jié)構(gòu)紅利，從而提高農(nóng)民的農(nóng)業(yè)收入，縮小城鄉(xiāng)收入差距；另一方面，農(nóng)村勞動力轉(zhuǎn)移能夠增加農(nóng)民的非農(nóng)收入，農(nóng)民非農(nóng)收入的增加在一定程度也會縮小城鄉(xiāng)收入差距。

第二，分析農(nóng)業(yè)機(jī)械化對農(nóng)村勞動力轉(zhuǎn)移的影響。由表5的估計結(jié)果可以看出，農(nóng)業(yè)機(jī)械化對農(nóng)村勞動力轉(zhuǎn)移的影響為正，并且通過了顯著性水平檢驗；農(nóng)業(yè)機(jī)械化滯后項對農(nóng)村勞動力轉(zhuǎn)移的影響為正，但沒有通過顯著性水平檢驗。結(jié)果表明農(nóng)業(yè)機(jī)械化能夠促進(jìn)農(nóng)村勞動力轉(zhuǎn)移[11]，但農(nóng)業(yè)機(jī)械化對農(nóng)村勞動力轉(zhuǎn)移的影響不存在滯后效應(yīng)。可能的解釋是，由邊際技術(shù)替代效應(yīng)可知，資本要素與勞動要素之間存在著替代關(guān)系，而農(nóng)業(yè)機(jī)械化水平提高意味在農(nóng)業(yè)內(nèi)資本要素投入的增加。由勞動要素與資本要素的替代關(guān)系可知，農(nóng)業(yè)部門在一定程度上會減少勞動要素的投入量，進(jìn)而促進(jìn)農(nóng)村勞動力向非農(nóng)部門轉(zhuǎn)移。

八、穩(wěn)健性檢驗

前文的分析表明，推動農(nóng)業(yè)機(jī)械化在一定程度上能夠縮小城鄉(xiāng)收入差距，但更重要的是需要進(jìn)一步分析農(nóng)業(yè)機(jī)械化水平提升對城鄉(xiāng)收入差距的影響是否穩(wěn)健。本部分從以下四個方面對前文的分析結(jié)果進(jìn)一步討論：

第一，重新構(gòu)建農(nóng)業(yè)機(jī)械化指標(biāo)。由于農(nóng)業(yè)機(jī)械化水平和城鄉(xiāng)收入差距均為結(jié)果變量，因此模型內(nèi)可能存在內(nèi)生性問題，本部分將縣域農(nóng)業(yè)機(jī)械化總動力取自然對數(shù)作為縣域農(nóng)業(yè)機(jī)械化的衡量指標(biāo)，進(jìn)一步討論農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距是否存在影響。空間計量模型的鄰接關(guān)系、模型選取和樣本分類方式等均不變，模型的估計結(jié)果如表6所示。

表6 重新構(gòu)建農(nóng)業(yè)機(jī)械化指標(biāo)

表6報告了重新構(gòu)建農(nóng)業(yè)機(jī)械化指標(biāo)的估計結(jié)果。從表6的估計結(jié)果可以看出，農(nóng)業(yè)機(jī)械化對全國、東部縣域和西部縣域農(nóng)民收入和城鄉(xiāng)收入差距的影響為正，并且通過了顯著性水平檢驗；對中部農(nóng)民收入和城鄉(xiāng)收入差距的影響為正，但沒有通過顯著性水平檢驗。與表3的基準(zhǔn)回歸不同的是，西部縣域重新構(gòu)建了農(nóng)業(yè)機(jī)械化的變量后，農(nóng)業(yè)機(jī)械化水平提升促進(jìn)了西部縣域農(nóng)民收入水平的提升，同時也縮小了與城市居民之間的收入差距。同時空間誤差項、空間自回歸和擾動項的系數(shù)也與前文的估計結(jié)果基本一致。說明重新構(gòu)建農(nóng)業(yè)機(jī)械化指標(biāo)后，農(nóng)業(yè)機(jī)械化水平提升對城鄉(xiāng)收入差距影響的估計結(jié)果穩(wěn)健。

第二，重新定義縣域之間的相鄰關(guān)系。前文的分析過程中，空間權(quán)重矩陣的設(shè)定是基于縣域之間的距離。由地理經(jīng)濟(jì)學(xué)的相關(guān)理論可知地理位置越近的縣域，其交流也就越方便。因此，本部分借鑒孫學(xué)濤等(2018)[18]的研究成果，嘗試采用簡單權(quán)重矩陣，即兩個縣域之間是否存在空間上的相鄰，如果兩個縣域存在空間上的相鄰，權(quán)重設(shè)定為1，反之設(shè)定為0。空間計量模型的農(nóng)業(yè)機(jī)械化指標(biāo)、模型選取和樣本分類方式等均不變，模型的估計結(jié)果如表7所示。

表7 重新定義相鄰關(guān)系

從表7的估計結(jié)果可以看出，農(nóng)業(yè)機(jī)械化對全國和東部縣域農(nóng)民收入和城鄉(xiāng)收入差距的影響為正，并且通過了顯著性水平檢驗；對中部縣域農(nóng)民收入和城鄉(xiāng)收入差距的影響為負(fù)，但沒有通過顯著性水平檢驗；對西部縣域農(nóng)民收入和城鄉(xiāng)收入差距的影響顯著為負(fù)，并且通過了顯著性水平檢驗?？臻g誤差項、空間自回歸和擾動項的系數(shù)也與前文的估計結(jié)果基本一致。上述結(jié)果表明重新定義空間計量模型相鄰關(guān)系后，農(nóng)業(yè)機(jī)械化水平提升對農(nóng)民收入和城鄉(xiāng)收入差距的估計結(jié)果穩(wěn)健。

第三，重新構(gòu)建空間計量模型。前文的分析過程中空間計量模型采用SARAR模型分析農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響，在該模型內(nèi)空間自回歸項和空間誤差項同時存在，本部分嘗試采用空間自回歸模型和空間誤差模型分別研究農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的影響?？臻g計量模型的農(nóng)業(yè)機(jī)械化指標(biāo)、鄰接關(guān)系和樣本分類方式等均不變，模型的估計結(jié)果如表8所示。

表8 重新構(gòu)建空間計量模型

從表8的估計結(jié)果可以看出，農(nóng)業(yè)機(jī)械化水平對全國和東部縣域農(nóng)民收入和城鄉(xiāng)收入差距的影響為正，并且通過了顯著性水平檢驗；對中部縣域農(nóng)民收入和城鄉(xiāng)收入差距的影響為正，但沒有通過顯著性水平檢驗；對西部縣域農(nóng)民收入的影響為負(fù)，但對西部縣域城鄉(xiāng)收入差距的影響為正，并且通過了顯著性水平檢驗。同時空間誤差項和擾動項的系數(shù)也與前文基本一致。即重新構(gòu)建空間計量模型后，農(nóng)業(yè)機(jī)械化水平提升對農(nóng)民收入和城鄉(xiāng)收入差距影響的估計結(jié)果穩(wěn)健。

第四，將樣本進(jìn)一步細(xì)分。前文的分析過程中樣本的數(shù)據(jù)采用全國縣域的數(shù)據(jù)以及將全國縣域數(shù)據(jù)按照東、中、西部劃分為三個地區(qū)。本部分嘗試借鑒國務(wù)院發(fā)展研究中心的區(qū)劃方法，將縣域劃分為八大經(jīng)濟(jì)區(qū)，進(jìn)一步討論農(nóng)業(yè)機(jī)械化對農(nóng)民收入和城鄉(xiāng)收入差距的影響?？臻g計量模型的農(nóng)業(yè)機(jī)械化指標(biāo)、鄰接關(guān)系和模型選取等均不變，模型的估計結(jié)果如表9所示。

表9 進(jìn)一步細(xì)分樣本

限于文章的篇幅，本文只報告了農(nóng)業(yè)機(jī)械化對農(nóng)民收入的估計結(jié)果，農(nóng)業(yè)機(jī)械化對城鄉(xiāng)收入差距的估計未報告。從表9的估計結(jié)果可以看出，農(nóng)業(yè)機(jī)械化水平對南部沿海和東北綜合經(jīng)濟(jì)區(qū)農(nóng)民收入的影響為正，并且通過了顯著性水平檢驗；對北部沿海、大西北綜合經(jīng)濟(jì)區(qū)農(nóng)民收入的影響顯著為負(fù)；對其余經(jīng)濟(jì)區(qū)農(nóng)民收入的影響沒有通過顯著性水平檢驗。由東、中、西部地區(qū)與八大經(jīng)濟(jì)區(qū)的對應(yīng)關(guān)系可知，農(nóng)業(yè)機(jī)械化對農(nóng)民收入的影響與前文的估計結(jié)果基本一致?？臻g誤差項和擾動項的系數(shù)與前文的估計結(jié)果也基本一致。這說明將縣域細(xì)分為八大經(jīng)濟(jì)區(qū)后，農(nóng)業(yè)機(jī)械化水平提升對農(nóng)民收入的估計結(jié)果穩(wěn)健。

九、結(jié)論與政策建議

研究發(fā)現(xiàn)，中國縣域農(nóng)業(yè)機(jī)械化水平存在著顯著的空間負(fù)相關(guān)關(guān)系，即本縣域的農(nóng)業(yè)機(jī)械化水平提升會抑制相鄰縣域農(nóng)業(yè)機(jī)械化水平的提升；農(nóng)業(yè)機(jī)械化能夠提高農(nóng)民收入，最終實現(xiàn)縮小城鄉(xiāng)收入差距的目標(biāo)；農(nóng)業(yè)機(jī)械化水平提升能夠提高東部地區(qū)農(nóng)民收入，進(jìn)而縮小其城鄉(xiāng)收入差距，但農(nóng)業(yè)機(jī)械化水平提升會抑制西部地區(qū)的農(nóng)民增收，進(jìn)而導(dǎo)致其收入差距加大；農(nóng)業(yè)機(jī)械化不僅會直接提高農(nóng)民收入，進(jìn)而縮小城鄉(xiāng)收入差距，而且農(nóng)業(yè)機(jī)械化還會替代農(nóng)業(yè)勞動要素，促進(jìn)農(nóng)村勞動力轉(zhuǎn)移，最終實現(xiàn)縮小收入差距的目標(biāo)。

根據(jù)研究結(jié)論，為促進(jìn)農(nóng)民增收和縮小城鄉(xiāng)收入差距，本文提出如下政策建議：

第一，推動農(nóng)業(yè)機(jī)械化，提高農(nóng)民收入。農(nóng)業(yè)機(jī)械化水平提升不僅能夠提高農(nóng)民收入，而且還會促進(jìn)城鄉(xiāng)收入差距的縮小。因此，在人口紅利即將消失和農(nóng)村勞動力轉(zhuǎn)移的背景下，農(nóng)業(yè)機(jī)械化成為保證糧食安全和促進(jìn)農(nóng)民增收的重要方式。在適合農(nóng)業(yè)機(jī)械的地區(qū)推進(jìn)農(nóng)業(yè)機(jī)械化，同時鼓勵農(nóng)業(yè)機(jī)械跨區(qū)作業(yè)，形成有中國特色的農(nóng)業(yè)機(jī)械化模式，進(jìn)而將農(nóng)業(yè)機(jī)械化與小農(nóng)經(jīng)營相耦合，促進(jìn)農(nóng)民增收，最終實現(xiàn)縮小收入差距的目標(biāo)。

第二，推動農(nóng)業(yè)機(jī)械化，促進(jìn)農(nóng)村勞動力轉(zhuǎn)移。由邊際技術(shù)替代效應(yīng)和前文結(jié)論可知，促進(jìn)農(nóng)民增收和縮小城鄉(xiāng)收入差距的根本出路在于推動農(nóng)村勞動力轉(zhuǎn)移，只有改善了農(nóng)業(yè)部門的勞動要素和土地要素的比例，進(jìn)而通過提高農(nóng)業(yè)機(jī)械化水平促進(jìn)農(nóng)民增收和縮小城鄉(xiāng)收入差距才能成為可能。但需要處理好農(nóng)業(yè)機(jī)械化在農(nóng)村勞動力轉(zhuǎn)移過程中的作用，即農(nóng)業(yè)機(jī)械對農(nóng)村勞動力的替代是以農(nóng)業(yè)生產(chǎn)未受到損失為前提，而農(nóng)村勞動力轉(zhuǎn)移并不能通過農(nóng)業(yè)機(jī)械化促進(jìn)農(nóng)民增收和縮小城鄉(xiāng)收入差距。

第三，充分利用縣域資源，推動農(nóng)業(yè)機(jī)械化發(fā)展。由于中國縣域經(jīng)濟(jì)發(fā)展水平、自然資源環(huán)境和科學(xué)技術(shù)水平等方面存在較大的差異，因此中國縣域經(jīng)濟(jì)之間需要根據(jù)自身的比較優(yōu)勢，推動農(nóng)業(yè)機(jī)械化。東部地區(qū)應(yīng)該充分發(fā)揮資本和科技優(yōu)勢，走資本或技術(shù)密集型的農(nóng)業(yè)機(jī)械化道路；中、西部地區(qū)應(yīng)該根據(jù)本縣域的地形和地勢等資源推動農(nóng)業(yè)機(jī)械化，即地形復(fù)雜和地勢起伏相對較大的地區(qū)實施退耕還林還草，在地形平坦和地勢起伏相對較小的地區(qū)推進(jìn)農(nóng)業(yè)機(jī)械化。

第四，加強縣域之間的交流與合作?？h域農(nóng)業(yè)機(jī)械化存在顯著的負(fù)向空間相關(guān)關(guān)系，即農(nóng)業(yè)機(jī)械化水平較高的縣域與農(nóng)業(yè)機(jī)械化水平較低的縣域相鄰。農(nóng)業(yè)機(jī)械水平較低的縣域可以借鑒農(nóng)業(yè)機(jī)械化水平較高縣域的先進(jìn)技術(shù)和管理經(jīng)驗來提升本地區(qū)的農(nóng)業(yè)機(jī)械化水平，從而縮小縣域之間農(nóng)業(yè)機(jī)械化水平的差距，進(jìn)一步促進(jìn)縣域農(nóng)業(yè)機(jī)械化均衡發(fā)展。