薛 平，韓昕儒，張麗麗，王秀東

（1.中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)經(jīng)濟(jì)與發(fā)展研究所，北京 100081；2.山東省煙臺(tái)市農(nóng)業(yè)科學(xué)院，山東 煙臺(tái) 264010）

一、引 言

改革開放以來(lái)，我國(guó)畜牧業(yè)發(fā)展成果顯著，逐漸成為農(nóng)業(yè)農(nóng)村經(jīng)濟(jì)重要支柱產(chǎn)業(yè)（劉剛等，2018；陳偉生等，2019）。肉羊產(chǎn)業(yè)作為畜牧業(yè)重要產(chǎn)業(yè)之一，存欄量、出欄量和羊肉產(chǎn)量均增長(zhǎng)較快。其中，羊肉產(chǎn)量增幅最為明顯，由1990年的106.8萬(wàn)噸增至2019年的487.5萬(wàn)噸，增長(zhǎng)3.6倍①數(shù)據(jù)來(lái)源：國(guó)家統(tǒng)計(jì)局。。隨著我國(guó)居民收入水平不斷提高，食物消費(fèi)結(jié)構(gòu)也隨之改變（Gao等，2020），肉類食物消費(fèi)由以豬肉消費(fèi)為主轉(zhuǎn)向豬牛羊肉及禽肉多元化消費(fèi)。而且《肉牛肉羊優(yōu)勢(shì)區(qū)域發(fā)展規(guī)劃（2003～2007年）》《全國(guó)肉羊優(yōu)勢(shì)區(qū)域布局規(guī)劃（2008～2015年）》等政策均提出充分發(fā)揮各地區(qū)肉羊生產(chǎn)優(yōu)勢(shì)，科學(xué)合理的建設(shè)肉羊生產(chǎn)優(yōu)勢(shì)區(qū)域已成為我國(guó)肉羊產(chǎn)業(yè)發(fā)展的重要方向。因此，研究我國(guó)肉羊生產(chǎn)的空間集聚特征及其影響因素對(duì)肉羊生產(chǎn)向優(yōu)勢(shì)產(chǎn)區(qū)集聚，充分發(fā)揮優(yōu)勢(shì)產(chǎn)區(qū)資源優(yōu)勢(shì)、滿足城鄉(xiāng)居民肉類多元化消費(fèi)需求、促進(jìn)肉羊產(chǎn)業(yè)可持續(xù)健康發(fā)展等具有重要意義。

產(chǎn)業(yè)空間集聚對(duì)資源利用率（季書涵等，2016）、生產(chǎn)效率（Tao等，2019）、產(chǎn)業(yè)專業(yè)化（李二玲等，2019）、產(chǎn)業(yè)升級(jí)（袁鵬等，2021）、經(jīng)濟(jì)發(fā)展（Wei等，2019；黃慶華等，2020）等方面具有重要影響。眾多國(guó)內(nèi)外學(xué)者對(duì)農(nóng)業(yè)領(lǐng)域集聚問(wèn)題展開研究，包括農(nóng)產(chǎn)品和畜禽產(chǎn)品兩個(gè)領(lǐng)域，指標(biāo)主要選取產(chǎn)量、出欄量等，方法主要采用生產(chǎn)集中度指數(shù)、區(qū)位熵指數(shù)、區(qū)域重心分析法、全局Moran's I指數(shù)等。在農(nóng)產(chǎn)品領(lǐng)域，國(guó)內(nèi)外學(xué)者對(duì)糧食（趙穎文等，2020）、蔬菜（吳建寨等，2016）、水果（張聰穎等，2018）、棉花（肖海峰等，2018）等產(chǎn)業(yè)的集聚問(wèn)題進(jìn)行研究。如Erling等（2017）選取產(chǎn)量指標(biāo)，運(yùn)用重心模型、基尼系數(shù)、全局Moran's I指數(shù)對(duì)11種農(nóng)作物產(chǎn)業(yè)集聚的演化模式探討，發(fā)現(xiàn)中國(guó)種植業(yè)在空間分布上呈集聚趨勢(shì)，區(qū)域?qū)I(yè)化和多樣化程度逐步加強(qiáng)。在畜禽產(chǎn)品領(lǐng)域，國(guó)內(nèi)外學(xué)者對(duì)肉牛（楊春等，2013）、肉雞（劉愛民等，2011）、生豬（趙俊偉等，2019）、肉羊（丁存振等，2018）等產(chǎn)業(yè)的集聚問(wèn)題進(jìn)行研究。如Csonka等（2019）采用養(yǎng)豬場(chǎng)規(guī)模為指標(biāo)，分析匈牙利生豬產(chǎn)業(yè)的集聚效應(yīng)，認(rèn)為在匈牙利豬肉行業(yè)中個(gè)體農(nóng)場(chǎng)在生豬養(yǎng)殖上具有較明顯的產(chǎn)業(yè)集聚效應(yīng)，而企業(yè)農(nóng)場(chǎng)在生豬養(yǎng)殖上的集聚效應(yīng)并不明顯；王歡等（2017）選取生豬出欄量為指標(biāo)，運(yùn)用生產(chǎn)集中度指數(shù)及區(qū)域重心分析法對(duì)生豬產(chǎn)業(yè)布局問(wèn)題進(jìn)行研究，認(rèn)為生豬產(chǎn)業(yè)存在集聚現(xiàn)象。

