基于三階段數據包絡分析模型的竹林生產效率研究

楊水生，徐秀英,2，符椒燕，徐堇寒

（1.浙江農林大學 經濟管理學院，浙江 杭州311300；2.浙江省農民發展研究中心 浙江省哲學社會科學重點研究基地，浙江 杭州311300）

基于三階段數據包絡分析模型的竹林生產效率研究

楊水生1，徐秀英1,2，符椒燕1，徐堇寒1

（1.浙江農林大學 經濟管理學院，浙江 杭州311300；2.浙江省農民發展研究中心 浙江省哲學社會科學重點研究基地，浙江 杭州311300）

林業生產效率是關系林業經濟發展的重要問題。對林業生產效率進行科學地測算與分析，真實反映林業生產效率的實際狀況，以期為提高林業生產效率提供科學依據。以竹林生產為例，基于浙江省安吉縣6個行政村110戶農戶調查數據，應用三階段數據包絡分析（DEA）模型對農戶竹林生產效率進行測算與評價分析。結果表明：環境變量中，戶主受教育年限、竹林收入比例、家庭從事竹林生產人數、家庭經營竹林地塊數對竹林生產效率有顯著的負向影響；戶主年齡、參與竹林生產技術培訓、林權證的獲得對竹林生產效率有顯著的正向影響。在剔除環境效應和隨機誤差影響后，樣本戶竹林生產平均技術效率由0.537下降為0.436，平均純技術效率由0.602上升至0.803，而平均規模效率則由0.891下降至0.556；竹林生產規模報酬均處于遞增階段。基于上述分析，提出降低竹林地細碎化程度，引導農戶合理配置竹林生產投入要素，進一步明晰竹林地產權，加大對農戶竹林生產技術培訓力度等建議。表6參27

林業經濟學；竹林生產效率；三階段數據包絡分析模型；環境效應；隨機誤差；安吉縣

20世紀80年代初期的林業 “三定”（穩定山權林權、劃定自留山、確定林業生產責任制）以來，農戶逐漸成為重要的林業生產經營主體。林業生產效率關系到林業經濟的發展，如何優化林業生產要素的配置，提高現有資源的利用效率，使實際林產品的產出接近潛在產出，并厘清哪些因素會影響林業生產效率，是進一步深化集體林權制度改革的重要舉措。因此研究如何提高林業生產效率對于促進農民增產增收、發展地區經濟具有十分重要的意義。近年來，林業生產效率的研究一直受到國內外學者的重視，并展開了大量的研究。部分學者基于宏觀統計數據，采用數據包絡分析（DEA）模型或隨機前沿分析（SFA）模型對中國各省份林業生產效率進行了測算與評價［1-4］，也有學者對國外木材采伐經營效率進行了分析［5-6］。部分學者基于農戶調查數據，主要采用DEA方法對農戶經營林產品（以杉木Cunninghamia lanceolata或果林為主）的生產效率進行實證分析［7-9］。學者們在對林業生產效率進行測算的基礎上，對林業生產效率的影響因素進行了分析［10-12］，徐秀英等［11］利用改進的C-D（COBB-DOUGLAS）生產函數模型實證分析了林地細碎化對竹林產出的影響。然而傳統DEA方法不能剔除環境效應和隨機誤差的影響［13-14］，難以真實地反映林業生產效率的實際狀況。FRIED等［15］提出的三階段DEA模型可以剔除環境效應和隨機誤差的影響，能夠客觀真實地測算生產效率。為此，部分學者采用三階段DEA模型對福建省林業生產效率進行了評價分析［16-17］，研究發現規模無效是導致效率低下的主要原因。從現有文獻來看，以竹林作為研究對象，基于農戶微觀層面數據進行生產效率分析的較為薄弱，采用三階段DEA模型對竹林生產效率進行分析的則更為少見。竹林資源是世界上重要的森林資源，而中國是世界上最主要的產竹國，浙江省作為中國竹林資源的主要分布區之一，其竹林年總產值居于全國前列，但竹林生產效率并不高［18］。鑒于此，以浙江省安吉縣農戶竹林生產的微觀數據為例，采用三階段DEA模型，期望對竹林生產效率進行更為準確的測算；另外，找出主要影響竹林生產效率的因素，為政府部門決策提供依據。

