王振波,王新明

1 中國科學(xué)院地理科學(xué)與資源研究所,北京 1001012 中國科學(xué)院大學(xué),北京 100049

自2007年7月1日實施《中華人民共和國農(nóng)民專業(yè)合作社法》以來,中國農(nóng)民專業(yè)合作社進入了依法發(fā)展的新時期,“走到了一個新的十字路口”,從橫向合作向縱向合作深化,從單一功能向多種功能拓展,從傳統(tǒng)合作向新型合作演變,從農(nóng)戶間合作向社際協(xié)作邁進,中國農(nóng)民專業(yè)合作社取得了引人矚目的發(fā)展,其發(fā)展形式更加組織化、產(chǎn)業(yè)化與綜合化[1-2]。2013年中央1號文件首次提出“農(nóng)民合作社”的概念,這是根據(jù)近年來我國農(nóng)民合作組織蓬勃發(fā)展的實際情況提出的,農(nóng)民專業(yè)合作社涉及的領(lǐng)域越來越多,由種植養(yǎng)殖業(yè)向加工、勞務(wù)、運輸、信息、資金、技術(shù)和銷售等眾多領(lǐng)域延伸。隨著多種類型合作社的出現(xiàn),單純的“專業(yè)合作”已經(jīng)不能涵蓋農(nóng)民合作的內(nèi)容了,“農(nóng)民合作社”這個概念更科學(xué)、更全面,除了單純從事種養(yǎng)業(yè)外,農(nóng)機服務(wù)合作社、鄉(xiāng)村旅游合作社、手工業(yè)合作社等從事二三產(chǎn)業(yè)的合作社快速發(fā)展,資金互助合作社、土地合作社、社區(qū)股份合作社、富民合作社等提供金融、土地、租賃等創(chuàng)新服務(wù)的新型合作社也方興未艾[3]。截止到2017年7月底,全國依法登記的農(nóng)民專業(yè)合作社達193.3萬家,是2007年底農(nóng)民專業(yè)合作社法頒布施行初期的74倍,年均增長60%。實有入社農(nóng)戶超過1億戶,約占全國農(nóng)戶總數(shù)的46.8%,平均每個行政村有近3家合作社[4]。由中華人民共和國農(nóng)業(yè)部、國家發(fā)展改革委、財政部、國土資源部、人民銀行、稅務(wù)總局于2017年聯(lián)合印發(fā)的《關(guān)于促進農(nóng)業(yè)產(chǎn)業(yè)化聯(lián)合體發(fā)展的指導(dǎo)意見》中明確指出組建農(nóng)業(yè)產(chǎn)業(yè)化聯(lián)合體要著力提升農(nóng)民合作社服務(wù)能力[5]；“十三五”脫貧攻堅規(guī)劃亦提出要培育壯大貧困地區(qū)農(nóng)民專業(yè)合作社這一新型經(jīng)營主體。農(nóng)民合作社作為農(nóng)村經(jīng)濟發(fā)展中的重要組織載體和經(jīng)濟力量之一,將以新型農(nóng)業(yè)經(jīng)營主體的身份成為大政策和行動的承載主體、實施主體,成為推進鄉(xiāng)村振興、精準扶貧戰(zhàn)略以及實現(xiàn)小農(nóng)戶與農(nóng)業(yè)現(xiàn)代化發(fā)展有機銜接的重要力量[6-7]。

基于中國知網(wǎng)檢索(設(shè)置篇名為“農(nóng)業(yè)合作社”或含“農(nóng)民專業(yè)合作社”或含“農(nóng)民合作社”,期刊為CSSCI和CSCD)發(fā)現(xiàn)共有關(guān)于此主題的論文874篇,從合作社法實施年份開始論文發(fā)表量呈波動增加趨勢,至2014年達到峰值126篇(圖1)。基于Web of Science Core Cllection用同樣的檢索策略進行搜索,結(jié)果呈現(xiàn)29篇文獻。對以往研究分析發(fā)現(xiàn),目前對農(nóng)民合作社的主要集中于其立法探討[8-10]、組織發(fā)展[11-12]、財務(wù)管理[13]、發(fā)展模式[14]、建設(shè)治理[15-17]、制度安排[18-19]、融資困境與解決之道[20- 22]、供應(yīng)鏈管理[23- 24]、對合作社發(fā)展過程中呈現(xiàn)問題的反思以及發(fā)展方向的思考[25- 29]、合作社發(fā)展與精準扶貧[30-31]、合作社發(fā)展的影響因素分析[32-34]、合作社效率或績效分析[35-39]、對聯(lián)合合作社的探討[40]、國際經(jīng)驗的梳理[41- 43]等。研究方法主要有案例分析[44]、實地調(diào)研[45]、計量模型[46]、問卷調(diào)查[47]等。研究尺度有省域[48-49]、村鎮(zhèn)或具體合作社[50-52]、或是某一類合作社的研究[53]等。綜上,農(nóng)民專業(yè)合作社的相關(guān)研究成果是理論研究多于實證研究,定性分析多于定量分析,微觀研究多于宏觀研究,基本是小尺度研究,缺乏某一區(qū)域合作社分布時空動態(tài)特征的研究資料,在進行合作社發(fā)展影響因素分析時更是少有資料將時間與空間因素納入回歸模型。此外,京津冀協(xié)同發(fā)展為合作社發(fā)展帶來了新機遇,協(xié)同發(fā)展有利于充分發(fā)揮河北的農(nóng)業(yè)資源優(yōu)勢,從而形成冀與京津的供需暢通格局。新型城鎮(zhèn)化、產(chǎn)業(yè)結(jié)構(gòu)升級和產(chǎn)業(yè)布局調(diào)整,新修訂的《農(nóng)民專業(yè)合作社法》增設(shè)“農(nóng)民專業(yè)社聯(lián)合社”一章,這些都為京津冀現(xiàn)代農(nóng)業(yè)發(fā)展提供了更為廣闊的空間。因此將“過程-格局-機理”的地理研究思維作為本文的分析框架,不僅可以填補以上研究空白,而且可以充分發(fā)揮地理學(xué)科的綜合性與交叉性優(yōu)勢,促進地理學(xué)科與其他學(xué)科的融合發(fā)展,同時,也可以為京津冀農(nóng)民合作社向綠色生態(tài)產(chǎn)業(yè)合作社轉(zhuǎn)型提供思路、貢獻力量。

