楊　丹，黃　鳳，劉自敏

(1.西南大學(xué) 經(jīng)濟管理學(xué)院，重慶 400716；2.西南大學(xué) 農(nóng)業(yè)教育發(fā)展研究中心，重慶 400716)

?

基于三階段Bootstrap EBM的異質(zhì)性合作社效率評價研究：重慶例證*

楊丹1，2，黃鳳1，劉自敏1

(1.西南大學(xué) 經(jīng)濟管理學(xué)院，重慶 400716；2.西南大學(xué) 農(nóng)業(yè)教育發(fā)展研究中心，重慶 400716)

基于對傳統(tǒng)DEA模型距離函數(shù)特征、決策單元異質(zhì)性及評價結(jié)果穩(wěn)健性等方面的缺陷，構(gòu)建三階段Bootstrap EBM方法。利用重慶市合作社調(diào)查的隨機抽樣調(diào)查數(shù)據(jù)，分析了合作社所處的外部環(huán)境變量對各類投入值的不同影響，對比分析了基于EBM距離函數(shù)的環(huán)境變量調(diào)整前后評價效率以及Bootstrap EBM調(diào)整前后的效率值及其置信區(qū)間的差異，并對比分析了不同類型合作社的效率值差異，最后基于更為準確與客觀的效率評價值，有針對性地提出準確衡量及提高合作社效率的政策建議。

異質(zhì)性；合作社；EBM距離函數(shù)；三階段DEA后一個；Bootstrap抽樣；置信區(qū)間

一、引　言

隨著我國農(nóng)村合作經(jīng)濟的發(fā)展，尤其是2007年《農(nóng)民專業(yè)合作社法》的實施，合作社在我國取得了長足的發(fā)展。2013年、2014年連續(xù)兩年的中央一號文件更是強調(diào)大力發(fā)展多種形式的新型農(nóng)民合作組織。當前合作社的數(shù)量迅猛發(fā)展，根據(jù)國家工商總局的最新統(tǒng)計，截至2013年底，全國依法登記注冊的專業(yè)合作、股份合作等農(nóng)民合作社達98.24萬家，同比增長42.6%；實際入社農(nóng)戶7 412萬戶，約占農(nóng)戶總數(shù)的28.5%，同比增長39.8%。各級示范社超過10萬家，聯(lián)合社達到6 000多家。西部及重慶市的合作社發(fā)展速度更為迅猛，2013年重慶市農(nóng)民專業(yè)合作社達到19 271個，同比增長62%。

合作社從數(shù)量及規(guī)模上逐漸發(fā)展壯大，實現(xiàn)了規(guī)模擴張。但另一方面，一部分合作社發(fā)展過程中存在規(guī)模小、發(fā)展不規(guī)范、管理松散、效率低下等問題，相對于公司、協(xié)會、農(nóng)場等其他類型的農(nóng)村經(jīng)濟組織，合作社發(fā)展時間較短，發(fā)展水平參差不齊，合作社需要增強自身競爭能力，提高運行效率。但當前合作社的評價存在模式不系統(tǒng)、指標不統(tǒng)一、數(shù)據(jù)采集不完整、主觀性強、設(shè)定不合理等問題，難以在各個不同的合作社之間進行有效比較， 繼而提出有針對性的質(zhì)量提升方式。如何有效地對合作社的發(fā)展質(zhì)量進行評價，并規(guī)范合作社的良性有序發(fā)展，在當前合作社的組織形式逐漸多樣化的趨勢下，如何對不同區(qū)域、不同形式、不同經(jīng)營內(nèi)容的合作社進行客觀比較，我們需要對合作社的發(fā)展效率進行評價，并在評價的基礎(chǔ)上提出相應(yīng)的建議。

本文在理論分析的基礎(chǔ)上，利用重慶市的隨機抽樣調(diào)查微觀數(shù)據(jù)，借鑒三階段數(shù)據(jù)包絡(luò)分析(Data Envelopment Analysis，DEA)方法，采用徑向與非徑向結(jié)合的EBM(Epsilon-Based Measure)方向距離函數(shù)，采用隨機前沿分析(Stochastic Frontier Analysis，SFA)方法過濾掉環(huán)境因素及隨機因素，得到準確的合作社效率評價數(shù)據(jù)，在此基礎(chǔ)上，利用自助DEA(Bootstrap DEA)方法，最終得到穩(wěn)健的合作社效率評價結(jié)果。本文的研究結(jié)論不僅對重慶市的合作社發(fā)展能起到積極的指導(dǎo)作用，對西部以及全國的合作社的發(fā)展也有一定的指導(dǎo)意義。本文的結(jié)構(gòu)安排如下：第二部分對所采用的EBM方向距離函數(shù)、三階段DEA及Bootstrap方法進行簡介，并對比其與傳統(tǒng)的DEA方法、徑向及非徑向距離函數(shù)的優(yōu)勢；第三部分是采用的數(shù)據(jù)及實證分析過程，并對不同區(qū)域、不同類型及不同經(jīng)營內(nèi)容的合作社效率進行對比及分解；最后得出相應(yīng)的結(jié)論并提出建議。