1 數據來源及研究方法

1.1 數據來源

安吉縣隸屬于浙江省湖州市，常住人口46萬人，總面積為1 886 km2，森林覆蓋率71.1%，林地面積13.83萬hm2，其中竹林面積7.20萬hm2，有 “中國竹鄉”的美譽。安吉縣有毛竹Phyllostachys edulis林5.66萬hm2，占林地總面積40.93%，竹林面積的78.61%。毛竹現存量為1.8億株，采伐3 000萬株·a-1，生產竹筍5.6萬t·a-1。因此，選擇安吉縣的農戶調查數據來研究竹林生產效率具有代表性。

本研究所使用的數據來源于2015年7月對浙江省安吉縣3個鄉鎮6個行政村的農戶調查，抽取的3個鄉鎮（孝豐鎮、杭垓鎮、報福鎮）為縣內竹林資源較為豐富。樣本村和樣本農戶均采取隨機抽樣的方法抽取，選取行政村2個·鄉鎮-1，共6個行政村，選取農戶20個·村-1，共調查農戶120戶，剔除無經營竹林生產的樣本和無效樣本后，共獲得有效問卷110份，有效率91.67%。調查內容主要包括農戶在2014年的家庭基本情況、農戶家庭經營竹林地的資源特征及竹林地塊2013-2014年的投入產出狀況等。

1.2 研究方法

數據包絡分析（DEA）模型［19］自 1978年提出以來，已經由最初的 CCR（CHARNES，COOPER和RHODES）模型發展到現在的幾十種擴展模型。其中學者們針對傳統DEA模型沒有考慮環境因素影響的缺陷，提出了包括兩階段DEA法、三階段DEA法以及四階段DEA法的調整方法［20］。其中三階段DEA模型由于剔除了外部環境等非經營性因素對效率的影響，能夠更加準確地測算生產效率［15］。

第1階段：傳統DEA模型（BCC模型）。選取農戶作為決策單元，由于農戶更加容易控制生產過程中的投入要素，因此，該階段采用以一定產出水平下投入最小化為目標的投入導向BCC（BANKER，CHARNES和COOPER）模型。該模型已經非常成熟，不再進行詳細描述。

第2階段：類似SFA模型。第1階段BCC模型計算得到竹林面積、資本投入、勞動力投入的松弛變量值，由于3種投入松弛變量受環境效應、隨機誤差和管理無效率的影響，通過建立類似SFA模型分別觀察上述3種因素對松弛變量的影響，以此調整投入要素，使得決策單元處于相同的環境和運氣下。建立類似SFA模型：

式（1）中：sij表示第j個決策單元第i項投入的松弛變量；fi（zj;βi）表示環境效應對投入松弛變量的影響；vij+uij是混合誤差項，其中vij表示隨機誤差，假設其服從N（0，σ2vi）分布，而uij表示管理無效率，假設其服從截斷正態分布 N+（ui， σ2ui）， uij與 vij相互獨立不相關。 另外， 公式 γ＝σ2ui/（σ2vi+σ2ui）中， γ 為待估參數，表示混合誤差項中管理無效率所占的比例，其取值介于0到1之間。當γ趨向于1時，表明生產函數的誤差主要是由管理因素引起；而當γ趨向于0時，表明生產函數的誤差主要是由隨機誤差引起。

其次，在 “成本函數”隨機前沿模型中，計算管理無效率的估計公式［21］如下：

式（2）中： εj＝vij+uij， σ2＝σu2+σv2， λ＝σu+σv。 結合式（1）和式（2）， 推導出隨機誤差的計算公式：

式（3）中： sij可由第 1 階段 DEA 模型求出， fi（zj;是可觀測的環境變量，uij｜vij+uij］可由式（2）求得。

在此基礎上，對各決策單元的投入進行調整，剔除環境效應和隨機誤差的影響，使其面臨相同的環境和運氣。基于最有效的決策單元，給出計算調整各投入量的公式：