圖1 中國知網(wǎng)有關(guān)農(nóng)民專業(yè)合作社的論文發(fā)表量(CSSCI與CSCD：1999—2018)Fig.1 The amount of papers published on China Knowledge Network about farmer professional cooperatives (CSSCI and CSCD: 1999—2018)

農(nóng)民專業(yè)合作社與農(nóng)民合作社是不同時期對同一經(jīng)營組織的稱呼,本質(zhì)上無差別,由于本文的研究年限為2007—2013年,故仍采用舊稱農(nóng)民專業(yè)合作社以迎合時代特征。本文將京津冀13個地市的農(nóng)民專業(yè)合作社依據(jù)合作社名稱分為種植業(yè)、養(yǎng)殖業(yè)、供產(chǎn)銷與農(nóng)機專業(yè)合作社四類,之后基于ArcGIS平臺,先進行可視化表達、統(tǒng)計分析與空間自相關(guān)性分析,發(fā)現(xiàn)其時空演化規(guī)律,再運用空間計量模型進行影響因素的解析,最后得出結(jié)論并指出不足。

1 研究區(qū)域、數(shù)據(jù)說明及研究方法

1.1 研究區(qū)域概況

圖2 京津冀區(qū)域范圍Fig.2 The region of Beijing-Tianjin-Hebei

京津冀(Beijing-Tianjin-Hebei Urban Agglomeration,BTHUA) 是中國的“首都經(jīng)濟圈”,五大國家級城市群之一,包括北京、天津和河北省石家莊、保定、滄州、承德、邯鄲、衡水、廊坊、秦皇島、唐山、邢臺、張家口等13個地級以上城市,總面積21.72 km2,占全國總面積的2.26%(圖2)。京津地區(qū)人多地少,耕地資源嚴重匱乏,相當部分的農(nóng)產(chǎn)品靠外地輸入,而區(qū)域內(nèi)河北地區(qū)可耕地面積達600 hm2,居全國第4位,農(nóng)業(yè)資源豐富,勞動力價格低廉,具備優(yōu)越的發(fā)展基礎(chǔ)。

1.2 數(shù)據(jù)來源及處理

本研究的數(shù)據(jù)主要分為兩部分,一是每個農(nóng)民專業(yè)合作社地理坐標的獲取,二是京津冀各市各年份各類農(nóng)民專業(yè)合作社的數(shù)量及影響因素指標數(shù)據(jù)的收集。為了進行GIS可視化表達,需要獲取各類農(nóng)民專業(yè)合作社的經(jīng)緯度,該工作是借助地理空間數(shù)據(jù)分析工具Python3.6編寫經(jīng)緯度拾取程序代碼從百度地圖上獲取。各類農(nóng)民專業(yè)合作社的數(shù)量借助ArcGIS10.2操作平臺的“相交”工具統(tǒng)計；影響因素指標包括人均GDP、職工平均工資、第一產(chǎn)業(yè)從業(yè)人員、第一產(chǎn)業(yè)產(chǎn)值、農(nóng)林水務(wù)財政支出,其中前4個指標來源于《城市統(tǒng)計年鑒》,農(nóng)林水務(wù)財政支出來源于各市統(tǒng)計年鑒。為消除價格因素影響,以2007年為基期采用GDP平減指數(shù)法對人均GDP、職工平均工資、第一產(chǎn)業(yè)產(chǎn)值作預(yù)處理。由于2007年《中華人民共和國農(nóng)民專業(yè)合作社法》開始實施,2013年中央1號文件首次提出“農(nóng)民合作社”的概念,故選取兩個具有代表性、轉(zhuǎn)折性的時間點為研究年限(2007—2013),以期為農(nóng)民專業(yè)合作社提供階段性成果資料。

1.3 研究方法

1.3.1 空間自相關(guān)性分析