二、研究思路與方法

經(jīng)典 DEA 模型假定 決策單元(Decision Making Unit，DMU) 具有同質(zhì)性，即同類型 DMU 應(yīng)具相同的目標和任務(wù)、外部環(huán)境、輸入輸出指標，現(xiàn)有DEA 研究成果大都基于該假設(shè)。但現(xiàn)實中，學(xué)者發(fā)現(xiàn)評價對象并不必然同質(zhì)。人們常忽略異質(zhì)性的存在，如果仍套用經(jīng)典DEA 模型，就會導(dǎo)致結(jié)果出現(xiàn)偏差。同質(zhì)性假設(shè)下，效率分數(shù)能反映DMU 轉(zhuǎn)化投入為產(chǎn)出的能力，異質(zhì)性假設(shè)下，效率分數(shù)夾雜著環(huán)境等異質(zhì)因素對這一能力的影響。因此 DMU 異質(zhì)背景下，需要對經(jīng)典 DEA 模型進行改進以符合現(xiàn)實需要，從而提高 DEA 的適用性、有效性、合理性。

為便于具有異質(zhì)性的各DMU之間的效率比較，我們從距離函數(shù)設(shè)定、DMU所處環(huán)境調(diào)整、效率真實值的測量等角度構(gòu)建三階段Bootstrap EBM模型。

(一)EBM距離函數(shù)

經(jīng)典的 DEA 模型可以劃分為兩類，第一類以徑向測算(radial measure)為基礎(chǔ)的規(guī)模報酬不變(Constant Returns to Scale，CRS) 模型為代表，第二類以非徑向測算(non-radial measure) 為基礎(chǔ)的 SBM(slack-based measure) 模型。但不論傳統(tǒng)DEA徑向模型還是 SBM 模型都存在一定缺陷。為了有效解決 CRS 模型和 SBM 模型測算效率分值存在的問題，Tone(2010b)構(gòu)建了一個綜合徑向和非徑向特點的 EBM(Epsilon-Based Measure)模型[1-2]。對于具有m個投入要素(x)和s個產(chǎn)出(y)的n個決策單元，EBM模型可以表示為：

s.t.θx0-Xλ-s-=0，

λY≥y0，λ≥0，s-≥0.

(二)BootstrapDEA

Kniep(2003)指出，標準的DEA效率值相對于絕對效率水平的估計來說，是有偏的、不一致的估計量[3]。因此，有必要進一步引進Bootstrap方法[4]。Bootstrap方法是由Efron(1979)發(fā)展起來的一種基于數(shù)據(jù)的模擬方法。其基本思想是僅利用樣本數(shù)據(jù)并借助計算機來進行重抽樣，每一次抽樣可以得到一個Bootstrap樣本，然后根據(jù)Bootstrap樣本來計算統(tǒng)計量和估計總體的分布特征，而不對模型做任何假設(shè)[5]。

SimarandWilson(1998；2000)將Bootstrap方法引入到DEA分析中[6]，其基本步驟為：對于第二階段SFA調(diào)整后的投入與產(chǎn)出變量，利用有放回的抽樣技術(shù)，從初始樣本中抽取一個規(guī)模為N的naivebootstrap樣本，并進行平滑化處理，計算DEA模型，得到BootstrapDEA估計值，并重復(fù)m(m>500)次抽樣，最終得到每個DMU的初始效率值偏誤與偏誤修正后的效率得分，并可利用Bootstrap值構(gòu)建真實值的置信區(qū)間。BootstrapDEA的詳細算法步驟見Simar&Wilson(2000)[7]、Kneipet.al(2008)[8]、龔鋒(2008)[9]及扶玉枝(2011)[10]。Wilson(2006)指出，Bootstrap方法是彌補傳統(tǒng)DEA方法缺陷的唯一可行方法。

(三)三階段BootstrapEBM

為保留DEA法的優(yōu)點，同時又能分理出環(huán)境變量等的影響，從而更精確地估計出不同生產(chǎn)單元的效率，F(xiàn)ried等(2002)提出了一種新的效率評價模型——三階段DEA法[11]。該方法最大的特點是能夠?qū)⑼獠凯h(huán)境與隨機誤差對效率的影響去除，使得所計算出來的效率值能更加真實地反映實際情況。三階段DEA法的基本思路與步驟詳見Fried等(2002)。