式（4）中：xij為生產過程中的初始投入值，x＾ij為調整后的投入值， β＾i和v＾ij分別表示環境變量參數和隨機干擾項的估計值。

第3階段：調整后的DEA模型。針對調整后的投入值與原產出值，使用BCC模型，重新測算各決策單元的生產效率。

2 變量選取及描述統計

2.1 投入產出變量選取及描述統計

本研究以調查農戶為決策單元，考慮到竹林生產存在明顯的大小年現象，投入產出變量均采用樣本農戶經營竹林地的2013-2014年的平均值。其中竹林產出變量為竹材產值和竹筍產值的加總；竹林投入變量選取包括土地投入、資本投入和勞動力投入［17］。土地投入指農戶經營的竹林地總面積；資本投入包括化肥、灌溉、竹筍覆蓋物費用及其他物質費用；勞動力投入包括竹林地日常管護、挖筍及竹材采伐的自投工及雇工，按8 h·d-1折算成工日。農戶家庭竹林生產投入產出情況見表1。

DEA模型要求各投入變量與產出變量之間必須符合 “同向性”原則，一般采用Pearson相關性檢驗方法對其進行檢驗［22］。檢驗結果如表2所示。結果顯示：竹林面積、資本和勞動力的投入與竹林總產值之間的相關系數均在1%的水平上顯著為正，表明各投入變量與產出變量之間滿足 “同向性”原則，具有合理性。

2.2 環境變量選取及描述統計

環境變量是指影響生產效率的外部客觀環境因素［22］。結合已有研究［23-24］，環境變量選取分別為農戶人力資本特征變量、家庭經營特征變量、竹林地細碎化變量以及相關政策因素等。考慮到生產經營決策直接影響生產效率的高低，農戶家庭的決策多數是由戶主來決定，因此農戶人力資本特征變量選取戶主年齡及受教育年限來衡量［16］。家庭經營特征變量選取竹林收入比例、家庭從事竹林生產人數來衡量。一般竹林收入占比及從事竹林生產人數越多的農戶，其生產的專業性越強。根據斯密的分工理論，專業化可以帶來效率。竹林地細碎化變量選取農戶家庭竹林地塊數來衡量。一般來說竹林地細碎化會降低竹林生產的規模效應，影響竹林生產效率［11］。選取是否獲得林權證來衡量政府對林業發展的政策因素［16］。是否獲得林權證主要表現為集體林權制度改革后林地的地權穩定性，地權越穩定，越有利于農戶對竹林生產進行長期投資。此外，選取是否參與竹林生產技術培訓作為配套集體林權制度改革服務變量［25］。是否參與竹林生產技術培訓表現為政府在指導林業經濟發展過程中提供的技術支持，技術培訓一般有利于竹林生產技術的更新，是促進生產效率提高的重要因素。環境變量的說明及描述統計見表3。

表1 農戶家庭竹林生產投入產出情況Table 1 Status of bamboo production of households’input and output

表2 竹林投入與產出變量的Pearson相關系數Table 2 Pearson correlation coefficient of household input and output variables

表3 環境變量的說明及描述統計Table 3 Interpretation and statistical description of environment variables

由表3可知：農戶戶主平均年齡為53.63歲，表明從事竹林生產的農戶戶主年齡較大；戶主平均受教育年限為6.38 a，文化程度主要集中在小學、初中；農戶家庭竹林收入占家庭總收入的比例平均為19%，平均每個家庭從事竹林生產的人數為1.37個，表明調研地區農戶竹林生產在家庭經濟中占有一定的地位；農戶家庭經營的竹林地塊數戶均為3.61塊，竹林生產的細碎化程度較高。另外，樣本農戶竹林地林權證獲得率為91%，還有少部分農戶竹林地沒有拿到林權證，表明調研地區林權主體改革尚未徹底。家庭有成員參與竹林生產技術培訓的農戶僅占樣本總數的25%，表明重要的林權配套改革——技術服務開展程度較低。

3 實證結果分析

3.1 第1階段DEA實證結果