空間自相關(guān)性檢驗通常有Moran′s I 指數(shù)和Geary′s C 指數(shù)兩種方法,本文采用Moran′s I指數(shù)測度空間關(guān)聯(lián)性,可以有效的檢驗區(qū)域內(nèi)單元的空間依賴性、空間異質(zhì)性或者相互獨立性[54-55]。其原理如下：

(1)

(2)

式中,Z為檢E(I)驗統(tǒng)計量,為期望值,Var(I)為方差。

1.3.2 空間計量模型

空間計量模型可以有效解決普通回歸分析無法處理空間依賴性的問題,常用的回歸模型包括空間誤差(SEM)、空間滯后(SLM)等[56-58]。當模型的誤差項存在空間相關(guān)時,即為SEM模型,可表達如下：

Y=Xβ+λWμ+ε,ε∈N(0,δ2)

(3)

式中,Y表示被解釋變量,為(n×1)的向量；X表示解釋變量,假設(shè)解釋變量共K個,則為(n×k)的矩陣；β表示回歸系數(shù),為(k×1)向量；μ為隨機誤差向量；λ為回歸殘差之間的空間相關(guān)性系數(shù)；W表示(n×n)空間權(quán)重矩陣,ε是獨立分布的隨機項。

當被解釋變量間的空間依賴性顯著且會影響模型效果時,即為空間自回歸模型,也稱作空間滯后模型,可表達如下：

Y=ρWY+Xβ+ε,ε∈N(0,δ2)

(4)

式中,ρ表示內(nèi)生交互效應(yīng)(WY)的系數(shù),其大小反應(yīng)空間擴散與空間依賴的程度,如其顯著,表示被解釋變量存在一定的空間依賴性。

2 京津冀農(nóng)民專業(yè)合作社時空演化

依據(jù)其詳細名稱將農(nóng)民專業(yè)合作社分為種植業(yè)、養(yǎng)殖業(yè)、供銷、農(nóng)機專業(yè)合作社共4類,利用ArcGIS平臺對2007—2013年京津冀各類農(nóng)民專業(yè)合作社的空間位置分別進行可視化表達,以揭示其時空演化特征。

2.1 總體特征

四類農(nóng)民專業(yè)合作社數(shù)量總體上呈逐年增長趨勢,且合作社數(shù)量為種植業(yè)>養(yǎng)殖業(yè)>供銷>農(nóng)機,增長速度亦是如此(圖3)。那是因為種植業(yè)與養(yǎng)殖業(yè)受資源約束相對較弱,京津冀乃至中國大部分農(nóng)戶都從事種植業(yè)、養(yǎng)殖業(yè)生產(chǎn),此類產(chǎn)業(yè)的從業(yè)人員較多,因此較為容易發(fā)展。相對而言,供銷與農(nóng)機類需要一定的初始資本與專業(yè)管理、技術(shù)人力資源積累,因此制約了此類農(nóng)民專業(yè)合作社的發(fā)展；加之農(nóng)機專業(yè)合作社屬于服務(wù)性合作經(jīng)濟組織,其服務(wù)的對象和服務(wù)的范圍也具有一定的區(qū)域特征(如對于土地的要求應(yīng)該是平整連片等),進一步限制了農(nóng)機類合作社的發(fā)展。另外市場需求影響供給,京津冀乃至全國的消費者對于必需品種植類產(chǎn)品(如糧食、水果、蔬菜等)的需求大于養(yǎng)殖類產(chǎn)品(如生豬、肉牛羊等)及其他類產(chǎn)品。

圖3 京津冀各類農(nóng)民專業(yè)合作社總數(shù)(2007—2013)Fig.3 Total amount of farmers′ professional cooperatives in Beijing-Tianjin-Hebei (2007—2013)

2.2 時空演化

2.2.1 種植業(yè)專業(yè)合作社時空演化

京津冀種植業(yè)專業(yè)合作社數(shù)量整體上呈現(xiàn)逐年增加趨勢,增加速度為北京>天津>石家莊>其他市(圖4)。北京、天津與石家莊的種植業(yè)專業(yè)合作社有明顯的集聚分布效應(yīng),雙核特征顯著,經(jīng)濟發(fā)展水平上呈現(xiàn)的“省會城市一家獨大”現(xiàn)象在農(nóng)業(yè)領(lǐng)域亦有所體現(xiàn),隨著時間遷移,合作社分布由兩核向東南部擴散,秦皇島、唐山、滄州、衡水、邢臺、邯鄲組成東南發(fā)展帶,聯(lián)通兩核,因此京津冀種植業(yè)專業(yè)合作社呈現(xiàn)“兩核一帶”空間格局,合作社數(shù)量東南部多于西北部,廊坊、保定夾在“兩核一帶”中間形成中部塌陷區(qū)(圖5)。

圖4 京津冀各市種植業(yè)專業(yè)合作社數(shù)量(2007—2013)Fig.4 Number of planting industry professional cooperatives in each city of Beijing-Tianjin- Hebei (2007—2013)

圖5 京津冀種植業(yè)專業(yè)合作社時空演化(2007—2013)Fig.5 Time and space evolution of planting industry professional cooperatives in Beijing-Tianjin-Hebei (2007—2013)