產(chǎn)業(yè)集聚影響因素的研究方面，國(guó)內(nèi)外學(xué)者主要建立計(jì)量經(jīng)濟(jì)學(xué)模型和空間計(jì)量經(jīng)濟(jì)學(xué)模型進(jìn)行實(shí)證研究，并歸納影響因素。如羅萬(wàn)純等（2005）運(yùn)用固定效應(yīng)模型對(duì)影響糧食區(qū)域格局的耕地資源、非農(nóng)收入、糧食單產(chǎn)、畜牧業(yè)發(fā)展?fàn)顩r等因素實(shí)證分析；葉志標(biāo)等（2019）運(yùn)用空間計(jì)量經(jīng)濟(jì)學(xué)模型對(duì)影響小麥空間集聚的氣候、技術(shù)進(jìn)步、生產(chǎn)投入和其他社會(huì)因素定量分析；Rosenthal等（2004）認(rèn)為自然優(yōu)勢(shì)、國(guó)內(nèi)市場(chǎng)、消費(fèi)機(jī)會(huì)和尋租等因素對(duì)產(chǎn)業(yè)集聚均具有正向影響，產(chǎn)業(yè)集聚在產(chǎn)業(yè)、空間和時(shí)間三個(gè)維度上存在相互作用。

學(xué)術(shù)界對(duì)產(chǎn)業(yè)集聚問(wèn)題的研究相對(duì)豐富，為本文提供借鑒，但仍存在一些角度值得深入探索。首先，產(chǎn)業(yè)集聚研究在農(nóng)產(chǎn)品領(lǐng)域相對(duì)豐富，涵蓋糧食、蔬菜、水果等產(chǎn)業(yè)，而在畜禽產(chǎn)品領(lǐng)域，主要集中在生豬產(chǎn)業(yè)，而對(duì)肉羊、肉牛等草食性動(dòng)物研究相對(duì)較少；其次，已有研究選取的產(chǎn)業(yè)集聚特征分析指標(biāo)多樣且方法豐富，但較少關(guān)注產(chǎn)業(yè)在空間上的集聚特征；最后，學(xué)術(shù)界對(duì)產(chǎn)業(yè)集聚影響因素的研究較為全面，主要包括自然、經(jīng)濟(jì)、市場(chǎng)、社會(huì)、政策等因素，但較少關(guān)注空間因素。鑒于此，本文將以肉羊?yàn)橹饕芯繉?duì)象，選取羊肉產(chǎn)量為主要指標(biāo)，運(yùn)用全局Moran's I指數(shù)和局部Moran's I指數(shù)，并建立空間計(jì)量經(jīng)濟(jì)學(xué)模型，將空間因素納入肉羊生產(chǎn)集聚研究之中，以期完善肉羊產(chǎn)業(yè)和產(chǎn)業(yè)集聚兩方面研究。

二、研究方法與數(shù)據(jù)說(shuō)明

（一）研究方法

1.全局Moran's I指數(shù)

