張 斌，程長(zhǎng)林，任愛(ài)勝

（1.北京市農(nóng)林科學(xué)院農(nóng)業(yè)信息與經(jīng)濟(jì)研究所，北京 100097；2.北京物資學(xué)院經(jīng)濟(jì)學(xué)院，北京 101149；3.中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)經(jīng)濟(jì)與發(fā)展研究所，北京 100081）

當(dāng)前，我國(guó)畜牧業(yè)生產(chǎn)制約因素多，產(chǎn)業(yè)關(guān)聯(lián)性弱，種養(yǎng)結(jié)合效率低，亟需提高草畜轉(zhuǎn)化效率，優(yōu)化畜牧業(yè)種養(yǎng)結(jié)構(gòu)，提升最終產(chǎn)品附加值[1-2]。草畜高效轉(zhuǎn)化有助于充分利用技術(shù)、人才、資金等要素整合，提高市場(chǎng)要素配置效率，優(yōu)化產(chǎn)業(yè)布局，建立多個(gè)利益相關(guān)主體聯(lián)結(jié)機(jī)制，激活農(nóng)戶、家庭農(nóng)場(chǎng)、合作社等經(jīng)營(yíng)主體，促進(jìn)農(nóng)牧民增收[3-5]。青藏高原地區(qū)人類(lèi)活動(dòng)較少，高海拔、高嚴(yán)寒、降水少等生態(tài)條件給畜牧業(yè)發(fā)展帶來(lái)巨大挑戰(zhàn)。此外，青藏高原地區(qū)傳統(tǒng)畜牧業(yè)發(fā)展依草而牧，但高寒地區(qū)天然草地生產(chǎn)力恢復(fù)緩慢，草地載畜量持續(xù)增加，而人工種草意識(shí)薄弱，草畜矛盾激增，極大威脅青藏高原地區(qū)的天然生態(tài)。因此，青藏高原牧區(qū)草畜轉(zhuǎn)化效率是直接影響高寒地區(qū)畜牧業(yè)發(fā)展效率與質(zhì)量的關(guān)鍵因素。在此背景下，青藏高原社區(qū)畜牧業(yè)選取西藏墨竹工卡與羊八井、青海河南縣與稱多縣、四川紅原與甘孜、甘肅夏河與云南迪慶等8個(gè)社區(qū)示范點(diǎn)，集科研人才、實(shí)用技術(shù)及資金等多項(xiàng)投入，在科研幫扶式模式下帶動(dòng)牧民進(jìn)行天然草地保護(hù)、人工種草、飼草加工、優(yōu)化放牧等活動(dòng)，通過(guò)轉(zhuǎn)變藏區(qū)牧民生產(chǎn)生活方式，探索促進(jìn)青藏高原地區(qū)畜牧業(yè)發(fā)展與牧民增收的產(chǎn)業(yè)發(fā)展模式。本研究旨在對(duì)青藏高原社區(qū)畜牧業(yè)草畜生產(chǎn)轉(zhuǎn)化中各產(chǎn)業(yè)要素的投入產(chǎn)出及影響因素進(jìn)行研究，分析青藏高原社區(qū)畜牧業(yè)發(fā)展效率，為當(dāng)?shù)匦竽翗I(yè)發(fā)展提供借鑒。

1 方法與指標(biāo)選擇

1.1 DEA-Malmquist指數(shù)

1.1.1 數(shù)據(jù)包絡(luò)分析法（DEA） DEA隸屬于多投入多產(chǎn)出的非參數(shù)投入產(chǎn)出效率分析方法體系[6-9]，建立在投入或產(chǎn)出不變的前提下，通過(guò)建立決策單元與有效生產(chǎn)前沿面偏離度結(jié)構(gòu)方程，在相關(guān)線性規(guī)劃約束基礎(chǔ)上，分析決策單元在生產(chǎn)前沿面的投入產(chǎn)出效率[10-12]。

上述方程所得最優(yōu)解λ*、S*-、θ*的經(jīng)濟(jì)學(xué)內(nèi)涵：