本文將EBM距離函數(shù)和Bootstrap方法嵌入到三階段DEA法中，構(gòu)建三階段BootstrapEBM法。一階段運用EBM模型得到各樣本單位的初始效率值與投入(或產(chǎn)出)松弛值；二階段運用隨機前沿分析(SFA)方法分解第一階段的投入松弛值，調(diào)整投入或產(chǎn)出值，以排除外部環(huán)境和隨機誤差因素的影響；三階段則將第二階段得出的調(diào)整后的投入值與原始產(chǎn)出值，再次利用EBM距離函數(shù)，利用Bootstrap方法計算各決策單位的效率。

因此，相對于傳統(tǒng)DEA方法，三階段BootstrapEBM具有以下優(yōu)勢：(1)利用參數(shù)的SFA方法調(diào)整了外部環(huán)境與隨機因素的影響，使得各決策單元均調(diào)整至同樣的環(huán)境下，這樣得到的效率值反映了決策單元自身可控的因素所得的效率，能更加真實公正地反映各決策單位的實際運行狀況，也使得各決策單元自身有能力逐步調(diào)整至前沿面*很顯然，外部環(huán)境與隨機因素不是決策單元能夠控制和調(diào)整的。。(2)利用徑向與非徑向距離函數(shù)結(jié)合的EBM距離函數(shù)，避免了傳統(tǒng)DEA方法與SBM方法所得效率值改進方向與現(xiàn)實不吻合的缺點，利用EBM方向距離函數(shù)可以得到更為貼近現(xiàn)實的改進方向與效率評價值。(3)利用BootstrapDEA方法通過重復(fù)抽樣來模擬數(shù)據(jù)生成過程，并且在模擬樣本中應(yīng)用原始估計值，可以近似地得到原始估計值的樣本分布，克服了極端值對效率值變化的影響，使得所得效率值具有穩(wěn)健性，而傳統(tǒng)DEA值極容易受到某個極端值是否進入決策單元的影響。(4)傳統(tǒng)DEA方法極易得到處于相同前沿面的多個單元，這樣無法進行完整排序，而BootstrapDEA法由于進行了大量的重復(fù)抽樣，兩個決策單元效率相同的概率極小，

這樣也有利于對各決策單元進行完整評價。

三、數(shù)據(jù)和變量

(一)數(shù)據(jù)來源與說明

本文的實證數(shù)據(jù)來自于重慶市農(nóng)委2011年進行隨機抽樣調(diào)查的共計95家合作社，數(shù)據(jù)來源準確和完整，樣本地域分布合理，具有代表性。利用三階段BootstrapEBM方法，準確地對合作社效率進行評價，這樣所挑選出的合作社一方面具有較高的執(zhí)行能力與生產(chǎn)效率，可以很好地利用資源，發(fā)展成為更為高效的合作社。另一方面，通過效率比較與分析，可以得到重慶市不同類型、不同區(qū)域的合作社效率差異，通過分析其松弛變量及進行效率分解，有利于效率不高的合作社尋找到提高效率的途徑。同時，這95家合作社代表了典型的重慶市合作社發(fā)展水平，以重慶市的合作社發(fā)展作為評價，對于我國西部地區(qū)及全國的合作社發(fā)展也能起到積極的參考。

(二)投入與產(chǎn)出變量

借鑒扶玉枝(2011)[10]并結(jié)合Ariyaratneetal.(2000)[11]及Galdeano-Gomez(2008)[12]等對合作社效率評價時投入與產(chǎn)出指標的選擇，我們采取入社員數(shù)、資產(chǎn)總額、基地面積作為合作社的各項投入。合作社社員人數(shù)作為合作社的人事投入，合作社資產(chǎn)總額作為合作社的財產(chǎn)投入，基地面積作為合作社的物質(zhì)投入。這三個指標綜合反映了合作社人財物的綜合投入[13]。

合作社的產(chǎn)出績效不僅包括合作社的直接經(jīng)濟效率，趙佳榮(2010)指出，合作社還承擔了相應(yīng)的社會責(zé)任，尤其是帶動其他農(nóng)民致富的重要責(zé)任[14]。因此，經(jīng)濟績效指標上，我們采用年銷售收入與年凈利潤來刻畫合作社的經(jīng)濟規(guī)模與經(jīng)營績效[15]。采用帶動農(nóng)戶數(shù)及增加農(nóng)民純收入兩個指標來刻畫社會績效[16]。 各變量的描述性統(tǒng)計如表1所示。

表1　合作社投入產(chǎn)出描述性統(tǒng)計表

(三)環(huán)境變量