全局Moran's I指數(shù)一方面可檢驗(yàn)地區(qū)間是否存在空間相關(guān)性，另一方面可為建立空間計(jì)量經(jīng)濟(jì)學(xué)模型提供依據(jù)。目前，常用的檢驗(yàn)方法包括Moran's I指數(shù)檢驗(yàn)、Rao's Score檢驗(yàn)和拉格朗日乘數(shù)（LM）檢驗(yàn)（Anselin，1988）。本文選用學(xué)術(shù)界使用較為廣泛的Moran's I指數(shù)（王釗等，2013），選取各省羊肉產(chǎn)量占全國(guó)羊肉產(chǎn)量的比重為指標(biāo)，進(jìn)行空間相關(guān)性檢驗(yàn)。公式如下。

其中，I為全局Moran's I值，E(I)為期望值，VAR(I)為方差。

2.局部Moran's I指數(shù)

考慮到全局Moran's I指數(shù)判斷的是整體的空間相關(guān)性，但未對(duì)某省與鄰省的具體空間相關(guān)性特征進(jìn)行判斷。因此，引入局部Moran′s I指數(shù)進(jìn)一步分析鄰省間肉羊生產(chǎn)空間相關(guān)性。計(jì)算公式如下。

式中，zi和zj是觀測(cè)值i和j的標(biāo)準(zhǔn)化值，其余參數(shù)與全局Moran′s I指數(shù)公式定義保持一致。從上式可知：局部Moran′s I指數(shù)由兩部分組成，一是zi，二是空間滯后向量高值為H，低值為L(zhǎng)，可識(shí)別出高高集聚區(qū)（H—H）、高低集聚區(qū)（H—L）、低高集聚區(qū)（L—H）和低低集聚區(qū)（L—L）四種空間關(guān)聯(lián)模式。

3.空間面板模型

若地區(qū)間存在空間相關(guān)性，將構(gòu)建空間面板模型進(jìn)行估計(jì)，并采用極大似然估計(jì)（MLE）得到可信的參數(shù)估計(jì)值。本文將建立空間滯后面板模型或空間誤差面板模型對(duì)肉羊生產(chǎn)的集聚效應(yīng)進(jìn)行研究（Anselin，2008）。兩個(gè)模型分別表示如下。

在以上兩個(gè)模型中，yit為t年i地區(qū)羊肉產(chǎn)量占全國(guó)羊肉產(chǎn)量比重，βi為系數(shù)，xit為解釋變量，ρ為空間自回歸系數(shù)，λ為空間自相關(guān)系數(shù)，εit和?it為隨機(jī)誤差項(xiàng)；β0為截距項(xiàng)；n為地區(qū)個(gè)數(shù)；Wij為空間鄰接矩陣。當(dāng)ρ=0時(shí)式（5）為普通面板模型。

（二）數(shù)據(jù)說(shuō)明

本文研究區(qū)域包括北京、天津、河北、山西、內(nèi)蒙古、遼寧、吉林、黑龍江、上海、江蘇、浙江、安徽、福建、江西、山東、河南、湖北、湖南、廣東、廣西、四川、貴州、云南、西藏、陜西、甘肅、青海、寧夏、新疆，共29個(gè)?。ㄖ陛犑?、自治區(qū)）。為保持統(tǒng)計(jì)口徑相同，將重慶劃入四川，海南劃入廣東。采用數(shù)據(jù)為1990～2017年公開統(tǒng)計(jì)數(shù)據(jù)，主要來(lái)源于《中國(guó)統(tǒng)計(jì)年鑒》《中國(guó)農(nóng)村統(tǒng)計(jì)年鑒》《全國(guó)農(nóng)產(chǎn)品成本收益資料匯編》、國(guó)家統(tǒng)計(jì)局等。

三、肉羊生產(chǎn)空間集聚特征分析

（一）肉羊生產(chǎn)的全局相關(guān)特征分析