利用DEAP 2.1軟件對110個樣本農戶的竹林生產效率進行測算［26］，竹林生產各效率值測算結果如表4所示。

由表4可以看出：樣本地區竹林生產效率整體水平不高，平均技術效率為0.537，離最優效率值1還有較大的差距。如果消除技術效率損失，竹林生產存在46.30%的提升空間。而竹林生產平均純技術效率和平均規模效率分別為0.602和0.891，竹林生產規模效率更接近生產前沿面，表明純技術效率無效是導致生產效率低下的主要原因。在110戶農戶樣本中，僅有9戶農戶處于竹林生產有效率狀態（技術效率為1，且純技術效率與規模效率均為1），僅占到了樣本農戶總數的8.18%。竹林生產技術效率、純技術效率處于無效率程度中等及以下的戶數分別有84戶、77戶，占樣本農戶總數的76.36%和70.00%。而規模效率為無效率程度輕微和有效率的有99戶，占樣本總數的90.00%。

表4 農戶竹林生產效率分布Table 4 The distribution of the bamboo production efficiency of households

3.2 第2階段SFA回歸結果

由第1階段得到了各松弛變量值，通過運用Frontier 4.1軟件，將其分解為環境效應、隨機誤差和管理無效率。各松弛量作為因變量，將7個環境變量作為自變量，建立類似SFA模型，得到回歸結果如表5所示。

表5 第2階段SFA估計結果Table 5 The result of SFA regression in the second stage

由表5可知：竹林面積、資本和勞動力3種投入松弛變量的γ值都接近于1，且均達到了1%的顯著水平，表明管理因素對3種投入要素的影響占據主導地位，對竹林生產效率存在顯著的影響；同時還說明管理無效率是混合誤差項變異的主要影響因素，驗證了運用SFA模型的合理性。

進一步分析環境因素對竹林生產投入松弛變量的影響，當回歸系數通過統計上的顯著性檢驗，且回歸系數為正（負）時，表明增加該變量值不利于（有利于）生產技術效率的提高。具體分析結果如下：

第一，農戶人力資本特征變量。戶主年齡對竹林面積松弛變量的系數為負，并在1%的統計水平上通過顯著性檢驗。表明隨著戶主年齡的增長，竹林面積松弛量會減少，可能是因為年齡較大的農民在干中學積累了豐富的生產經驗，有利于降低竹林面積投入要素的資源浪費，進而對竹林生產效率產生有利影響。戶主受教育年限對竹林面積松弛變量、勞動力投入松弛變量的系數均在1%的水平上顯著為正。結合樣本特點，88.18%的戶主受教育年限在9 a及以下，農戶戶主文化水平普遍較低，而受教育年限較多的戶主往往傾向于向城市和非農產業轉移，缺乏竹林生產的經驗和技術，并疏于竹林地的經營管理，抑制了農戶竹林生產效率的提高。

第二，家庭經營特征變量。竹林收入占比對資本投入松弛變量、勞動力投入松弛變量的系數為正，分別在1%，10%的水平上顯著，表明竹林收入比例增加時，資本投入松弛變量、勞動力投入松弛變量也會增加。家庭從事竹林生產人數對竹林面積松弛變量、勞動力投入松弛變量的系數為正，且均通過1%的顯著性檢驗。表明隨著農戶家庭從事竹林生產勞動力數量的增加，竹林面積松弛變量、勞動力投入松弛變量也會增加。可能是因為農戶家庭在竹林生產中投入越多勞動力以及竹林收入占比越大時，表明農戶生產經營較為單一、對竹林生產的依賴性越大，這就使得農戶存在較高的生產風險，高風險可能帶來農戶生產的過度投入［27］。從而造成投入資源的浪費，對竹林生產效率產生不利的影響。

