□文/王懷亮

（菏澤學(xué)院 山東·菏澤）

引言

在全面建設(shè)小康社會(huì)和實(shí)現(xiàn)中國夢背景下，農(nóng)民、農(nóng)村和農(nóng)業(yè)的三農(nóng)問題已成為政府和社會(huì)廣泛關(guān)注的社會(huì)焦點(diǎn)問題之一。從農(nóng)民角度看，三農(nóng)問題的核心是農(nóng)民收入問題。山東省是一個(gè)農(nóng)業(yè)大省，農(nóng)民收入問題是政府和各級學(xué)者高度關(guān)注和分析的問題。改革開放以來，隨著山東經(jīng)濟(jì)水平的快速發(fā)展，山東農(nóng)民收入也逐漸增加，但由于各地市資源條件、交通條件、經(jīng)濟(jì)發(fā)展水平和農(nóng)民自身素養(yǎng)的不同，各地市農(nóng)民收入增長呈現(xiàn)出顯著的差異性，這種不平衡已成為政府關(guān)注的熱點(diǎn)。然而，從空間的視角研究山東農(nóng)民收入分布特征及其空間關(guān)聯(lián)的文獻(xiàn)缺乏。本研究以空間視角選取山東省2013年的17個(gè)地區(qū)農(nóng)民人均純收入的空間差異，以探討分析山東17地區(qū)的農(nóng)民人均純收入分布情況，為相關(guān)的管理部門制定山東三農(nóng)問題政策提供科學(xué)支持。

一、研究方法與數(shù)據(jù)來源

（一）研究方法?？臻g自相關(guān)反映的是一個(gè)區(qū)域單元上的某種地理屬性與鄰近區(qū)域單元上的同一屬性值的相關(guān)程度，有全局指標(biāo)和局部指標(biāo)測度，主要包括空間相互作用與空間擴(kuò)散過程。本研究將基于R語言使用空間統(tǒng)計(jì)方法中的全局統(tǒng)計(jì)量Moran’s系數(shù)、Geary’sC比率和局部空間數(shù)據(jù)插值進(jìn)行山東省各市農(nóng)民人均純收入空間統(tǒng)計(jì)自相關(guān)分析。首先，對文中涉及的幾個(gè)指標(biāo)做簡要介紹。

表1 2013年山東各市區(qū)農(nóng)民人均純收入

1、全局指標(biāo)。全局Moran’sI系數(shù)反映空間鄰接區(qū)域的屬性值的近似程度，其值取值范圍是［-1，1］，0值表示所研究屬性在空間分布為隨機(jī)模式，-1表示反研究屬性在空間分布存在完全的負(fù)空間自相關(guān)，1表示所研究屬性在空間分布存在完全的正空間自相關(guān)，絕對值越大表示相關(guān)程序越高。Geary’sC的取值范圍為［0，2］，1值表示所研究屬性在空間分布為隨機(jī)模式，當(dāng)Geary’sC取值為0～1時(shí)，表示所研究屬性存在正的空間自相關(guān)，當(dāng)Geary’sC取值為1～2時(shí)，表示所研究屬性存在負(fù)的空間自相關(guān)。R包spdep可以計(jì)算這兩個(gè)全局指標(biāo)。

2、局部指標(biāo)。局部空間自相關(guān)統(tǒng)計(jì)分析方法常用空Moran散點(diǎn)圖和空間數(shù)據(jù)插值分析等方法。Moran散點(diǎn)圖四個(gè)象限分別對應(yīng)于區(qū)域單元與其鄰近單元之間四種類型的局部空間關(guān)系：第一象限表示高觀測值的區(qū)域單元被高值的區(qū)域包圍的空間分布關(guān)系；第二象限表示低觀測值的區(qū)域單元被高值的區(qū)域包圍的空間分布關(guān)系；第三象限表示高低測值的區(qū)域單元被低值的區(qū)域包圍的空間分布關(guān)系；第四象限表示高測值的區(qū)域單元被低值的區(qū)域包圍的空間分布關(guān)系；此局部指標(biāo)比局部Moran統(tǒng)計(jì)量系數(shù)能夠進(jìn)一步指出研究單元和其鄰近單元屬于高值和高值、低值和低值、高值和低值、低值和高值之中的哪種空間聯(lián)系方式。

（二）數(shù)據(jù)來源與處理。本研究以山東省地市級行政區(qū)劃為準(zhǔn)，總共包括17個(gè)地市級單位，各地級市的研究中心點(diǎn)假定在各地級市的市府所在經(jīng)緯度點(diǎn)，農(nóng)民人均收入統(tǒng)計(jì)數(shù)據(jù)來自于2014年山東省統(tǒng)計(jì)年鑒，數(shù)據(jù)如表1。指標(biāo)計(jì)算由R軟件中的包spdep計(jì)算得到，空間插值圖利用arcgis10.2分析計(jì)算。（表1）

表2 山東省2013年市域農(nóng)民人均純收入全局Moran’sI系數(shù)及檢驗(yàn)