1990～2017年各省全局Moran's I指數(shù)如表1所示。1990～2017年，全局Moran's I指數(shù)在0.2～0.5之間變動(dòng)，除2005年顯著性水平略大于5%外，其余年份均通過(guò)5%顯著性水平檢驗(yàn)。表明中國(guó)肉羊生產(chǎn)整體上呈顯著空間正相關(guān)趨勢(shì)，且空間相關(guān)程度相對(duì)較高。從全局Moran's I指數(shù)的變化趨勢(shì)來(lái)看，1990～1996年，全局Moran's I指數(shù)呈下降趨勢(shì)，由0.460降至0.227，降幅為50.7%；1997～2004年，全局Moran's I指數(shù)波動(dòng)較為平穩(wěn)；2005～2017年，全局Moran's I指數(shù)呈增長(zhǎng)趨勢(shì)，由0.197增至0.271，增幅為37.6%。表明我國(guó)肉羊生產(chǎn)空間相關(guān)程度自2005年以來(lái)不斷增加。

表1 1990～2017年中國(guó)肉羊生產(chǎn)全局Moran's I指數(shù)

（二）肉羊生產(chǎn)的局部空間相關(guān)性特征

因肉羊生產(chǎn)局部空間相關(guān)性特征變化趨勢(shì)大體相同，故僅將1990、2000、2010、2017四個(gè)年份局部空間關(guān)聯(lián)結(jié)果列于表2中。1990年有12個(gè)?。ㄖ陛犑?、自治區(qū)）位于高高集聚區(qū)，分別為河北、山西、內(nèi)蒙古、安徽、山東、河南、四川、西藏、陜西、甘肅、青海、新疆；有9個(gè)位于低低集聚區(qū)，分別為云南、湖北、貴州、浙江、廣西、江西、湖南、福建、廣東；有7個(gè)位于高低集聚區(qū)，只有1個(gè)位于低高集聚區(qū)。到2017年有12個(gè)省（直轄市、自治區(qū)）位于高高集聚區(qū)，分別為山西、內(nèi)蒙古、寧夏、黑龍江、山東、河南、四川、陜西、甘肅、青海、新疆、湖北；有9個(gè)位于低低集聚區(qū)，分別為北京、天津、貴州、浙江、廣西、江西、福建、廣東、江蘇；有4個(gè)位于高低集聚區(qū)；有4個(gè)位于低高集聚區(qū)。綜上可知，1990年和2017年的高高集聚區(qū)主要為肉羊生產(chǎn)集中度較高省，低低集聚區(qū)主要為肉羊生產(chǎn)集中度較低的省；1990和2017年分布在第一象限和第三象限的省占全部省比重均超過(guò)70%，表明我國(guó)肉羊生產(chǎn)具有較為顯著的空間集聚特征。

表2 部分年份的肉羊生產(chǎn)局部集聚結(jié)果

四、肉羊生產(chǎn)空間集聚影響因素實(shí)證分析

（一）變量設(shè)置

運(yùn)用各省羊肉產(chǎn)量占全國(guó)羊肉產(chǎn)量比重可較好反映肉羊生產(chǎn)集聚的基本特征，因此，選取各省羊肉產(chǎn)量占全國(guó)羊肉產(chǎn)量比重為被解釋變量，選取肉羊生產(chǎn)集聚的影響因素為解釋變量。具體如下。

1.自然資源條件

自然資源條件是畜牧業(yè)發(fā)展的前提和基礎(chǔ)，一方面肉羊生產(chǎn)傾向于向飼草或玉米等糧食作物資源豐富的地區(qū)集中（孫致陸等，2013），且當(dāng)?shù)剞r(nóng)牧民也會(huì)受自然資源條件影響形成肉羊養(yǎng)殖習(xí)慣，故選取糧食總產(chǎn)量衡量各區(qū)域自然資源條件豐裕程度；另一方面旱災(zāi)、水災(zāi)、冰凍等自然災(zāi)害通過(guò)影響飼草生長(zhǎng)、糧食生產(chǎn)及引發(fā)肉羊疫病等途徑阻礙肉羊產(chǎn)業(yè)發(fā)展，考慮到旱災(zāi)對(duì)各地區(qū)飼草生長(zhǎng)、糧食產(chǎn)量等肉羊產(chǎn)業(yè)發(fā)展所需生產(chǎn)資料的影響較大，故選取旱災(zāi)成災(zāi)面積衡量各區(qū)域自然災(zāi)害情況。

2.經(jīng)濟(jì)發(fā)展水平