2.2.2 養(yǎng)殖業(yè)專業(yè)合作社時空演化

京津冀養(yǎng)殖業(yè)專業(yè)合作社數(shù)量整體上呈現(xiàn)逐年增加趨勢,增加速度為北京>天津>其他市(圖6)。北京、天津與石家莊的養(yǎng)殖業(yè)業(yè)專業(yè)合作社有明顯的集聚分布效應(yīng),雙核特征顯著,經(jīng)濟發(fā)展水平上呈現(xiàn)的“省會城市一家獨大”現(xiàn)象在農(nóng)業(yè)領(lǐng)域亦有所體現(xiàn),隨著時間遷移,合作社分布由兩核向東南部擴散,秦皇島、唐山、滄州、衡水、邢臺、邯鄲組成東南發(fā)展帶,聯(lián)通兩核,因此京津冀種植業(yè)專業(yè)合作社呈現(xiàn)“兩核一帶”空間格局,合作社數(shù)量東南部多于西北部,廊坊、保定夾在“兩核一帶”中間形成中部塌陷區(qū)。養(yǎng)殖業(yè)專業(yè)合作社空間演化特征與種植業(yè)專業(yè)合作社近似(圖7)。

圖6 京津冀各市養(yǎng)殖業(yè)專業(yè)合作社數(shù)量(2007—2013)Fig.6 Number of aquaculture industry professional cooperatives in each city of Beijing-Tianjin- Hebei (2007—2013)

圖7 京津冀養(yǎng)殖業(yè)專業(yè)合作社時空演化(2007—2013)Fig.7 Time and space evolution of aquaculture industry professional cooperatives in Beijing-Tianjin-Hebei (2007—2013)

2.2.3 供銷專業(yè)合作社時空演化

京津冀供銷專業(yè)合作社數(shù)量整體上呈現(xiàn)逐年增加趨勢,但增長速度較緩慢(北京>天津>保定>其他市),明顯小于種植業(yè)與養(yǎng)殖業(yè)專業(yè)合作社(圖8)。京津冀供銷專業(yè)合作社分布梯度特征明顯,合作社數(shù)量表現(xiàn)為中部>南部>北部,北京、天津呈現(xiàn)明顯的極化效應(yīng)(圖9)。

圖8 京津冀各市供銷專業(yè)合作社數(shù)量(2007—2013)Fig.8 Number of supply & marketing professional cooperatives in each city of Beijing-Tianjin- Hebei (2007—2013)

圖9 京津冀供銷專業(yè)合作社時空演化(2007—2013)Fig.9 Time and space evolution of supply & marketing professional cooperatives in Beijing-Tianjin-Hebei (2007—2013)

2.2.4 農(nóng)機專業(yè)合作社時空演化

京津冀農(nóng)機專業(yè)合作社數(shù)量整體上呈現(xiàn)逐年增加趨勢,增長速度北京、天津>石家莊>其他市,北京與天津呈現(xiàn)互相追趕態(tài)勢(圖10)。京津冀農(nóng)機專業(yè)合作社分布梯度特征明顯,合作社數(shù)量表現(xiàn)為中部>南部>北部,2007—2011年北京、天津呈現(xiàn)明顯的極化效應(yīng),2012、2013年演化為北京、天津與石家莊雙核效應(yīng),2013年“雙核一帶”格局初現(xiàn),廊坊、保定夾在“兩核一帶”中間形成中部塌陷區(qū)。西、北部張家口、承德等城市的農(nóng)機專業(yè)合作社數(shù)量較低,那是因為京津冀西北向為燕山-太行山山系構(gòu)造,東、南面相對與西北面較為平坦,對于農(nóng)機專業(yè)合作社的發(fā)展具有地形上的先天優(yōu)勢(圖11)。

圖10 京津冀各市農(nóng)機專業(yè)合作社數(shù)量(2007—2013)Fig.10 Number of agricultural machinery professional cooperatives in each city of Beijing-Tianjin- Hebei (2007—2013)

圖11 京津冀農(nóng)機專業(yè)合作社時空演化(2007—2013)Fig.11 Time and space evolution of agricultural machinery professional cooperatives in Beijing-Tianjin-Hebei (2007—2013)

2.3 空間自相關(guān)性分析

空間自相關(guān)結(jié)果顯示(表1),種植業(yè)專業(yè)合作社2007與2013年的Moran′ I 指數(shù)分別為-0.062、-0.206,且均在5%的顯著性水平上通過了檢驗,表明京津冀種植業(yè)專業(yè)合作社分布呈現(xiàn)明顯的空間負相關(guān)性,且在2007—2013年間相關(guān)性呈增強趨勢；養(yǎng)殖業(yè)專業(yè)合作社2007與2013年的Moran′I 指數(shù)分別為-0.052、-0.056,且均在10%的顯著性水平上通過了檢驗,表明京津冀養(yǎng)殖業(yè)專業(yè)合作社呈現(xiàn)明顯的空間負相關(guān)性,且在2007—2013年間相關(guān)性呈增強趨勢；供銷專業(yè)合作社2007與2013年的Moran′I 指數(shù)分別為-0.120、-0.078,且分別在1%、5%的顯著性水平上通過了檢驗,表明京津冀供銷專業(yè)合作社呈現(xiàn)明顯的空間負相關(guān)性,但在2007—2013年間相關(guān)性呈減弱趨勢；農(nóng)機專業(yè)合作社2007與2013年的Moran′I 指數(shù)分別為-0.165、-0.193,且分別在1%、5%的顯著性水平上通過了檢驗,表明京津冀農(nóng)機專業(yè)合作社呈現(xiàn)明顯的空間負相關(guān)性,且在2007—2013年間相關(guān)性呈增強趨勢。總之,京津冀農(nóng)民專業(yè)合作社分布存在明顯的空間負相關(guān)性。