①當(dāng)且僅當(dāng)θ*=1，s*+=0，s*-=0 時(shí)，表明決策單元有效且其技術(shù)效率與規(guī)模效率均為有效。

②當(dāng)且僅當(dāng)線性約束方程存在最優(yōu)解且值為1，技術(shù)指標(biāo)值大于0時(shí)，決策單元的規(guī)模效率與技術(shù)效率存在結(jié)構(gòu)性矛盾，DEA效率僅能表達(dá)較弱的有效性。

③當(dāng)θ*＜1時(shí)，決策單元為非DEA有效，表明其生產(chǎn)活動(dòng)技術(shù)效應(yīng)和規(guī)模效應(yīng)均未達(dá)到最佳。

1.1.2 Malmquist指數(shù)及其分解 1953年，Malmquist全要素生產(chǎn)率指數(shù)算法提出后，在Fare、Ross等人的改進(jìn)基礎(chǔ)上[13-15]，形成了從時(shí)期T至T+1的全要素生產(chǎn)率增長(zhǎng)率指數(shù)，其表達(dá)式：

式（1）中，(xt+1,yt+1) 和(xt,yt)分別表示為（t+1）和t 2個(gè)時(shí)期的投入和產(chǎn)出向量；D0t和D0

(2)病毒入侵。在大數(shù)據(jù)時(shí)代，網(wǎng)絡(luò)信息數(shù)據(jù)具有很強(qiáng)的共享性和開(kāi)放性特征，這為病毒入侵提供了良好的條件，并且計(jì)算機(jī)病毒自身大都具有隱藏性、信息儲(chǔ)存性以及可執(zhí)行性等基本特征，它們?nèi)绻行ЫY(jié)合了計(jì)算機(jī)網(wǎng)絡(luò)程度，那么就會(huì)自動(dòng)地釋放出具有破壞性、潛伏性、觸發(fā)性以及感染性等特征的程序或數(shù)據(jù)，進(jìn)而會(huì)影響網(wǎng)絡(luò)數(shù)據(jù)的安全性。在當(dāng)下的病毒入侵當(dāng)中，病毒傳播介質(zhì)主要包括光盤(pán)、硬盤(pán)與軟盤(pán)三種類(lèi)型，在這三種介質(zhì)中運(yùn)行的數(shù)據(jù)可能會(huì)受到病毒入侵，且危害性越強(qiáng)的病毒也會(huì)相應(yīng)地對(duì)網(wǎng)絡(luò)數(shù)據(jù)信息造成更嚴(yán)重的破壞。比如，CIH病毒以及“熊貓燒香”等均是危害計(jì)算機(jī)網(wǎng)絡(luò)運(yùn)行安全性的病毒。

t+1表示在以Tt為參照的t和t+1的距離函數(shù)。

將式（1）變形可知，生產(chǎn)效率可分解為技術(shù)效率與生產(chǎn)前沿2部分：

同理，技術(shù)效率進(jìn)一步分解，可轉(zhuǎn)換為規(guī)模效率（SC）和純技術(shù)效率（PC），即：

換言之，式（3）中技術(shù)效率又可進(jìn)一步分解為技術(shù)效率和規(guī)模效率：

1.2 Tobit模型 考慮到技術(shù)效率具有刪截性質(zhì)，受限因變量的OLS估計(jì)結(jié)果有偏，該情況下采用Tobit模型可以截取掉不在這個(gè)區(qū)間段的因變量值，從而表現(xiàn)出一致性和漸進(jìn)正態(tài)性，得出無(wú)偏的估計(jì)結(jié)果[16]。

式中，TEi*為技術(shù)效率的潛在變量，xi為解釋變量， δ～N(0,σ2)，TEi為實(shí)際技術(shù)效率。

1.3 指標(biāo)選取