隨著人們生活水平不斷提高，膳食營(yíng)養(yǎng)標(biāo)準(zhǔn)隨之提高，羊肉作為高品質(zhì)肉類食物，深受人們喜愛。根據(jù)經(jīng)濟(jì)學(xué)中供給需求定理可知，居民對(duì)羊肉需求增加，會(huì)促使羊肉供給增加，進(jìn)而促進(jìn)羊肉價(jià)格上漲，對(duì)農(nóng)戶養(yǎng)殖肉羊產(chǎn)生良好的激勵(lì)作用；從另一個(gè)角度來(lái)看，肉羊產(chǎn)業(yè)發(fā)展會(huì)占用大量土地并產(chǎn)生大量排泄物污染環(huán)境，而隨著經(jīng)濟(jì)的發(fā)展，商業(yè)用地需求不斷增加和對(duì)環(huán)境關(guān)注度逐漸升高阻礙肉羊產(chǎn)業(yè)發(fā)展?；诖?，選取城鎮(zhèn)化水平衡量各區(qū)域經(jīng)濟(jì)發(fā)展水平。

3.比較效益

肉羊生產(chǎn)的比較效益主要從生產(chǎn)者的選擇行為及從事不同行業(yè)的機(jī)會(huì)成本兩個(gè)方面衡量。一方面，養(yǎng)殖肉羊的決策會(huì)受到養(yǎng)殖肉羊與養(yǎng)殖其他畜禽的比較效益影響，只有當(dāng)養(yǎng)殖肉羊的效益相對(duì)高于養(yǎng)殖其他畜禽時(shí)，農(nóng)民才會(huì)選擇養(yǎng)殖肉羊；另一方面，隨著農(nóng)業(yè)生產(chǎn)效率和非農(nóng)收入提高，越來(lái)越多的農(nóng)民選擇進(jìn)城務(wù)工以提高家庭收入，農(nóng)民進(jìn)城務(wù)工與養(yǎng)殖肉羊的收入差距在一定程度上影響農(nóng)民肉羊養(yǎng)殖行為?；诖?，選取肉羊產(chǎn)業(yè)內(nèi)比較優(yōu)勢(shì)，即羊肉產(chǎn)值占豬牛羊肉總產(chǎn)值比重，非農(nóng)收入比重兩個(gè)指標(biāo)衡量肉羊生產(chǎn)的產(chǎn)業(yè)內(nèi)和產(chǎn)業(yè)間比較優(yōu)勢(shì)。

4.農(nóng)業(yè)勞動(dòng)力

在我國(guó)，肉羊生產(chǎn)方式相對(duì)較為傳統(tǒng)，主要靠投入大量勞動(dòng)力養(yǎng)殖肉羊，屬于勞動(dòng)密集型產(chǎn)業(yè)。在肉羊養(yǎng)殖上，勞動(dòng)力投入數(shù)量增多，肉羊養(yǎng)殖精細(xì)程度及肉羊飼草料配比水平相對(duì)提高，有利于肉羊胴體重增加和羔羊存活率提高?；诖耍x取農(nóng)林牧漁從業(yè)人員數(shù)量衡量農(nóng)業(yè)勞動(dòng)力數(shù)量對(duì)肉羊生產(chǎn)集聚的影響。

Lifestyle and diet are highly influential factors for kidney stone formation in the United States, especially lifestyle,in comparison to the world (Figure 4).

5.技術(shù)進(jìn)步

在肉羊生產(chǎn)中，技術(shù)進(jìn)步包括科學(xué)飼養(yǎng)方式的提出、疫病防控流程的完善、優(yōu)良肉羊品種的選育、有效生產(chǎn)組織模式的提出等。隨著農(nóng)業(yè)技術(shù)的進(jìn)步、專業(yè)技術(shù)人才投入的增加，肉羊生產(chǎn)技術(shù)出現(xiàn)大幅提升，對(duì)肉羊產(chǎn)業(yè)發(fā)展產(chǎn)生明顯促進(jìn)作用?；诖?，選取羊肉單產(chǎn)為衡量技術(shù)進(jìn)步對(duì)肉羊生產(chǎn)的貢獻(xiàn)。

6.政策因素