第三，竹林地細碎化變量。竹林地塊數對3種投入松弛變量的系數均為正，且均通過1%的顯著性檢驗，表明農戶家庭經營竹林地塊數增多時，3種投入松弛量將會增加，導致3種投入要素的過度浪費，對竹林生產效率產生不利的影響。一方面是因為竹林地細碎化程度越高，就越不便于先進機械和技術的推廣應用，難以實現規模經營；另一方面，家庭經營竹林地塊數越多，越不利于統一經營管理，造成農戶在不同地塊之間奔波時間的增加，在一定程度上造成資源浪費，從而影響竹林生產效率的提高。

第四，相關政策因素。家庭成員是否有人參與竹林生產技術培訓對3種投入松弛變量的系數均為負，且均通過1%的顯著性檢驗；是否獲得林權證對3種投入松弛變量的系數也均為負，竹林面積松弛變量、資本投入松弛變量均在1%的水平上顯著。表明農戶家庭有成員參與竹林生產技術培訓或者竹林地有林權證，有利于3種投入松弛變量的減少。說明林業生產技術培訓作為集體林權制度改革的配套服務有利于竹林生產效率的提高，林權證的發放使得產權更加明晰和安全，增強了農戶進行竹林生產的積極性，并進行生產要素的合理配置，進而減少了3種投入要素的浪費。

3.3 調整后的DEA實證結果

根據式（4）得到調整后的竹林生產投入要素值，將初始產出值與調整后的投入值再次進行BCC模型測算，最后得到竹林生產各效率值及其區間分布，如表6所示。

表6 調整的農戶竹林生產效率分布Table 6 Distribution of the bamboo production efficiency of households

對比表4和表6可以看出：在剔除環境效應和隨機誤差對效率影響后的同質環境下，第3階段的竹林生產平均技術效率從第1階段的0.537下降為0.436，下降幅度為18.81%；平均純技術效率由0.602上升為0.803，上升幅度為33.39%；平均規模效率由0.891下降為0.556，下降幅度達37.6%。調整之后的純技術效率更為接近效率前沿面，規模效率低下成為農戶竹林生產效率低下的主要制約因素。在110戶農戶樣本中，僅有3戶農戶的竹林生產處于有效率狀況，比調整前降低66.67%。竹林生產技術效率、純技術效率處于無效率程度中等及以下的農戶比例分別為84.54%和29.09%，分別比第1階段高8.17個百分點和低40.91個百分點，而規模效率為無效率程度輕微和有效率占比35.46%，比第1階段低54.54個百分點。調整后的規模效率大幅度下降，表明農戶竹林生產的實際規模效率并沒有那么高；而調整后的純技術效率大幅度上升，表明調整前農戶竹林生產的純技術效率較低是由于較差的外部環境導致的。另外，農戶竹林生產規模報酬均處于遞增階段，農戶擴大竹林生產規模，各種生產要素投入量增加1倍將獲得大于1倍的產出增加，表明樣本地區農戶竹林生產效率著重改進的方向為提高規模效率。

4 結論及建議

本研究基于浙江省安吉縣農戶竹林生產微觀數據，采用三階段DEA模型，對竹林生產效率及其影響因素進行了分析。得出以下主要結論：①第1階段測算結果顯示：竹林生產平均技術效率為0.537，平均純技術效率為0.602，平均規模效率為0.891。表明農戶竹林生產效率整體水平較低，純技術效率無效是導致生產效率低下的主要原因。在現有生產水平下，若消除技術無效率，竹林的產出水平還能提高46.3%。②第3階段測算結果顯示：竹林生產平均技術效率為0.436，平均純技術效率為0.803，平均規模效率為0.556。表明當剔除環境效應和隨機誤差的影響后，即所有農戶面臨同質的環境和運氣時，竹林產出水平還能提高56.4%。同時，對比調整前后的技術效率、純技術效率及規模效率，都有較大幅度的變化，技術效率的下降主要受規模效率的大幅度下降的影響，而純技術效率有一定程度的上升。③通過第2階段的SFA回歸分析發現，環境變量中，戶主受教育年限、竹林收入比例、家庭從事竹林生產人數、家庭經營竹林地塊數是竹林生產效率提高的不利因素；戶主年齡、參與竹林生產技術培訓、林權證的獲得是竹林生產技術效率提高的有利因素，有利于生產要素的優化配置。