表1 京津冀農(nóng)民專業(yè)合作社空間自相關(guān)性分析

*、**、***分別表示變量在10%、5%、1%水平下通過顯著性檢驗

3 京津冀農(nóng)民專業(yè)合作社分布的影響因素解析

前文分析顯示京津冀合作社數(shù)量及分布的影響因素包含地區(qū)經(jīng)濟發(fā)展水平、市場需求、人力資本。另外第一產(chǎn)業(yè)發(fā)展水平是合作社發(fā)展的關(guān)鍵,政府對農(nóng)業(yè)的支持也對合作社發(fā)展起到推動作用。再者,京津冀是我國三大城市群之一,作為高度同城化、一體化的城市群體,其城市間的發(fā)展是緊密聯(lián)系的[59],農(nóng)民專業(yè)合作社的發(fā)展也不例外,此外,前文分析表明京津冀農(nóng)民專業(yè)合作社分布存在明顯的空間負相關(guān)性,因此需要將空間因素納入影響因素分析?；诖藰?gòu)建兩種空間計量模型(空間誤差、空間滯后),以京津冀各市各類專業(yè)合作社數(shù)量為被解釋變量(y),種植業(yè)專業(yè)合作社數(shù)量為y1、養(yǎng)殖業(yè)專業(yè)合作社數(shù)量為y2、供銷專業(yè)合作社數(shù)量為y3、農(nóng)機專業(yè)合作社數(shù)量為y4；以人均GDP(x1)、職工平均工資(x2)、第一產(chǎn)業(yè)從業(yè)人員數(shù)(x3)、第一產(chǎn)業(yè)產(chǎn)值(x4)、農(nóng)林水務(wù)財政支出(x5)為解釋變量。人均GDP代表地區(qū)經(jīng)濟發(fā)展水平,預(yù)期經(jīng)濟發(fā)展水平越高的地方合作社數(shù)量越多；職工平均工資反映了收入水平,收入水平越高,消費水平相對也會越高,消費水平的上升會擴大市場需求,預(yù)期市場需求對合作社發(fā)展的影響為正；第一產(chǎn)業(yè)從業(yè)人員數(shù)代表人力資本,管理學(xué)認為人是組織發(fā)展的核心,預(yù)期人力資本對合作社發(fā)展的影響為正；第一產(chǎn)業(yè)產(chǎn)值代表第一產(chǎn)業(yè)發(fā)展水平,預(yù)期其對合作社發(fā)展的影響為正；為更好的符合三農(nóng)政策和支持以水利為重點的農(nóng)業(yè)基礎(chǔ)設(shè)施,我國對于農(nóng)業(yè)方面的支出主要按照農(nóng)林水務(wù)財政支出計算,因此農(nóng)林水務(wù)財政支出可以很好的代表政府對農(nóng)業(yè)的支持,預(yù)期其對合作社發(fā)展的影響。對時間序列數(shù)據(jù)取自然對數(shù)能夠使數(shù)據(jù)趨于線性化并在一定程度上消除異方差,而且不會改變原始數(shù)據(jù)屬性,故模型采用取自然對數(shù)(In)后的指標值,其中空間權(quán)重矩陣采用queen鄰接法得到。

先運用LM(Lagrange Multiplier)及Robust LM檢驗對空間滯后模型與空間誤差模型進行選擇。四類農(nóng)民專業(yè)合作社的LM與RLM檢驗結(jié)果見下表。結(jié)果表明種植業(yè)與供銷專業(yè)合作社的空間誤差模型的LM與RLM均大于空間滯后模型,且均在1%水平下通過了顯著性檢驗,因此預(yù)期選用空間誤差模型；而養(yǎng)殖業(yè)與農(nóng)機專業(yè)合作社的空間滯后模型的LM與RLM均大于空間誤差模型,且均在1%與5%水平下通過了顯著性檢驗,因此預(yù)期選用空間滯后模型。綜上,本文預(yù)期種植業(yè)與供銷專業(yè)合作社選用空間誤差模型,養(yǎng)殖業(yè)與農(nóng)機專業(yè)合作社選用空間滯后模型。

表2 空間滯后與空間誤差模型的LM與RLM檢驗