政府產(chǎn)業(yè)政策對(duì)產(chǎn)業(yè)發(fā)展具有一定積極作用，如畜牧業(yè)良種補(bǔ)貼、肉羊生產(chǎn)補(bǔ)貼。在2003與2008年分別發(fā)布的政策規(guī)劃對(duì)肉羊生產(chǎn)的優(yōu)勢(shì)區(qū)域明確劃分，并提出了促進(jìn)肉羊生產(chǎn)向優(yōu)勢(shì)產(chǎn)區(qū)集聚的規(guī)劃?；诖?，選取政策虛擬變量分析2003年與2008年政策的影響。

7.空間互動(dòng)效應(yīng)

肉羊生產(chǎn)集聚會(huì)受到相鄰地區(qū)相關(guān)因素影響，在經(jīng)濟(jì)地理上存在一定空間關(guān)聯(lián)性。具體而言，肉羊生產(chǎn)集聚具有一定區(qū)域特征，相鄰地區(qū)因自然條件、資源稟賦、肉羊養(yǎng)殖習(xí)慣等具有一定相似性，從而可能具有較為相似的肉羊生產(chǎn)集聚特征，也可能某一地區(qū)肉羊生產(chǎn)規(guī)模較大，比較優(yōu)勢(shì)較強(qiáng)，抑制相鄰地區(qū)肉羊生產(chǎn)，從而使肉羊生產(chǎn)集聚特征相反。同時(shí)，相鄰地區(qū)肉羊生產(chǎn)的互動(dòng)作用也可能影響肉羊生產(chǎn)集聚，如區(qū)域間養(yǎng)殖技術(shù)擴(kuò)散可能促進(jìn)肉羊生產(chǎn)共同發(fā)展，也可能對(duì)其產(chǎn)生負(fù)向作用；因此，將空間滯后項(xiàng)引入模型之中，分析肉羊生產(chǎn)的空間互動(dòng)效應(yīng)對(duì)肉羊生產(chǎn)集聚的影響。

（二）變量描述性統(tǒng)計(jì)分析

各變量的描述性統(tǒng)計(jì)分析結(jié)果如表3所示。各省羊肉產(chǎn)量占全國(guó)羊肉產(chǎn)量比重均值為3.4%，標(biāo)準(zhǔn)差為4.2%；糧食總產(chǎn)量均值為1 767.8萬(wàn)噸，標(biāo)準(zhǔn)差為1 400.8萬(wàn)噸；旱災(zāi)成災(zāi)面積均值為381.8千公頃，標(biāo)準(zhǔn)差為517.1千公頃，表明各年旱災(zāi)成災(zāi)面積變化較大；城鎮(zhèn)化水平均值為43.2%，標(biāo)準(zhǔn)差為18.4%，表明近年來(lái)我國(guó)城鎮(zhèn)化水平變化較大；肉羊產(chǎn)業(yè)內(nèi)比較優(yōu)勢(shì)均值為1.8%，標(biāo)準(zhǔn)差為2.6%；非農(nóng)收入比重均值為40.8%，標(biāo)準(zhǔn)差為19.8%；農(nóng)林牧漁從業(yè)人員數(shù)量均值為1 002.4萬(wàn)人，標(biāo)準(zhǔn)差為831.7萬(wàn)人；羊肉單產(chǎn)均值為15.2千克，標(biāo)準(zhǔn)差為3.3千克；2003年和2008年政策虛擬變量不再進(jìn)行描述。從各變量描述性統(tǒng)計(jì)分析結(jié)果來(lái)看，各變量取值較為合理。

表3 描述性統(tǒng)計(jì)分析

（三）實(shí)證分析

通過(guò)肉羊生產(chǎn)空間相關(guān)性分析，得出空間因素不可被忽略，否則將產(chǎn)生估計(jì)偏差。從極大似然估計(jì)檢驗(yàn)結(jié)果可知（見表4），空間滯后面板模型與空間誤差面板模型拉格朗日乘數(shù)與穩(wěn)健的拉格朗日乘數(shù)的P值顯著性相同，均在1%水平下顯著，根據(jù)已有研究，在空間滯后面板模型與空間誤差面板模型的拉格朗日乘數(shù)與穩(wěn)健的拉格朗日乘數(shù)顯著性相同的情況下，選擇拉格朗日乘數(shù)值較大的模型，故選擇空間滯后面板模型（楊春等，2010）。

表4 兩種模型的拉格朗日乘數(shù)檢驗(yàn)結(jié)果