根據(jù)Simar和Wilson(2007) 的“分離假設(shè)”[3]，環(huán)境變量應(yīng)選取對農(nóng)民合作社投入產(chǎn)出效率產(chǎn)生影響但不在樣本主觀可控范圍的因素[12]。具體到影響合作社的環(huán)境變量選擇，這主要包括國家的宏觀經(jīng)濟環(huán)境、政府對農(nóng)業(yè)發(fā)展的相關(guān)政策、人力資源因素等。結(jié)合前述文獻綜述中的研究，本文所采用的環(huán)境變量有：

(1)第一產(chǎn)業(yè)總產(chǎn)值。當?shù)氐牡谝划a(chǎn)業(yè)總產(chǎn)值總量對合作社的生產(chǎn)和銷售產(chǎn)生影響，當?shù)剞r(nóng)業(yè)總產(chǎn)值越大，合作社的原材料來源越為容易。另一方面，較大規(guī)模的本地農(nóng)業(yè)產(chǎn)值為合作社提供了本地化市場，合作社越能夠方便地進行本地生產(chǎn)和銷售，節(jié)約交易與運輸成本，加快資金流轉(zhuǎn)等，提高合作社經(jīng)營效率。

(2)農(nóng)業(yè)商品率。本地的農(nóng)業(yè)商品率越高，說明本地農(nóng)村的市場氛圍越濃，作為農(nóng)民合作主體的合作社參與到市場交換的可能性越大，合作社產(chǎn)品的生產(chǎn)、加工、銷售的容易程度也隨之增加。

(3)財政支農(nóng)支出。在政府對農(nóng)業(yè)發(fā)展的相關(guān)政策方面，財政政策是對合作社發(fā)展的重要政策之一，而財政支農(nóng)支出的增加對合作社生產(chǎn)效率的提高具有積極作用，因此本文選擇各區(qū)縣財政支農(nóng)支出來衡量政府對農(nóng)業(yè)發(fā)展的政策。財政支農(nóng)能夠促進合作社生產(chǎn)的積極性，進而促進農(nóng)業(yè)發(fā)展，提升合作社效率。

(4)農(nóng)村居民人均純收入。隨著農(nóng)民收入的提高，其對合作社的認識也將逐步加深，同時，多項研究表明，收入也是個人能力的表征，由于我們未能獲取分區(qū)縣農(nóng)民的平均受教育程度的信息，因此，我們使用農(nóng)村居民人均純收入作為人力資源質(zhì)量因素的代理變量。

(5)農(nóng)業(yè)人口比例。農(nóng)民專業(yè)合作社作為勞動密集型的機構(gòu)，生產(chǎn)過程需要大量的勞動力資源，因此，是否擁有豐富的本地勞動力資源是制約很多專業(yè)合作社發(fā)展壯大的因素，我們使用農(nóng)業(yè)人口比例來度量人力資源的數(shù)量。但同時，農(nóng)業(yè)人口比例過大，說明該區(qū)域城鎮(zhèn)化與經(jīng)濟發(fā)展水平較低，可能不利于合作社籌集資本。

環(huán)境變量是對合作社這一微觀主體所不能主動調(diào)整和控制的宏觀經(jīng)濟社會環(huán)境的刻畫，為盡量保證對合作社所處環(huán)境的準確描繪，我們使用了縣級數(shù)據(jù)*根據(jù)我們的調(diào)查，當前合作社的生產(chǎn)范圍很少超過縣域。，以上各環(huán)境變量值均來自于《重慶統(tǒng)計年鑒》中對重慶38個區(qū)縣的統(tǒng)計數(shù)據(jù)。

表2　合作社環(huán)境變量描述性統(tǒng)計表

四、實證分析

(一)一階段EBM

傳統(tǒng)DEA模型分析不考慮環(huán)境變量與隨機誤差，以投入與產(chǎn)出變量為基礎(chǔ)，對重慶市95家合作社的效率進行評估，本文考慮產(chǎn)出水平一定的情況下，合作社如何有效地利用所投入的資源，即采用投入導(dǎo)向型的DEA方法，結(jié)合EBM距離函數(shù)，得到以下合作社效率的分布表。

表3　合作社第一階段EBM效率值分布表

在規(guī)模報酬不變情況下，重慶市合作社的技術(shù)效率均值為0.222 2。這說明，合作社的平均技術(shù)效率水平很低，有接近40%的合作社技術(shù)效率低于0.1，8%左右的合作社效率大于0.6，說明在不考慮環(huán)境因素的情況下，合作社之間的效率差異非常大。當前合作社的發(fā)展仍處于起步階段，許多合作社由于管理不善、信息不充分、技術(shù)落后等原因約束了對合作社所投入諸要素的有效利用。