1.3.1 投入產(chǎn)出指標(biāo)選取 本文草畜轉(zhuǎn)化是按畜牧業(yè)產(chǎn)業(yè)鏈環(huán)節(jié)從飼草生產(chǎn)、加工、牲畜飼養(yǎng)、育肥、畜產(chǎn)品加工、市場(chǎng)等，實(shí)現(xiàn)產(chǎn)品價(jià)值由草向畜及最終產(chǎn)品的轉(zhuǎn)化。而草畜轉(zhuǎn)化效率屬于技術(shù)效率范疇，是草畜轉(zhuǎn)化過(guò)程中要素投入產(chǎn)出效率。因此，根據(jù)相關(guān)文獻(xiàn)及DEA分析對(duì)效率分析中投入產(chǎn)出指標(biāo)的擇取依據(jù)，考慮指標(biāo)的代表性和相關(guān)性，社區(qū)畜牧業(yè)投入指標(biāo)選取：①科研經(jīng)費(fèi)內(nèi)部支出是項(xiàng)目實(shí)施過(guò)程中配套科研投入經(jīng)費(fèi)；②科技活動(dòng)人數(shù)是項(xiàng)目人力資源的投入，主要是科技人力資源；③R&D經(jīng)費(fèi)內(nèi)部支出額是項(xiàng)目實(shí)施過(guò)程中相關(guān)研發(fā)經(jīng)費(fèi)投入；④R&D人數(shù)是參與重大研發(fā)人力投入；⑤草籽投入額、機(jī)械投入（馬力*時(shí)）是社區(qū)畜牧業(yè)基本生產(chǎn)資料投入。科技產(chǎn)出指標(biāo)選?。孩賹?zhuān)利、發(fā)表論文是主要科技產(chǎn)出；②家畜掉膘、家畜育肥、家畜死亡率是主要畜牧生產(chǎn)產(chǎn)出。采用2012—2016年青藏高原社區(qū)畜牧業(yè)發(fā)展過(guò)程中要素投入和畜牧產(chǎn)出相關(guān)數(shù)據(jù)，通過(guò)DEAP 2.1軟件，得出8個(gè)示范社區(qū)畜牧業(yè)草畜轉(zhuǎn)化效率結(jié)果及分解結(jié)果。

本文投入產(chǎn)出指標(biāo)與相關(guān)文獻(xiàn)中所涉及指標(biāo)體系構(gòu)建差異主要體現(xiàn)在：青藏高原社區(qū)畜牧業(yè)項(xiàng)目直接投入產(chǎn)出指標(biāo)體現(xiàn)在基本生產(chǎn)資料投入、科技投入及科技產(chǎn)出、畜牧業(yè)產(chǎn)出；社區(qū)畜牧業(yè)特征及指標(biāo)選取可量化與代表性等。此外，相關(guān)數(shù)據(jù)來(lái)自于青藏高原社區(qū)畜牧業(yè)各專(zhuān)項(xiàng)基于社區(qū)基線調(diào)查及跟蹤調(diào)查統(tǒng)計(jì)，其統(tǒng)計(jì)特征見(jiàn)表1。

表1 投入產(chǎn)出指標(biāo)統(tǒng)計(jì)描述

1.3.2 Tobit模型變量選取 在模型變量選取上，考慮到因變量效率值的限制，采用受限制的Tobit回歸方法?；谏鐓^(qū)畜牧業(yè)生產(chǎn)活動(dòng)特點(diǎn)，分析人、財(cái)、物、技術(shù)轉(zhuǎn)化對(duì)畜牧業(yè)生產(chǎn)技術(shù)效率的影響，從科研人員投入數(shù)量（RI）、科研經(jīng)費(fèi)投入總額（RFI）、生產(chǎn)資料投入（PI）、牧民培訓(xùn)（ST）4個(gè)方面進(jìn)行Tobit計(jì)量分析。首先，青藏高原社區(qū)畜牧業(yè)發(fā)展離不開(kāi)科研人員和實(shí)用科技成果的支持，是高原草畜業(yè)關(guān)鍵技術(shù)重要載體，反映科研進(jìn)駐對(duì)邊緣高寒地區(qū)發(fā)展畜牧業(yè)的潛力。其次，在科技投入方面，科研人員及資金是實(shí)現(xiàn)科技成果轉(zhuǎn)化的重要要素驅(qū)動(dòng)力，結(jié)合牧民培訓(xùn)情況，能針對(duì)性地反映出青藏高原社區(qū)畜牧業(yè)關(guān)鍵技術(shù)集成推廣成效。