進一步對各類專業(yè)合作社的空間滯后模型與空間誤差模型進行擬合分析,探討經(jīng)濟地理因素對農(nóng)民專業(yè)合作社發(fā)展的影響(表3)。種植業(yè)、供銷專業(yè)合作社的LogL是空間誤差模型大于空間滯后模型,AIC與SC是空間誤差模型小于空間滯后模型,表明擬合效果為空間誤差優(yōu)于空間滯后,因此確定選用空間誤差模型；養(yǎng)殖業(yè)、農(nóng)機專業(yè)合作社的LogL是空間滯后模型大于空間誤差模型,AIC與SC是空間滯后模型小于空間誤差模型,表明擬合效果為空間滯后優(yōu)于空間誤差,因此確定選用空間滯后模型。綜上,種植業(yè)與供銷專業(yè)合作社選用空間誤差模型,養(yǎng)殖業(yè)與農(nóng)機專業(yè)合作社選用空間滯后模型。

表3 空間滯后與空間誤差模型的擬合分析

因此,對種植業(yè)與供銷專業(yè)合作社著重于空間誤差模型的分析,對養(yǎng)殖業(yè)與農(nóng)機專業(yè)合作社則著重于空間滯后模型的分析(表4、5),結(jié)果顯示：

表4 種植業(yè)與供銷專業(yè)合作社空間誤差模型估計結(jié)果

X1:人均國內(nèi)生產(chǎn)總值 Per capita Gross Domestic Product；X2:職工平均工資 Average wage of employees；X3:第一產(chǎn)業(yè)從業(yè)人員數(shù) Number of employees in the primary industry；X4:第一產(chǎn)業(yè)產(chǎn)值 Output value of primary industry；X5: 農(nóng)林水務(wù)財政支出 Financial expenditure of agriculture, forestry and water；λ:地理因素 Geographical factor

表5 養(yǎng)殖業(yè)與農(nóng)機專業(yè)合作社空間滯后模型估計結(jié)果

WlnY：地理因素 Geographical factor

(1)地理因素對京津冀種植業(yè)、養(yǎng)殖業(yè)、供銷與農(nóng)機專業(yè)合作社的發(fā)展具有顯著的正向影響。地理因素對京津冀種植業(yè)、養(yǎng)殖業(yè)、供銷與農(nóng)機專業(yè)合作社的影響系數(shù)分別為0.4969、0.5799、0.5950、0.5280,且均在1%的水平上通過了顯著性檢驗,說明鄰近區(qū)域相似的經(jīng)濟地理環(huán)境可以促進京津冀農(nóng)民專業(yè)合作社的交流與集聚發(fā)展,再一次驗證了考慮空間影響因素的必要性。根據(jù)回歸系數(shù)大小可知影響強度排序為供銷>養(yǎng)殖業(yè)>農(nóng)機>種植業(yè)。

(2)經(jīng)濟發(fā)展水平對京津冀種植業(yè)、養(yǎng)殖業(yè)、供銷與農(nóng)機專業(yè)合作社的發(fā)展具有顯著的正向影響。經(jīng)濟發(fā)展水平對京津冀種植業(yè)、養(yǎng)殖業(yè)、供銷與農(nóng)機專業(yè)合作社的影響系數(shù)分別為0.0115、0.1321、0.1073、0.2038,且均在10%的水平上通過了顯著性檢驗,說明經(jīng)濟發(fā)展水平對合作社的發(fā)展具有拉動效應(yīng),與預(yù)期相符,但拉動效應(yīng)較弱。根據(jù)回歸系數(shù)大小可知影響強度排序為農(nóng)機>養(yǎng)殖業(yè)>供銷>種植業(yè)。

(3)市場需求對京津冀種植業(yè)、養(yǎng)殖業(yè)、供銷專業(yè)合作社的發(fā)展具有顯著正向影響,但對農(nóng)機專業(yè)合作社的發(fā)展無顯著影響。市場需求對京津冀種植業(yè)、養(yǎng)殖業(yè)、供銷專業(yè)合作社的影響系數(shù)分別為3.5474、0.8807、1.6750,且均在1%的水平上通過了顯著性檢驗,但對農(nóng)機專業(yè)合作社的影響未通過檢驗,那是因為本文的市場需求變量實質(zhì)是人民的生活水平,很大程度上代表生活上的需求,對和食物等常用生活品密切相關(guān)的種植業(yè)、養(yǎng)殖業(yè)、供銷合作社的影響自然大于農(nóng)機等科技產(chǎn)品,再者,農(nóng)業(yè)機械購買成本較高,且大多為農(nóng)機公司承包用以實現(xiàn)農(nóng)業(yè)的機械化生產(chǎn),故而影響不顯著。根據(jù)回歸系數(shù)大小可知影響強度排序為種植業(yè)>供銷>養(yǎng)殖業(yè)。