根據以上結論，建議：第一，降低竹林地細碎化程度。鼓勵農戶通過林地的轉讓、出租、入股等流轉方式，將林地向林業專業大戶、家庭林場等現代林業經營主體集中，實現林地規模經營，并降低竹林地的細碎化程度，進一步提高竹林生產效率。第二，引導農戶合理配置竹林生產要素。在增加竹林生產規模的同時，應注意不能盲目地擴大生產投入，當地政府應該積極引導農戶合理分配竹林生產各投入要素，特別是應采用機械等先進技術設備來代替林業勞動力，合理優化資源配置，減少資源浪費。第三，進一步明晰竹林地產權。新一輪集體林權制度改革的林權主體改革尚未徹底，當地政府應加強勘界發證工作，確保林權證發放到農戶，提高農戶林業生產的積極性。第四，加大對農戶竹林生產的技術培訓力度。通過引進先進的竹林生產技術和管理方法，加強技術培訓和指導，提高農戶林業生產的技術水平，促進林業生產效率的提高。

［1］ 李春華，李寧，駱華瑩，等.基于DEA方法的中國林業生產效率分析及優化路徑［J］.中國農學通報，2011， 27（19）： 55-59.LI Chunhua,LI Ning,LUO Huaying,et al.The efficiency analysis and path optimization of forestry input-output in China based on data envelopment analysis ［J］.Chin Agric Sci Bull,2011,27（19）:55-59.

［2］ 田淑英，許文立.基于DEA模型的中國林業投入產出效率評價［J］.資源科學，2012，34（10）：1944-1950.TIAN Shuying,XU Wenli.Evaluation of China’s forestry input-output efficiency based on DEA modeling ［J］.Resour Sci,2012,34（10）:1944-1950.

［3］ 宋長鳴，向玉林.林業技術效率及其影響因素研究：基于隨機前沿生產函數［J］.林業經濟，2012，34（2）：66-70.SONG Changming， XIANG Yulin.Study of forestry industry’s technical efficiency and its influencing factor—based on stochastic frontier production function model［J］.For Econ,2012,34（2）:66-70.

［4］ 田杰，姚順波.中國林業生產的技術效率測算與分析［J］.中國人口·資源與環境，2013，23（11）：66-72.TIAN Jie,YAO Shunbo.Research on technical efficiency of forestry production in China ［J］.China Popul Resour Environ,2013,23（11）:66-72.

［5］ CARTER D R,CUBBAGE F W.Stochastic frontier estimation and sources of technical efficiency in southern timber harvesting ［J］.For Sci,1995,41（3）:576-593.

［6］ SALEHIRAD N,SOWLATI T.Productivity and efficiency assessment of the wood industry:a review with a focus on Canada ［J］.For Prod J,2006,56（11/12）:25-32.

［7］ 臧良震，支玲，齊新民.天保工程區農戶林業生產技術效率的影響因素：以重慶武隆縣為例［J］.北京林業大學學報（社會科學版）， 2011， 10（4）： 59-64.ZANG Liangzhen,ZHI Ling,QI Xinmin.Influencing factors of farmers forestry production efficiency in natural forest protection project areas:a case of Wulong County in Chongqing City ［J］.J Beijing For Univ Soc Sci,2011,10（4）:59-64.

［8］ 翟秋，李樺，姚順波.后林權改革視角下家庭林地經營效率研究［J］.西北農林科技大學學報（社會科學版），2013， 13（2）： 64-69.ZHAI Qiu,LI Hua,YAO Shunbo.Research on efficiency of household forestland operating based on forest tenure reform:the case of collective forestland in countries of Sha and Shunchuang of Fujian Province ［J］.J Northwest Aamp;F Univ Soc Sci Ed,2013,13（2）:64-69.

［9］ 張春霞，許佳賢，黃森慰，等.基于木材供給生產目標下林農經營規模效率研究：以福建省杉木用材林為例［J］.中南林業科技大學學報（社會科學版）， 2010， 4（2）： 5-7， 31.ZHANG Chunxia,XU Jiaxian,HUANG Senwei,et al.A research into forestry farmer’s management scale efficiency based on the production aim of timber supply:a case study of fir timber in Fujian Province ［J］.J Centr South Univ of For Technol Soc Sci,2010,4（2）:5-7,31.