在規(guī)模報酬可變情況下將技術(shù)效率分解為純技術(shù)和規(guī)模效率，從上表中可以看出，純技術(shù)效率大大低于規(guī)模效率。因此，重慶市合作社平均技術(shù)效率較低主要是由于純技術(shù)低造成的?？梢姾献魃鐚ο冗M的農(nóng)業(yè)生產(chǎn)、農(nóng)業(yè)經(jīng)營管理的先進技術(shù)采用不夠，導(dǎo)致純技術(shù)效率低下。

表4　一階段EBM按類型分類的合作社效率值

從分類型的合作社效率分析來看，糧食和一般經(jīng)濟作物類合作社技術(shù)效率相對最高，畜禽水產(chǎn)養(yǎng)殖類的技術(shù)效率值與純技術(shù)效率值相對最低，但畜禽水產(chǎn)養(yǎng)殖類的規(guī)模效率最高。這與扶玉枝(2012)利用浙江省合作社數(shù)據(jù)分析的結(jié)果不完全一致，這主要是因為我們考慮了合作社的社會效益，扶玉枝(2012)的產(chǎn)出指標僅為合作社總收入這一個經(jīng)濟指標*我們使用扶玉枝(2012)的投入產(chǎn)出指標，所得效率值排序與其所研究結(jié)論基本一致。[5]，另一方面，我們還沒有排除環(huán)境因素和隨機因素對效率的干擾?，F(xiàn)實中，畜禽水產(chǎn)養(yǎng)殖類合作社在技術(shù)采用存在上相對較高的進入門檻，而蔬菜水果類的合作社在擴大規(guī)模獲取規(guī)模經(jīng)濟上存在更大障礙。

(二)二階段SFA

由于第一階段的分析是將任何與效率前沿的偏離都看作是管理無效率導(dǎo)致的，沒有考慮不同地區(qū)外部環(huán)境以及隨機誤差的影響，而外部環(huán)境因素和外生隨機因素往往會對一個合作社的產(chǎn)出水平高低有著重要的影響，在這種情況下計算得到的DEA效率得分可能難以真實反映各個合作社的管理效率，以此作為評價標準會有失客觀性與公正性。為了更好地了解掌握各個合作社的真實效率水平，以第一階段得到的各地區(qū)產(chǎn)出松弛值作為被解釋變量、以前文選定的外部環(huán)境因素為解釋變量，運用SFA方法基于半正態(tài)分布假設(shè)對產(chǎn)出松弛值進行分解。

表5　第二階段SFA估計表

注：*、**、***分別表示在10%、5%、1%的水平上顯著。

投入松弛變量是各合作社可以通過改善其運營管理水平，從而達到效率前沿可能節(jié)約的投入量。因而，可以將投入松弛變量視為各地區(qū)的機會成本，解釋變量與投入松弛變量正相關(guān)表明該解釋變量不利于效率的提高；反之，二者負相關(guān)，則表明該解釋變量有利于效率的提高。因此，從環(huán)境因素對投入松弛變量的回歸系數(shù)來看，可得出如下結(jié)論：

(1)農(nóng)業(yè)總產(chǎn)值的擴大有利于效率的提升，且對資產(chǎn)總額與入社社員有顯著影響。隨著農(nóng)業(yè)產(chǎn)值的增加，增加了合作社的各方對農(nóng)業(yè)投資的興趣，對合作社的諸要素投入可以帶來效率的提升，因此需要農(nóng)業(yè)管理部門積極引導(dǎo)各類資本進入農(nóng)業(yè)領(lǐng)域，增加農(nóng)業(yè)投資。

(2)農(nóng)業(yè)商品化率的提高推動了合作社向前沿面邁進，且對資產(chǎn)總額與入社社員顯著有效。隨著商品化率的提高，可以增加農(nóng)民收入，提高了農(nóng)民自身及把其土地入股參加合作社的積極性，也促使農(nóng)民積極加入合作社。

(3)財政支農(nóng)支出對合作社減少投入松弛變量、提升效率起到積極作用，且均顯著有效。這說明重慶市當前對農(nóng)業(yè)生產(chǎn)及合作社的積極支農(nóng)財政政策取得了良好的效果，合作社管理部門通過政策、項目及資金支持對合作社進行幫扶，從人財物各個方面均對合作社產(chǎn)生了重要的促進作用。

(4)農(nóng)村居民人均純收入的增加有利于減少投入松弛變量，提高效率。農(nóng)民收入的增加激勵了農(nóng)民將更多的人財物資源投入到合作社生產(chǎn)經(jīng)營中去，尤其是將資金與土地資本投入到合作社中，推動合作社發(fā)展壯大。