(4)人力資本對京津冀種植業(yè)、養(yǎng)殖業(yè)專業(yè)合作社發(fā)展具有顯著的正向影響,對供銷與農(nóng)機專業(yè)合作社發(fā)展的影響不顯著。人力資本對京津冀種植業(yè)、養(yǎng)殖業(yè)專業(yè)合作社的影響系數(shù)分別為1.5314、0.0629,且均在5%的水平上通過了顯著性檢驗,而對供銷、農(nóng)機專業(yè)合作社的影響未通過檢驗,那是因為供銷與農(nóng)機專業(yè)合作社屬于農(nóng)業(yè)服務(wù)業(yè),對人力資本的綜合素質(zhì)要求較高,第一產(chǎn)業(yè)從業(yè)人員中能夠從事此類管理與技術(shù)工作的人數(shù)占比較少,因而對供銷與農(nóng)機專業(yè)合作社的影響不顯著。根據(jù)回歸系數(shù)大小可知影響強度排序為種植業(yè)>養(yǎng)殖業(yè)。

(5)第一產(chǎn)業(yè)發(fā)展水平對京津冀種植業(yè)、養(yǎng)殖業(yè)、供銷與農(nóng)機專業(yè)合作社的發(fā)展具有顯著的正向影響。第一產(chǎn)業(yè)發(fā)展水平對京津冀種植業(yè)、養(yǎng)殖業(yè)、供銷與農(nóng)機專業(yè)合作社的影響系數(shù)分別為1.9737、1.2977、0.5751、2.3206,且均在10%及以上水平上通過了顯著性檢驗,與預(yù)期相符。根據(jù)回歸系數(shù)大小可知影響強度排序為農(nóng)機>種植業(yè)>養(yǎng)殖業(yè)>供銷。

(6)政府支持對京津冀養(yǎng)殖業(yè)、供銷、農(nóng)機專業(yè)合作社的發(fā)展具有顯著的正向影響,對種植業(yè)專業(yè)合作社的發(fā)展影響不顯著。政府支持對養(yǎng)殖業(yè)、供銷、農(nóng)機專業(yè)合作社的影響系數(shù)分別為0.5121、0.2145、0.4425,且均在10%及以上水平上通過了顯著性檢驗,對種植業(yè)的影響未通過檢驗,可能是因為政府對農(nóng)業(yè)的財政支持資金的分配上存在差異,具體原因還需整理更明細的財政支出情況。根據(jù)回歸系數(shù)大小可知影響強度排序為養(yǎng)殖業(yè)>農(nóng)機>供銷。

(7)根據(jù)回歸系數(shù)大小可知各影響因素對種植業(yè)專業(yè)合作社的影響強度排序為市場需求>第一產(chǎn)業(yè)發(fā)展水平>人力資本>地理因素>地區(qū)經(jīng)濟發(fā)展水平；對養(yǎng)殖業(yè)專業(yè)合作社的影響強度排序為市場需求>地理因素>第一產(chǎn)業(yè)發(fā)展水平>政府支持>地區(qū)經(jīng)濟發(fā)展水平>人力資本；對供銷專業(yè)合作社影響強度排序為第一產(chǎn)業(yè)發(fā)展水平>市場需求>地理因素>政府支持>地區(qū)經(jīng)濟發(fā)展水平；對農(nóng)機專業(yè)合作社影響強度排序為第一產(chǎn)業(yè)發(fā)展水平>地理因素>政府支持>地區(qū)經(jīng)濟發(fā)展水平。且以上影響均為正向。

4 結(jié)論與討論

4.1 主要結(jié)論

(1)京津冀種植業(yè)、養(yǎng)殖業(yè)、供銷與專業(yè)合作社數(shù)量總體上呈逐年增長趨勢,且合作社數(shù)量為種植業(yè)>養(yǎng)殖業(yè)>供銷>農(nóng)機,增長速度亦是種植業(yè)>養(yǎng)殖業(yè)>供銷>農(nóng)機。

(2)京津冀種植業(yè)與養(yǎng)殖業(yè)專業(yè)合作社數(shù)量增長速度為北京>天津>石家莊>其他市,北京、天津與石家莊的種植業(yè)專業(yè)合作社有明顯的集聚分布效應(yīng),雙核特征顯著,合作社分布隨時間遷移由兩核向東南部擴散,秦皇島、唐山、滄州、衡水、邢臺、邯鄲組成東南發(fā)展帶,聯(lián)通兩核,呈現(xiàn)“兩核一帶”空間格局,合作社數(shù)量東南部多于西北部,廊坊、保定夾在“兩核一帶”中間形成中部塌陷區(qū)；供銷專業(yè)合作社數(shù)量增長速度較緩慢(北京>天津>保定>其他市),合作社分布的梯度特征明顯,合作社數(shù)量表現(xiàn)為中部>南部>北部,北京、天津呈現(xiàn)明顯的極化效應(yīng)；農(nóng)機專業(yè)合作社數(shù)量增長速度北京、天津>石家莊>其他市,北京與天津呈現(xiàn)互相追趕態(tài)勢,合作社分布的梯度特征明顯,合作社數(shù)量表現(xiàn)為中部>南部>北部,2007—2011年北京、天津呈現(xiàn)明顯的極化效應(yīng),2012、2013年演化為北京、天津與石家莊雙核效應(yīng),2013年“雙核一帶”格局初現(xiàn),廊坊、保定夾在“兩核一帶”中間形成中部塌陷區(qū),地形因素導(dǎo)致西、北部張家口、承德等城市的農(nóng)機專業(yè)合作社數(shù)量較低。京津冀農(nóng)民專業(yè)合作社分布存在明顯的空間負相關(guān)性。