［10］ MANAGI S.Productivity measures and effects from subsidies and trade:an empirical analysis for Japan’s forestry［J］.Appl Econ,2010,42（30）:3871-3883.

［11］ 徐秀英，付雙雙，李曉格，等.林地細碎化、規模經濟與竹林生產：以浙江龍游縣為例［J］.資源科學，2014， 36（11）： 2379-2385.XU Xiuying,FU Shuangshuang,LI Xiaoge,et al.Forestland fragmentation,economies of scale and bamboo production ［J］.Resour Sci,2014,36（11）:2379-2385.

［12］ 許佳賢，蘇時鵬，黃安勝，等.農戶林業經營效率及其影響因素分析：基于閩浙贛235個固定觀察點6年的調查數據［J］.農村經濟， 2014（11）： 42-46.XU Jiaxian,SU Shipeng,HUANG Ansheng,et al.Farmers forestry management efficiency and its influencing factors:a study based on survey data of 235 fixed observation point 6 years from Fujian,Zhejiang and Jiangxi ［J］.Rural E-con,2014（11）:42-46.

［13］ KORHONEN P J,SIITARI P A.A dimensional decomposition approach to identifying efficient units in large-scale DEA models ［J］.Comput Oper Res,2009,36（1）:234-244.

［14］ TYAGI P,YADAV S P,SINGH S P.Relative performance of academic departments using DEA with sensitivity analysis ［J］.Evalu Progr Plan,2009,32（2）:168-177.

［15］ FRIED H O,LOVELL C A K,SCHMIDT S S,et al.Accounting for environmental effects and statistical noise in data envelopment analysis ［J］.J Prod Anal,2002,17（1）:157-174.

［16］ 李樺，姚順波，劉璨，等.集體林分權條件下不同經營類型商品林生產要素投入及其效率：基于三階段DEA模型及其福建、江西農戶調研數據［J］.林業科學，2014，50（12）：122-130.LI Hua,YAO Shunbo,LIU Can,et al.Input of production factors for different operational types of commercial forests and the technical efficiency in the reform of collective forestry property right system:based on three-stage DEA model and household data of Fujian and Jiangxi Province ［J］.Sci Silv Sin,2014,50（12）:122-130.

［17］ 申津羽，韓笑，侯一蕾，等.基于三階段DEA模型的南方集體林區不同林業經營形式效率研究［J］.南京林業大學學報（自然科學版）， 2015， 39（2）： 104-110.SHEN Jinyu,HAN Xiao,HOU Yilei,et al.Measuring the technical efficiency of different forestry management model in southern community forest area using three-stage DEA analysis ［J］.J Nanjing For Univ Nat Sci Ed,2015,39（2）:104-110.

［18］ 付雙雙，徐秀英，吳偉光.林農毛竹生產效率及其影響因素［J］.浙江農林大學學報，2015，32（4）：596-602.FU Shuangshuang,XU Xiuying,WU Weiguang.Analysis on the determinants of Phyllostachys edulis’s input-output efficiency ［J］.J Zhejiang Aamp;F Univ,2015,32（4）:596-602.

［19］ CHARNES A,COOPER W W,RHODES E.Measuring the efficiency of decision making units ［J］.Europ J Oper Res,1978,2（6）:429-444.

［20］ 羅登躍.三階段DEA模型管理無效率估計注記［J］.統計研究，2012，29（4）：104-107.LUO Dengyue.A note on estimating managerial inefficiency of three-stage DEA model ［J］.Stat Res,2012,29（4）:104-107.

［21］ KUMBHAKAR S C.Decomposition of technical change into input-specific components:a factor augmenting approach［J］.Jpn World Econ,2002,14（3）:243-264.