(5)農(nóng)業(yè)人口比例的增加并不利于合作社各類投入松弛變量的減小，且對入社社員松弛變量的增加有顯著的促進作用，說明傳統(tǒng)的農(nóng)業(yè)社會形態(tài)并不一定帶來合作社效率的提升。與之對應(yīng)，城鎮(zhèn)化水平的提高有利于增加合作社的效率值，因此，隨著我國城鎮(zhèn)化水平的進一步提高，合作社的發(fā)展擁有更為廣闊的前景。

(三)三階段BootstrapEBM

在得到經(jīng)過調(diào)整后的各地區(qū)投入項后，把調(diào)整后的投入和原始產(chǎn)出帶入DEA模型，與一階段一致，此處仍然選擇投入導(dǎo)向與EBM距離函數(shù)，利用BootstrapDEA方法，重復(fù)抽樣500次，我們可得剔除了環(huán)境因素和隨機因素影響的重慶市合作社Bootstrap效率值。

表6　合作社第三階段Bootstrap EBM效率值分布表

通過與一階段對比可以看到，在經(jīng)過第二階段的調(diào)整后，合作社的效率評價值發(fā)生了很大變化?？傮w上，合作社的技術(shù)效率與純技術(shù)效率均值有所下降，規(guī)模效率均值所有上升，這說明合作社自身不可控的環(huán)境因素確實對合作社的效率值產(chǎn)生了影響，剔除這些影響因素的效率評價值更能反映合作社自身的能力水平。

表7　三階段按類型分類的合作社效率值比較

從分類型的合作社效率分析來看，與一階段比較，三類合作社的效率值變化不一，說明環(huán)境因素對不同類型合作社的影響各不相同。經(jīng)過調(diào)整后，糧食和一般經(jīng)濟作物類合作社各類效率值較高，蔬菜水果類合作社的平均技術(shù)效率與純技術(shù)效率值均為最低，畜禽水產(chǎn)養(yǎng)殖類合作社的規(guī)模效率平均值最低。結(jié)合了Bootstrap方法后，各類效率值均有所調(diào)整，但各自調(diào)整的大小不一致。在剔除了環(huán)境因素與采用Bootstrap抽樣后，此時得到的效率值更能反映各個合作社的運行情況。

在95%的置信區(qū)間下，我們比較三階段未經(jīng)Bootstrap調(diào)整、經(jīng)過Bootstrap調(diào)整的技術(shù)效率、純技術(shù)效率及規(guī)模效率值，以及各類效率值置信區(qū)間的大小。

圖1　第三階段合作社Bootstrap效率變化圖

從圖1可以看出，Bootstrap方法使得大部分合作社的技術(shù)效率和純技術(shù)效率向下調(diào)整，規(guī)模效率向上調(diào)整，但各個合作社受自助抽樣影響的大小不一，如純技術(shù)效率值在效率值較大處受影響較大，而規(guī)模效率值在效率處于中間值時受影響較大。

同理，我們分析了三種類型的合作社的技術(shù)效率在未經(jīng)Bootstrap調(diào)整、經(jīng)過Bootstrap調(diào)整以及各類效率值置信區(qū)間的大小。

圖2　三類合作社的技術(shù)效率比較圖

從圖2可以看出，Bootstrap方法使得三類合作社的技術(shù)效率值向下調(diào)整，但各個合作社受自助抽樣影響的大小不一，糧食和一般經(jīng)濟作物類與蔬菜水果類合作社調(diào)整較小，而畜禽水產(chǎn)養(yǎng)殖類調(diào)整較大。且各個類型合作社內(nèi)部的調(diào)整幅度也不一致，前兩類合作社在效率值極大處調(diào)整最大，而第三類合作社在效率為中等及較大處調(diào)整最大。

圖3　三類合作社的純技術(shù)效率比較圖

我們同樣分析了三種類型的合作社的純技術(shù)效率在未經(jīng)Bootstrap調(diào)整、經(jīng)過Bootstrap調(diào)整以及各類效率值置信區(qū)間的大小。從圖3可以看出，Bootstrap方法使得三類合作社的純技術(shù)效率值向下調(diào)整，但各個合作社受自助抽樣影響的大小不一，糧食和一般經(jīng)濟作物類與蔬菜水果類合作社調(diào)整較小，而畜禽水產(chǎn)養(yǎng)殖類調(diào)整較大。在同一類型合作社內(nèi)部，前兩類合作社仍然在效率值大處調(diào)整較大，而第三類合作社則是較為均勻的調(diào)整。

圖4三類合作社的規(guī)模效率比較圖

圖4分析了三類合作社的規(guī)模效率在未經(jīng)Bootstrap調(diào)整、經(jīng)過Bootstrap調(diào)整以及各類效率值置信區(qū)間的大小，Bootstrap方法使得三類合作社的規(guī)模效率值向上調(diào)整，且各類合作社受自助抽樣的影響調(diào)整也較為一致，均為在效率處于中上值處調(diào)整較大，但置信區(qū)間的變化較大，三類合作社均在效率為中上值時置信區(qū)間的范圍較寬，說明此時合作社受效率評價樣本的影響較大。