2 結(jié) 果

2.1 示范社區(qū)草畜轉(zhuǎn)化效率評(píng)價(jià)

2.3.1 示范社區(qū)DEA效益整體評(píng)價(jià) 由表2可知，2012—2016年青藏高原社區(qū)畜牧業(yè)草畜轉(zhuǎn)化平均技術(shù)效率為0.863，純技術(shù)效率和規(guī)模效率分別為0.928和0.925。示范社區(qū)中DEA有效社區(qū)共6個(gè)，占75%，分別為夏河社區(qū)、香格里拉社區(qū)、墨竹工卡社區(qū)、河南社區(qū)、紅原社區(qū)以及玉樹(shù)社區(qū)，社區(qū)草畜轉(zhuǎn)化技術(shù)效率基本達(dá)到最優(yōu)，決策單元達(dá)到了技術(shù)效率和規(guī)模效率。羊八井社區(qū)與甘孜社區(qū)為非DEA有效社區(qū)，其中，羊八井社區(qū)DEA非有效是由純技術(shù)效率無(wú)效引起，甘孜社區(qū)則由純規(guī)模效率無(wú)效引起，因此可分別通過(guò)改善技術(shù)環(huán)境、提升技術(shù)創(chuàng)新能力及社區(qū)畜牧業(yè)生產(chǎn)要素投入破解社區(qū)草畜轉(zhuǎn)化技術(shù)效率無(wú)效的困境。同時(shí)，羊八井社區(qū)草畜轉(zhuǎn)化處于規(guī)模報(bào)酬遞增階段，增加要素投入可以提升草畜轉(zhuǎn)化效率水平；而甘孜社區(qū)處于規(guī)模報(bào)酬遞減階段，造成了規(guī)模不經(jīng)濟(jì)，需要合理調(diào)整要素投入。總體而言，由于部分社區(qū)純技術(shù)效率或規(guī)模效率指數(shù)較弱，拉低了8個(gè)社區(qū)草畜轉(zhuǎn)化綜合效率平均水平。

表2 2012—2016年示范社區(qū)草畜轉(zhuǎn)化效率整體結(jié)果

2.1.2 示范社區(qū)技術(shù)效率分解 根據(jù)前沿標(biāo)桿計(jì)算得出2012—2016年草畜轉(zhuǎn)化效率非DEA有效社區(qū)投入冗余量（表3），非DEA有效社區(qū)應(yīng)根據(jù)投入冗余量相應(yīng)減少相關(guān)指標(biāo)投入量。如羊八井社區(qū)應(yīng)相應(yīng)減少科技活動(dòng)人數(shù)、R&D 經(jīng)費(fèi)內(nèi)部支出額等指標(biāo)的投入，進(jìn)而實(shí)現(xiàn)社區(qū)草畜轉(zhuǎn)化效率有效。

2.1.3 示范社區(qū)草畜轉(zhuǎn)化效率動(dòng)態(tài)評(píng)價(jià) 由表4可知，2012—2016年青藏高原社區(qū)畜牧業(yè)草畜轉(zhuǎn)化效率Malmquist指數(shù)為1.095 3，說(shuō)明青藏高原社區(qū)畜牧業(yè)草畜轉(zhuǎn)化效率增長(zhǎng)了9.53%。從年平均變化來(lái)看，2012—2016年，青藏高原社區(qū)畜牧業(yè)草畜轉(zhuǎn)化效率變動(dòng)出現(xiàn)了先上升后穩(wěn)定的特征，2012—2014年為上升期，2014—2016為穩(wěn)定期。