［22］ 郭軍華，倪明，李幫義.基于三階段DEA模型的農業生產效率研究［J］.數量經濟技術經濟研究，2010（12）：27-38.GUO Junhua,NI Ming,LI Bangyi.Research on agricultural production efficiency based on three-stage DEA model［J］.J Quant Tech Econ,2010（12）:27-38.

［23］ CONWAY M C,AMACHER G S,SULLIVAN J,et al.Decisions nonindustrial forest landowners make:an empirical examination ［J］.J For Econ,2003,9（3）:181-203.

［24］ ZHANG Daowei,OWIREDU E A.Land tenure,market,and the establishment of forest plantations in Ghana ［J］.For Policy Econ,2007,9（6）:602-610.

［25］ 蘇時鵬，馬梅蕓，林群.集體林權制度改革后農戶林業全要素生產率的變動：基于福建農戶的跟蹤調查［J］.林業科學， 2012， 48（6）： 127-135.SU Shipeng,MA Meiyun,LIN Qun.Farmer forestry total factor productivity changes after the collective forestry property rights system reform:based on households surveys in Fujian Province ［J］.Sci Silv Sin,2012,48（6）:127-135.

［26］ 廖文梅，廖冰，金志農.林農經濟林經營效率及其影響因素分析：以贛南原中央蘇區為例［J］.農林經濟管理學報， 2014， 13（5）： 490-498.LIAO Wenmei,LIAO Bing,JIN Zhinong.Efficiency and factors influencing forest farmers’operated economic forest:a case study of Gannan Central Soviet Area ［J］.J Agro-For Econ Manage,2014,13（5）:490-498.

［27］ FOSTER W E,RAUSSER G C.Farmer behavior under risk of failure ［J］.Am J Agric Econ,1991,73（5）:276-288.

Research on bamboo production efficiency based on three-stage DEA model

YANG Shuisheng1,XU Xiuying1,2,FU Jiaoyan1,XU Jinhan1
（1.School of Economic and Management,Zhejiang A ＆ F University,Hangzhou 311300,Zhejiang,China； 2.Rural Development Research Center of Zhejiang,Key Research Base of Philosophy and Social Science of Zhejiang,Hangzhou 311300,Zhejiang,China）

Forestry production efficiency is an important issue in the development forestry economics.Measuring and analysing the efficiency of forestry production scientifically,which reflect truly the actual situation of forestry production efficiency to provide a scientific basis for improving forestry production efficiency.This paper forecasted and evaluated bamboo production efficiency using the three-stage Data Envelopment Analysis（DEA） model.The dataset adopted in the study was first hand data,collected based on 110 surveys from 6 administrative villages in Anji County in Zhejiang Province.The results showed:Among the environment variables,significant negative effects were found for factors such as the years of education for the householder,income ratio for bamboo production compare to total income,the number of family members engaged in the production of bamboo,the total size of land for bamboo production;while the following elements showed significant positive effects for bamboo production:age of household head,involvement in technical training in bamboo production,and certificate of forestland.After excluding the impact of environmental effects and random errors,this studies found that the sample household bamboo production average technical efficiency value showed asignificant decreased,from 0.537 to 0.436,while the average pure technical efficiency value from 0.602 to 0.803,and the average scale efficiency value decreased from 0.891 to 0.556.At the same time,increasing return to scale was witnessed for bamboo production.Based on the above analysis,suggestions,such as reduce the fragmentation of bamboo forestland,guide farmers to allocate bamboo production factors rationally,clarify the property rights of bamboo forestland,strengthen technical training for farmers’ bamboo production,are put forward. ［Ch,6 tab.27 ref.］

forest economics;bamboo production efficiency;three-stage DEA model;environment effect;random error;Anji County

S7-98

A

2095-0756（2017）06-1128-09

10.11833/j.issn.2095-0756.2017.06.022

2016-12-05；

2017-04-05

國家自然科學基金資助項目（71273245）；浙江省財政廳提升地方高校辦學水平專項（農民發展研究創新團隊）

楊水生，從事林業經濟理論與政策研究。E-mail：yangshuisheng@foxmail.com。通信作者：徐秀英，教授，博士，從事林業經濟理論與政策研究。E-mail：zjfcxxy@aliyun.com