當(dāng)確定選擇空間滯后面板模型實(shí)證分析時(shí)，需選擇固定效應(yīng)或隨機(jī)效應(yīng)進(jìn)行估計(jì)。豪斯曼檢驗(yàn)可以對(duì)選擇固定效應(yīng)還是隨機(jī)效應(yīng)進(jìn)行判斷。具體檢驗(yàn)結(jié)果如表5所示：豪斯曼檢驗(yàn)的結(jié)果為正且在1%水平下顯著，故選擇固定效應(yīng)模型進(jìn)行估計(jì)。

空間滯后面板模型的回歸結(jié)果顯示空間滯后項(xiàng)、羊肉單產(chǎn)、農(nóng)林牧漁從業(yè)人員數(shù)量、肉羊產(chǎn)業(yè)內(nèi)比較優(yōu)勢(shì)、糧食總產(chǎn)量、2003與2008年政策虛擬變量均顯著正向影響肉羊生產(chǎn)空間集聚，其中肉羊產(chǎn)業(yè)內(nèi)比較優(yōu)勢(shì)、農(nóng)林牧漁從業(yè)人員數(shù)量、糧食總產(chǎn)量影響相對(duì)較大；旱災(zāi)成災(zāi)面積雖然在5%水平上顯著，但系數(shù)較小，對(duì)肉羊生產(chǎn)集聚基本無(wú)影響；城鎮(zhèn)化水平影響不顯著。具體結(jié)果如下（見表5）。

表5 空間滯后面板模型估計(jì)結(jié)果

（1）我國(guó)肉羊生產(chǎn)集聚受到空間因素顯著正向影響，驗(yàn)證了采用空間計(jì)量經(jīng)濟(jì)學(xué)模型來(lái)分析肉羊生產(chǎn)集聚問(wèn)題是可靠的，同時(shí)表明相鄰地區(qū)的自然條件、資源稟賦、生產(chǎn)模式及肉羊生產(chǎn)技術(shù)較為相似，從而呈空間集聚態(tài)勢(shì)；

（2）肉羊產(chǎn)業(yè)內(nèi)比較優(yōu)勢(shì)在1%水平下顯著，其提高1%，各省羊肉產(chǎn)量在全國(guó)羊肉產(chǎn)量比重增加0.785%，表明在城鄉(xiāng)居民多元化消費(fèi)前提下，市場(chǎng)需求對(duì)羊肉產(chǎn)品供給產(chǎn)生一定影響，肉羊生產(chǎn)的比較優(yōu)勢(shì)影響農(nóng)戶養(yǎng)殖行為；

（3）農(nóng)林牧漁從業(yè)人員數(shù)量在1%水平下顯著，其提高1%，各省羊肉產(chǎn)量在全國(guó)羊肉產(chǎn)量比重增加0.679%。對(duì)肉羊這一勞動(dòng)密集型產(chǎn)業(yè)而言，勞動(dòng)力的豐裕程度對(duì)肉羊產(chǎn)業(yè)發(fā)展影響相對(duì)較大；

（4）糧食總產(chǎn)量在1%水平下顯著，其提高1%，各省羊肉產(chǎn)量在全國(guó)羊肉產(chǎn)量比重增加0.317%。驗(yàn)證了肉羊生產(chǎn)傾向于向糧食作物資源豐裕的地區(qū)集中的特征，同時(shí)表明在草原生態(tài)保護(hù)力度增強(qiáng)和牧草資源趨緊情況下，肉羊產(chǎn)業(yè)對(duì)糧食資源的依賴性逐漸增強(qiáng)。

（四）穩(wěn)健性檢驗(yàn)