五、結(jié)論及建議

使用隨機抽樣調(diào)查的重慶市95家合作社樣本數(shù)據(jù)，基于徑向與非徑向距離函數(shù)對效率改進測量的不足，我們結(jié)合EBM距離函數(shù)，利用三階段DEA方法剔除了異質(zhì)性合作社的環(huán)境因素與隨機因素對合作社效率評價的影響，在此基礎(chǔ)上，采用Bootstrap方法，我們分析了合作社的效率并對合作社效率進行了充分排序。通過對合作社的評價指標選取、距離函數(shù)識別、環(huán)境因素調(diào)整、真實效率評價，我們提出以下政策建議：

(1) 衡量合作社的效率應(yīng)考慮經(jīng)濟效益與社會效益相結(jié)合。作為農(nóng)村中帶動農(nóng)民致富的重要載體，合作社不僅應(yīng)考慮作為經(jīng)濟組織的經(jīng)濟效益，另一方面，合作社帶動其他農(nóng)戶致富也應(yīng)成為合作社考核中的重要指標。當前重慶市合作社效率不高的一個重要原因在于合作社發(fā)揮社會效益的功能仍不夠強大，因此，對合作社進行評價時應(yīng)加大對合作社社會功能的考核，引入帶動農(nóng)戶數(shù)及農(nóng)戶增收金額等合作社社會效益指標。

(2)環(huán)境因素對效率值有顯著影響，政府應(yīng)為合作社的發(fā)展創(chuàng)造良好的外部環(huán)境。通過二階段使用環(huán)境因素對投入松弛變量的影響分析，可以看到，各個環(huán)境因素對不同的投入變量影響不盡相同，第二階段對五類環(huán)境變量的SFA分析，我們得出，政府應(yīng)積極引導(dǎo)各類資本進入農(nóng)業(yè)領(lǐng)域，提高農(nóng)產(chǎn)品的商品化率，加大政府財政支農(nóng)支出的力度及對合作社的幫扶，并引導(dǎo)農(nóng)民在增收的同時理性進行農(nóng)業(yè)生產(chǎn)投入。

(3)政府進行效率評價應(yīng)采用多次重復(fù)抽樣的方式得到穩(wěn)健的效率評價值。對合作社效率評價時，應(yīng)充分認識到單次效率評價所得結(jié)果對樣本選擇的敏感性，由于效率評價的相對性，某個評價對象的進入或退出將可能對全體樣本產(chǎn)生影響，從而使得所得效率值不夠穩(wěn)健與真實。采用多次重復(fù)抽樣的評價方式，有效地避免了這一缺陷，通過Bootstrap效率評價結(jié)果，可以看出，重慶市合作社當前的整體效率值較低，尤其是對于蔬菜水果類合作社，通過環(huán)境變量調(diào)整后為效率最低類合作社，政府應(yīng)提高合作社的技術(shù)采用能力、加強技術(shù)利用效率、加大農(nóng)業(yè)技術(shù)推廣，擴大合作社規(guī)模、降低入社門檻，通過增大其純技術(shù)效率及規(guī)模效率，實現(xiàn)此類合作社效率的提升。

[1]TONE,K.,TSUTSUI,M.Anepsilon-basedmeasureofefficiencyinDEA-Athirdpoleoftechnicalefficiency[J].EuropeanJournalofOperationalResearch, 2010,207(1), 1554-1563.

[2] 師博，沈坤榮.城市化、產(chǎn)業(yè)集聚與EBM能源效率[J].產(chǎn)業(yè)經(jīng)濟研究，2012(6)：76-81.

[3]KNEIP,A.,et.al.AsymptoticsforDEAEstimatorsinNonparametricFrontierModels[EB/OL]. [2003-08-11].http://ideas.repec.org/a/cup/etheor/v24y2008i06p1663-1697_08.html.

[4]EfronB.Bootstrapmethods：anotherlookatthejackknife[J].TheAnnalsofStatistics，1979：1-26.

[5] 王亞華，吳凡，王爭.交通行業(yè)生產(chǎn)率變動的BootstrapMalmquist指數(shù)分析(1980—2005)[J].經(jīng)濟學(xué)(季刊)，2008(3)：15-29.

[6]SIMAR,L.,WILSONP.Estimationandinferenceintwo-stage,semi-parametricmodelsofproductionprocesses[J],JournalofEconometrics,2007, 136 (1) 31-64.

[7]SimarL，WilsonPW.Sensitivityanalysisofefficiencyscores：Howtobootstrapinnonparamtricfrontiermodels[J].Managementscience，1998，44(1)：49-61.