具體而言，2012—2014年，草畜轉(zhuǎn)化效率的Malmquist指數(shù)呈上升趨勢(shì)，其中技術(shù)效率上漲了8.5%，資源配置效率上升了14.7%。從結(jié)合效率區(qū)域差異看，2012—2016年，甘孜社區(qū)草畜轉(zhuǎn)化效率Malmquist指數(shù)為0.91，紅原社區(qū)為0.96，羊八井社區(qū)為0.871，夏河社區(qū)為0.971，墨竹工卡社區(qū)為0.959，河南社區(qū)為0.982，香格里拉社區(qū)為0.984，玉樹(shù)社區(qū)為0.98；各示范社區(qū)草畜轉(zhuǎn)化效率都呈穩(wěn)中有升之勢(shì)，但河南社區(qū)、夏河社區(qū)、紅原社區(qū)等草畜轉(zhuǎn)化效率降幅高于香格里拉社區(qū)及墨竹工卡社區(qū)。就配置效率而言，河南社區(qū)、夏河社區(qū)、紅原社區(qū)呈小幅上升，但總體仍呈下降之勢(shì)，表明要素配置效率下降是影響青藏高原示范社區(qū)草畜轉(zhuǎn)化效率下降的主要原因。墨竹工卡、香格里拉社區(qū)、玉樹(shù)社區(qū)則因規(guī)模效率與純技術(shù)效率下降導(dǎo)致技術(shù)效率下降。

表3 非DEA有效社區(qū)的投入冗余量

表4 示范區(qū)草畜轉(zhuǎn)化效率Malmquist指數(shù)及其分解

2.2 草畜轉(zhuǎn)化技術(shù)效率影響分析 由表6可知，科研人員投入量對(duì)青藏高原社區(qū)草畜轉(zhuǎn)化效率有積極的影響，且顯著性水平較高。科研人員對(duì)草畜轉(zhuǎn)化效率的作用主要體現(xiàn)在對(duì)社區(qū)牧民的陪伴式影響，是支撐高原畜牧業(yè)發(fā)展的核心之一，有助于增強(qiáng)青藏高原地區(qū)畜牧業(yè)發(fā)展軟實(shí)力。其次，科研資金投入與草畜轉(zhuǎn)化效率呈正相關(guān)關(guān)系，表明科研資金投入不僅能促進(jìn)社區(qū)畜牧業(yè)科技創(chuàng)新、技術(shù)成果轉(zhuǎn)化，同時(shí)科研資金投入對(duì)促進(jìn)青藏高原地區(qū)畜牧業(yè)生產(chǎn)中種養(yǎng)環(huán)節(jié)要素共享，提高產(chǎn)業(yè)關(guān)聯(lián)性具有積極影響。第三，生產(chǎn)資料投入與草畜轉(zhuǎn)化效率呈正向關(guān)聯(lián)，但影響力最小。生產(chǎn)資料投入是體現(xiàn)科技價(jià)值的重要途徑，青藏高原各示范社區(qū)在選址時(shí)充分考慮當(dāng)?shù)匦竽翗I(yè)生產(chǎn)資料供給狀況，因此在社區(qū)畜牧業(yè)發(fā)展過(guò)程中，盡管生產(chǎn)資料投入對(duì)社區(qū)畜牧業(yè)草畜轉(zhuǎn)化能起到積極作用，但相比其他科技投入要素而言，生產(chǎn)資料對(duì)草畜轉(zhuǎn)化效率的促進(jìn)作用較小。最后，牧民培訓(xùn)對(duì)促進(jìn)青藏高原社區(qū)畜牧業(yè)草畜轉(zhuǎn)化具有積極影響。社區(qū)牧民科技培訓(xùn)是高原牧區(qū)科技活動(dòng)的重要構(gòu)成，種植、放牧、養(yǎng)殖等技術(shù)培訓(xùn)可以加強(qiáng)牧民自我生產(chǎn)能力、培養(yǎng)獨(dú)立經(jīng)營(yíng)的牧民帶頭人、提高草畜轉(zhuǎn)化效率，從而使青藏高原社區(qū)畜牧業(yè)經(jīng)營(yíng)由幫扶式向獨(dú)立式轉(zhuǎn)型。