(3)空間誤差模型對種植業(yè)與供銷專業(yè)合作社有更好的擬合效果,而空間滯后模型對養(yǎng)殖業(yè)與農(nóng)機專業(yè)合作社有更好的擬合效果；地理因素對京津冀四類專業(yè)合作社的發(fā)展具有顯著的正向影響,與預(yù)期相符,且影響強度排序為供銷>養(yǎng)殖業(yè)>農(nóng)機>種植業(yè)；經(jīng)濟發(fā)展水平對四類專業(yè)合作社的發(fā)展具有顯著的正向影響,與預(yù)期相符,且影響強度排序為農(nóng)機>養(yǎng)殖業(yè)>供銷>種植業(yè)；市場需求對種植業(yè)、養(yǎng)殖業(yè)、供銷專業(yè)合作社的發(fā)展具有顯著正向影響,但對農(nóng)機專業(yè)合作社的發(fā)展無顯著影響。且影響強度排序為種植業(yè)>供銷>養(yǎng)殖業(yè)；人力資本對種植業(yè)、養(yǎng)殖業(yè)專業(yè)合作社發(fā)展具有顯著的正向影響,對供銷與農(nóng)機專業(yè)合作社發(fā)展的影響不顯著,且影響強度排序為種植業(yè)>養(yǎng)殖業(yè)；第一產(chǎn)業(yè)發(fā)展水平四類專業(yè)合作社的發(fā)展具有顯著的正向影響,與預(yù)期相符,且影響強度排序為農(nóng)機>種植業(yè)>養(yǎng)殖業(yè)>供銷；政府支持對養(yǎng)殖業(yè)、供銷、農(nóng)機專業(yè)合作社的發(fā)展具有顯著的正向影響,對種植業(yè)專業(yè)合作社的發(fā)展影響不顯著,且影響強度排序為養(yǎng)殖業(yè)>農(nóng)機>供銷。

(4)各影響因素對種植業(yè)專業(yè)合作社的影響強度排序為市場需求>第一產(chǎn)業(yè)發(fā)展水平>人力資本>地理因素>地區(qū)經(jīng)濟發(fā)展水平；對養(yǎng)殖業(yè)專業(yè)合作社的影響強度排序為市場需求>地理因素>第一產(chǎn)業(yè)發(fā)展水平>政府支持>地區(qū)經(jīng)濟發(fā)展水平>人力資本；對供銷專業(yè)合作社影響強度排序為第一產(chǎn)業(yè)發(fā)展水平>市場需求>地理因素>政府支持>地區(qū)經(jīng)濟發(fā)展水平；對農(nóng)機專業(yè)合作社影響強度排序為第一產(chǎn)業(yè)發(fā)展水平>地理因素>政府支持>地區(qū)經(jīng)濟發(fā)展水平。且以上影響均為正向。

4.2 討論

(1)數(shù)據(jù)尚存缺陷。依據(jù)合作社名稱里是否包含“種植”、“養(yǎng)殖”、“供銷”與“農(nóng)機”來對京津冀農(nóng)民專業(yè)合作社進行分類,首先會遺漏名稱里沒有相應(yīng)文字但本質(zhì)屬于相應(yīng)種類的合作社,還會遺漏其他類型的合作社；其次分類體系也存在差異,種植與養(yǎng)殖屬于生產(chǎn)型。供銷與農(nóng)機屬于服務(wù)型。但由于數(shù)據(jù)量的龐大,也只能運用馬克思原理抓住主要矛盾,需要說明的是本研究的農(nóng)民專業(yè)合作社并未涵蓋京津冀全部合作社,只是抽取幾類代表性的農(nóng)民專業(yè)合作社進行典型性分析。

(2)研究深度有待挖掘。影響因素分析部分主要是對京津冀四類農(nóng)民專業(yè)合作社整體的把握,缺乏地方思考,且在考慮各類型合作社影響因素時采用一刀切的方式欠妥。今后研究應(yīng)繼續(xù)深入挖掘影響京津冀城市群更小尺度地域合作社發(fā)展的主控因素,如可以對各變量對各地影響系數(shù)作GIS可視化表達,發(fā)現(xiàn)地區(qū)差異；供銷與農(nóng)機專業(yè)合作社是介于農(nóng)業(yè)與服務(wù)之間的經(jīng)濟組織,對其影響因素的分析需要綜合考慮,可適當增加服務(wù)業(yè)影響因子的比重。

(3)京津冀農(nóng)民專業(yè)合作社的發(fā)展可以從地區(qū)經(jīng)濟發(fā)展水平、市場需求、人力資本、第一產(chǎn)業(yè)發(fā)展水平、政府支持以及集聚規(guī)模發(fā)展等入手進行宏觀調(diào)控,但需“因社而異,各有側(cè)重”,如對于種植業(yè)與養(yǎng)殖業(yè)專業(yè)合作社而言,應(yīng)著力擴大市場需求,對于供銷與農(nóng)機專業(yè)合作社而言,應(yīng)大力提升第一產(chǎn)業(yè)水平等。