[8]KNEIP,A.,et.al.AsymptoticsandConsistentBootstrapsforDEAEstimatorsinNonparametricFrontierModels[J].EconometricTheory, 2008,24(1): 1663-1697.

[9] 龔鋒.地方公共安全服務(wù)供給效率評估——基于四階段DEA和BootstrappedDEA的實證研究[J].管理世界，2008(4)：80-90.

[10] 扶玉枝.農(nóng)業(yè)合作社效率研究——基于目標函數(shù)的靜態(tài)與動態(tài)分析[D].浙江大學(xué)，2012：1-30.

[11]ARIYARATNEC.B.,et.al.MeasuringX-EfficiencyandScaleEfficiencyforaSampleofAgriculturalCooperatives[J].AgriculturalandResourceEconomicsReview, 2000, 29(2):198-207.

[12]GALDEANOE.ProductivityEffectsofEnvironmentalPerformance:EvidencefromTFPAnalysisonMarketingCooperatives[J].AppliedEconomics, 2008,40(14): 1873-1888.

[13] 劉自敏，楊丹.基于成員異質(zhì)性的農(nóng)民股份合作社收益分配研究——雙邊專用性投資的視角[J].南京農(nóng)業(yè)大學(xué)學(xué)報：社會科學(xué)版，2014(1)：59-67.

[14] 趙佳榮.農(nóng)民專業(yè)合作社“三重績效”評價模式研究[J].農(nóng)業(yè)技術(shù)經(jīng)濟,2010(2):35-41.

[15] 黃季焜，鄧衡山，徐志剛.中國農(nóng)民專業(yè)合作經(jīng)濟組織的服務(wù)功能及其影響因素[J].管理世界，2010(5)：75-81.

[16] 楊丹，劉自敏，徐旭初.環(huán)境異質(zhì)性，合作社交叉效率與合作關(guān)系識別[J].農(nóng)業(yè)技術(shù)經(jīng)濟，2015(3)：51-60.

TheEfficiencyEvaluationofHeterogeneousCooperativesBasedonThreeStageBootstrapEBM：CaseofChongqing

YANGDan1，2，HUANGFeng1，LIUZi-min1

(1.CollegeofEconomicsandManagement,SouthwestUniversity，Chongqing400716，China；2.ResearchCenterofAgriculturalEducationandDevelopment,SouthwestUniversity，Chongqing400716，China)

BasedonthedefectsoftraditionalDEAmodelradialandnon-radialdistancefunctionandanalysisofheterogeneitybetweenDMU，webuildanEBMdistancefunctionwhichintegratesthecharacteristicsofradialandnon-radialtoresolvethedeficienciesoftraditionaldistancefunction.Thenweresolveheterogeneityproblembyusingthree-stageDEAmethod.Onthisbasis，weusebootstrapmethodtosolveproblemsofefficiencyvalueswhicharebiasedandinconsistent.Weusecooperativeempiricaldataandanalyzeenvironmentvariableswhichaffectthevalueofinputs，andcomparativeanalysisbasedonevaluationoftheefficiencybeforeandafteradjustingtheenvironmentvariablesandtheefficiencyvaluesandtheirconfidenceintervalsbeforeandafteradjustingbootstrapEBM.Finally，accordingtomoreaccurateandobjectiveefficiencyevaluationvalue，weproposeappropriatepolicyrecommendations.

Heterogeneity；Cooperatives；EBMDistanceFunction；ThreeStageDEA；Bootstrap；ConfidenceIntervals

[DOI]10.16164/j.cnki.22-1062/c.2015.06.014

2013-12-05

國家自然科學(xué)基金重點學(xué)科群項目(71333011)；國家自然科學(xué)基金項目(71373063)；教育部社科規(guī)劃青年項目(13YJC790101)；重慶市社科規(guī)劃青年項目(2012QNJJ011)；中央高校基本科研重點項目(SWU1509120)；重慶市人文社會科學(xué)重點研究基地重點項目(14SKB049)。

楊丹(1981-)，女，湖北宜昌人，西南大學(xué)經(jīng)濟管理學(xué)院副教授，西南大學(xué)農(nóng)業(yè)教育發(fā)展研究中心研究員；黃鳳(1990-)，女，貴州黔西人，西南大學(xué)經(jīng)濟管理學(xué)院碩士研究生；劉自敏(1981-)，男，四川德陽人，西南大學(xué)經(jīng)濟管理學(xué)院副教授，博士。

F224.0

A

1001-6201(2015)06-0073-09

[責(zé)任編輯：秦衛(wèi)波]

*感謝西南大學(xué)經(jīng)濟管理學(xué)院俞鋒同學(xué)優(yōu)秀的助研工作。