表5 2012—2016年示范社區(qū)草畜轉(zhuǎn)化效率Malmquist指數(shù)及分解結(jié)果

表6 Tobit模型回歸結(jié)果及統(tǒng)計(jì)量

3 結(jié) 論

本文首先運(yùn)用DEA-Malmquist分析法對(duì)2012—2016年青藏高原示范社區(qū)草畜轉(zhuǎn)化效率進(jìn)行了分析，并通過(guò)Tobit模型分析影響草畜轉(zhuǎn)化效率的關(guān)鍵因素，得出如下結(jié)論：

第一，青藏高原示范社區(qū)草畜轉(zhuǎn)化效率測(cè)算中，DEA有效社區(qū)6個(gè)，分別為紅原社區(qū)、香格里拉社區(qū)、墨竹工卡社區(qū)、河南社區(qū)、夏河社區(qū)、玉樹(shù)社區(qū)，2012—2016年技術(shù)效率均值為0.878。此外，DEA非有效社區(qū)主要由純技術(shù)效率變化引起，可通過(guò)提高技術(shù)創(chuàng)新力、加大技術(shù)推廣與應(yīng)用等提高示范社區(qū)草畜轉(zhuǎn)化效率。

第二，草畜轉(zhuǎn)化效率Malmquist指數(shù)為1.095 3，說(shuō)明草畜轉(zhuǎn)化效率在2012—2016年上升了9.53%，其中Malmquist指數(shù)穩(wěn)步增長(zhǎng)，說(shuō)明青藏高原社區(qū)畜牧業(yè)示范社區(qū)草畜轉(zhuǎn)化效率整體有效。在研究期間，Malmquist指數(shù)呈現(xiàn)先上升后趨穩(wěn)的態(tài)勢(shì)，邊際規(guī)模報(bào)酬先增后減，符合區(qū)域經(jīng)濟(jì)增長(zhǎng)階段規(guī)律。Malmquist指數(shù)分解可知，青藏高原社區(qū)畜牧業(yè)示范社區(qū)草畜轉(zhuǎn)化效率下降是由資源配置效率下降引起，因此，應(yīng)提高資源運(yùn)轉(zhuǎn)效率，加強(qiáng)項(xiàng)目管理，增強(qiáng)實(shí)用技術(shù)應(yīng)用效果。

第三，從示范社區(qū)結(jié)合效率來(lái)看，香格里拉社區(qū)草畜轉(zhuǎn)化效率最高，其次是玉樹(shù)社區(qū)、河南社區(qū)、夏河社區(qū)、紅原社區(qū)，再次是羊八井社區(qū)和甘孜社區(qū)。從Malmquist指數(shù)分解來(lái)看，示范社區(qū)資源配置效率均有下降，要素配置效率下降是影響示范社區(qū)草畜轉(zhuǎn)化效率下降的主因。墨竹工卡、香格里拉社區(qū)、玉樹(shù)社區(qū)技術(shù)效率下降，是由純技術(shù)效率和規(guī)模效率下降共同引起。

第四，從草畜轉(zhuǎn)化效率影響因素看，科研人員投入、科研資金投入、生產(chǎn)要素投入以及牧民培訓(xùn)對(duì)青藏高原社區(qū)畜牧業(yè)草畜轉(zhuǎn)化具有正向影響。從其影響因素優(yōu)先順序來(lái)看，科研人員投入影響最大，生產(chǎn)要素投入影響最小，反映科研人力資源投入對(duì)青藏高原社區(qū)、牧民、牧業(yè)的發(fā)展有重大作用。

基于此，在青藏高原畜牧業(yè)發(fā)展過(guò)程中，政府應(yīng)借鑒社區(qū)畜牧業(yè)發(fā)展模式，加大科研投入，促進(jìn)科研幫扶，有效解決青藏高原地區(qū)生態(tài)約束下的牧民增收與牧業(yè)高效發(fā)展問(wèn)題。