為進(jìn)一步檢驗(yàn)空間滯后面板模型實(shí)證分析結(jié)果的穩(wěn)健性，對(duì)比普通固定效應(yīng)面板模型與空間滯后面板模型的實(shí)證分析結(jié)果。穩(wěn)健性檢驗(yàn)結(jié)果如表6所示，普通固定效應(yīng)面板模型中，肉羊產(chǎn)業(yè)內(nèi)比較優(yōu)勢(shì)與糧食總產(chǎn)量對(duì)我國(guó)肉羊生產(chǎn)集聚變動(dòng)影響程度較大，均通過(guò)了1%顯著性水平檢驗(yàn)，與空間滯后面板模型檢驗(yàn)結(jié)果基本保持一致；而農(nóng)林牧漁從業(yè)人員數(shù)量未通過(guò)顯著性檢驗(yàn)，與空間滯后面板模型實(shí)證結(jié)果不一致；羊肉單產(chǎn)及政策因素未通過(guò)顯著性檢驗(yàn)，與空間滯后面板模型實(shí)證結(jié)果不一致；非農(nóng)收入比重具有顯著正向影響，與實(shí)際情況不相符。綜上可得出空間滯后面板模型與普通固定效應(yīng)面板模型在肉羊產(chǎn)業(yè)內(nèi)比較優(yōu)勢(shì)和糧食總產(chǎn)量等重要變量具有較為一致的估計(jì)，而空間滯后面板模型修正了羊肉單產(chǎn)、農(nóng)林牧漁從業(yè)人員數(shù)量、非農(nóng)收入比重及政策因素對(duì)肉羊生產(chǎn)集聚的影響，從而使模型估計(jì)結(jié)果更符合實(shí)際情況。

表6 空間滯后面板模型的穩(wěn)健性檢驗(yàn)

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

（一）結(jié)論

本文先運(yùn)用全局Moran's I指數(shù)和局部Moran's I指數(shù)對(duì)肉羊生產(chǎn)的空間相關(guān)性特征進(jìn)行分析；再通過(guò)構(gòu)建空間滯后面板模型實(shí)證分析我國(guó)肉羊生產(chǎn)空間集聚的主要影響因素。研究發(fā)現(xiàn)，我國(guó)肉羊生產(chǎn)存在顯著的空間正相關(guān)性，鄰省間也存在顯著的空間相關(guān)性，農(nóng)林牧漁從業(yè)人員數(shù)量、肉羊產(chǎn)業(yè)內(nèi)比較優(yōu)勢(shì)和糧食總產(chǎn)量對(duì)肉羊生產(chǎn)空間集聚影響相對(duì)較大，其余因素影響相對(duì)較弱。

（二）政策建議

根據(jù)以上結(jié)論，為促進(jìn)中國(guó)肉羊生產(chǎn)向優(yōu)勢(shì)產(chǎn)區(qū)集聚，合理調(diào)整我國(guó)肉羊生產(chǎn)結(jié)構(gòu)，進(jìn)一步增強(qiáng)我國(guó)肉羊產(chǎn)業(yè)優(yōu)勢(shì)，提出以下幾點(diǎn)政策建議。

第一，注重肉羊生產(chǎn)空間因素，因地制宜制定產(chǎn)業(yè)政策。在制定肉羊產(chǎn)業(yè)政策時(shí)，不宜“大水漫灌”，而應(yīng)針對(duì)不同地區(qū)的肉羊生產(chǎn)特征有針對(duì)地制定肉羊產(chǎn)業(yè)政策，使政策對(duì)當(dāng)?shù)厝庋虍a(chǎn)業(yè)的發(fā)展更具激勵(lì)性。

第二，減少非農(nóng)就業(yè)對(duì)肉羊生產(chǎn)集聚影響，提高肉羊生產(chǎn)比較效益。在肉羊生產(chǎn)方面應(yīng)給予肉羊養(yǎng)殖戶一定生產(chǎn)補(bǔ)貼，并建立合理的肉羊養(yǎng)殖保險(xiǎn)機(jī)制如肉羊死亡險(xiǎn)、羔羊死亡險(xiǎn)等，以增加肉羊生產(chǎn)的比較效益，激勵(lì)農(nóng)戶從事肉羊生產(chǎn)。

第三，充分發(fā)揮糧食主產(chǎn)省資源優(yōu)勢(shì)，倡導(dǎo)種養(yǎng)結(jié)合發(fā)展模式。在糧食主產(chǎn)區(qū)，應(yīng)充分利用當(dāng)?shù)丶Z食生產(chǎn)優(yōu)勢(shì)，建立合理的糧改飼機(jī)制，鼓勵(lì)種養(yǎng)結(jié)合的農(nóng)業(yè)經(jīng)營(yíng)模式，促進(jìn)當(dāng)?shù)胤N植業(yè)和肉羊產(chǎn)業(yè)協(xié)同發(fā)展。