表3 山東省2013年市域農(nóng)民人均純收入全局Geary’sC系數(shù)及檢驗(yàn)

圖1 2013年山東各地級市區(qū)農(nóng)民人均純收入空間插值分布圖

二、山東省地級市農(nóng)民人均純收入空間自相關(guān)分析

（一）全局空間自相關(guān)分析。利用全局Moran’sI系數(shù)和Geary’sC比率統(tǒng)計(jì)量，考察山東農(nóng)民人均純收入在整個(gè)研究區(qū)域內(nèi)空間相關(guān)的整體趨勢。在Moran’sI系數(shù)與Geary’sC比率檢驗(yàn)中，采用基于研究數(shù)據(jù)隨機(jī)檢驗(yàn)性和正態(tài)近似性檢驗(yàn)兩種方法，保證結(jié)果的可靠性。基于R包spdep計(jì)算出的Moran’sI系數(shù)、Geary’sC比率與檢驗(yàn)結(jié)果見表2和表3。（表2、表3）

從表 2 可知，Moran’sI系數(shù)為 0.347，表 3Geary’sC 系數(shù)為0.518，這說明從山東省空間全局來看，2013年山東省各個(gè)地級市的農(nóng)民人均純收入存在著一定的空間正相關(guān)，即農(nóng)民人均純收入高的地區(qū)與農(nóng)民人均純收入低的地區(qū)在山東省分布存在著一定的空間聚集，較高農(nóng)民人均純收入水平的地級市相對地趨于和較高農(nóng)民人均純收入水平的地級市相鄰；較低農(nóng)民人均純收入水平的地級市相對地趨于和較低農(nóng)民人均純收入水平的地級市相鄰?；?種檢驗(yàn)的條件下，Moran’sI與Geary’sC的統(tǒng)計(jì)量檢驗(yàn)的p值都小于0.05，說明全局自相關(guān)Moran’sI系數(shù)、Geary’sC系數(shù)都是顯著的，即所求的全局自相關(guān)Moran’sI系數(shù)、Geary’sC系數(shù)都是有意義的。

（二）局部空間自相關(guān)分析?？臻g全局統(tǒng)計(jì)量僅揭示了研究區(qū)域內(nèi)的總體情況，為了進(jìn)一步揭示山東省市域農(nóng)民人均純收入在局部尺度上的空間集聚程度，利用空間數(shù)據(jù)插值來探測和可視化局部的空間自相關(guān)和空間差異，如圖1所示。（圖1）可以看出，2013年山東各地級市農(nóng)民人均純收入確實(shí)存在空間正相關(guān)分布。具體來說，山東的膠東半島農(nóng)民人均純收入高，存在高值聚集，有著很好的帶動(dòng)輻射作用；魯西南地區(qū)農(nóng)民人均純收入低，存在著低值聚集。

三、結(jié)論

本文對2013山東省農(nóng)民人均純收入在地市級空間尺度上進(jìn)行了空間自相關(guān)分析，綜合上述研究得到如下主要結(jié)論：

（一）對區(qū)域經(jīng)濟(jì)差異的研究，普通統(tǒng)計(jì)分析方法，忽視了空間自相關(guān)性，與經(jīng)濟(jì)發(fā)展現(xiàn)實(shí)不符，通過引入空間地理屬性數(shù)據(jù)，解釋經(jīng)濟(jì)數(shù)據(jù)空間分布規(guī)律，增強(qiáng)了分析結(jié)果的可靠性。

（二）全局分析表明，2013年山東農(nóng)民人均純收入始終呈現(xiàn)出相似值間的空間集聚。局部分析表明魯西南地區(qū)農(nóng)民人均純收入較低，整個(gè)山東省各個(gè)地級市在農(nóng)民人均純收入水平上存在著失衡現(xiàn)象，這種現(xiàn)象不僅體現(xiàn)在整個(gè)山東省的農(nóng)民人均純收入的分布上，也體現(xiàn)在局部的地級市的農(nóng)民人均純收入水平分布上，所以在山東十三五規(guī)劃上，要繼續(xù)大力發(fā)展魯西南地區(qū)的經(jīng)濟(jì)水平，特別是菏澤地區(qū)，以增加當(dāng)?shù)氐霓r(nóng)民人均純收入，從而推動(dòng)整個(gè)山東省農(nóng)村經(jīng)濟(jì)的健康、和諧、持續(xù)發(fā)展，早日實(shí)現(xiàn)山東農(nóng)民的中國夢。

［1］楊中慶.基于R語言的空間統(tǒng)計(jì)分析研究與應(yīng)用.暨南大學(xué)，2006.

［2］張海軍.河南省市域農(nóng)民人均收入空間自相關(guān)分析［J］.南陽師范學(xué)院學(xué)報(bào)，2009.9.

［3］池建.精通ArcGIS地理信息系統(tǒng)［M］.北京：清華大學(xué)出版社，2011.