農作物病蟲害專業化防治服務對農藥施用強度的影響

應瑞瑤 徐斌

摘要農藥過量施用已造成嚴重的負外部性。推行農作物病蟲害專業化統防統治服務可以實現安全、科學、合理用藥，有效緩解農業面源污染問題和提升食品的質量安全水平。為了探究病蟲害專業化統防統治服務在緩解農業面源污染方面的效果，文章利用全國七省水稻病蟲害專業化防治服務調查數據，基于傾向得分匹配方法（PSM），剔除樣本自選擇內生性問題的影響，分析農戶將水稻病蟲害防治環節外包給植保專業化服務組織（“統防統治”）與“自防自治”戶相比使用農藥的種類和施用次數的差異，即植保專業化防治是否取得了較好的環境效應。研究結果表明，植保專業化服務顯著減少了農藥施用強度，提高了無公害低毒農藥的采用比例；并且其效果在小農戶和規模種植大戶之間存在明顯差異，與采納病蟲害統防統治服務的規模種植大戶相比，小規模種植戶在采納病蟲害統防統治服務后在降低農藥施用強度，提高無公害低毒農藥的應用比例方面效果更顯著。基于此，文章建議強化對植保專業化防治項目的財政扶持力度，提高病蟲害專業化防治覆蓋率；重點鼓勵小規模種植戶采納病蟲害統防統治服務；鼓勵農戶土地流轉，促進農業生產規模化經營，加強對種植專業戶的病蟲害防治方面的培訓和引導，并著力發揮專業大戶的示范帶頭作用，提高統防統治的效率，從而促進食品質量安全水平的提升和農業生態環境的保護。

關鍵詞統防統治；傾向得分匹配；農藥施用；內生性

中圖分類號S9

文獻標識碼A文章編號1002-2104（2017）08-0090-08DOI：10.12062/cpre.20170322

隨著我國農村勞動力的轉出，農業部門從業者老齡化的現象非常普遍，難以勝任技術要求高、勞動強度大的農作物病蟲害科學防治[1]。勞動者的病蟲害防治技術信息知識顯著影響農戶的農藥施用量[2]。一家一戶分散防治方式，由于大部分農戶無法掌握專業的植保知識，習慣于已有的經驗或農資商店的推薦，缺乏科學防治的依據[3]，難以及時、有效、合理使用農藥，導致農藥使用過量，造成農業面源污染，還容易造成農藥殘留超標所引發的食品安全問題。農作物病蟲害“統防統治”是農業社會化服務的主要形式之一。所謂統防統治即把一定地域內農作物病蟲害防治環節統一交給專業化服務組織，服務組織根據病蟲害疫情預警，統一防治時間、統一防治藥劑、統一組織實施進行防治。統防統治使得植保部門能夠對數量有限的專業化服務組織提供技術支持和常態化的人員技術培訓，使專業化防治人員的技術水平不斷得以提升，可以實現安全、科學、合理用藥，保障農產品質量安全，保護農業生態環境。本文利用全國七省水稻病蟲害專業化防治服務調查數據，通過計量模型研究統防統治是否能夠減少農藥施用強度，從而帶來環境保護效應；在研究方法方面應用傾向得分匹配方法（Propensity Score Matching，PSM），剔除了農戶的自選擇行為對估計結果產生影響，以避免評價的誤差。

1文獻綜述

有關農戶自選擇行為對農戶環境友好型生產行為或技術采納影響的實證文獻有兩種情況：一是忽略農戶的自選擇行為，認為其是外生的，二是將農戶自選擇行為作為內生變量來處理。大部分研究屬于第一種情況，即沒有解決農戶自選擇產生的內生性問題。如Nadine利用方差分析對美國蘋果種植戶的研究發現，由于規模種植戶更多的參與的政府的培訓項目，因而更傾向支持政府的禁用劇毒農藥計劃[4]。Yila等采用逐步線性回歸模型發現采納農業生產社會化服務對農戶的土壤保護性耕作有產生顯著負向影響[5]。Rahman利用Tobit模型研究決定農戶施藥的影響因素，發現農戶是否簽訂訂單會顯著改變農戶的施藥行為[6]。Tamini Lota D指出社會群體的網絡關系對農戶的清潔技術采納有顯著影響，但無法確定這種影響是通過群體的信息渠道還是影響農戶行為的主觀規范來產生作用[7]。還有作者如Hamilton J 等采用案例分析方法論述了北美山區農戶加入合作組織顯著減少了農藥的施用量[8]。國內也有一些學者探討有關農戶選擇變量對農戶環境友好型技術采納的影響，如張云華、馬九杰、孔祥智等以山西、陜西和山東15縣（市） 353個蘋果種植戶探討了與涉農農戶簽訂合同和加入農民專業技術協會對果農采納用無公害和綠色農藥行為的影響，發現農戶與涉農農戶和農業專業技術協會的聯系是影響農戶采用無公害及綠色農藥乃至其它綠色農業技術的主要因素[9]。褚彩虹、馮淑怡等利用聯立雙變量 Probit 模型研究了農戶是否為合作社成員和農業技術培訓經歷對農戶采用有機肥的影響，認為農戶參加合作社對農戶采用有機肥有顯著正向影響，而農民的農業技術培訓對農戶是否采用有機肥的影響不顯著[10]。

由于上述研究沒有考慮有關農戶選擇變量內生性問題，其潛在假設為農戶參加某個項目或活動的行為是一個外生或隨機選擇的過程。但事實上，農戶參與某個項目或活動與農戶的環境友好行為或技術采納往往受到很多共同因素的作用，例如農戶的環境認知、農戶的資源稟賦及社會資本等，這些錯綜復雜而又經常不可觀測的因素往往使農戶的環境友好型生產行為與農戶參加合作社、接受農技推廣活動變量由于自選擇效應而產生模型的內生性問題，比如參加合作社農戶本身的環境保護認知就比其它農戶強烈，因而參加合作社動機也就更強烈，農戶的環境友好行為是由于其自身比較強的環保認知而非參加合作社的影響[11]。為糾正內生性產生的計量偏誤，一些文獻利用了工具變量或面板數據倍差等方法，得到的研究結論更加準確。Feder G 等利用 1991—1999 年印度尼西亞的310個農戶調查面板數據，運用倍差法考察了農民培訓學校對入學的農戶及其周邊農戶采納IPM效果的影響，發現農民培訓學校對于其畢業農戶及其周邊農戶采納IPM的影響不顯著；通過擴展的三階段模型來降低樣本選擇偏誤，并發現在肯尼亞選擇出口的菜農的農藥施用量明顯低于以國內為市場菜農的農藥施用量[12]；Ali A 等采用配對方法（matching approach）估算了巴基斯坦農民培訓對于棉農采納IPM（病蟲害綜合防治）的影響效果，發現參加農民培訓學校的農民的平均IPM采納效果明顯好于沒有參加培訓的農民[13]。

從以上文獻梳理可知，從方法論的視角來看大部分的相關研究并未解決農戶自選擇變量的內生性問題，因而其結論的準確性因內生性的影響而備受質疑，同時筆者國內尚未檢索到有關統防統治對農作物病蟲害防治農藥施用強度和品種影響的實證分析。本文將在以上兩方面對現有研究進行延伸，利用傾向得分匹配方法，在解決自選擇偏誤問題的基礎上探索關統防統治對農作物病蟲害防治的影響，并嘗試揭示該影響的作用機制。

2研究方法與數據

2.1模型設定

2.1.1農戶統防統治服務選擇行為模型

為考察農戶選擇采納統防統治服務的影響機制，為后面的匹配分析做準備，我們首先對農戶病蟲害防治方式選擇行為模型做如下設定：農戶對病蟲害統防統治服務的選擇行為是二元離散選擇變量，即Gi=1表示農戶采納病蟲害統防統治服務；Gi=0，表示農戶自行防治病蟲害。為分析各變量對農戶是否采納統防統治服務的影響，本文建立Probit 模型如下：

P（Gi=1）=（β0+β′iZi+εi） （1）

其中，Zi代表影響農戶病蟲害統防統治服務選擇行為的外生變量向量，借鑒已有研究成果[14]，Zi主要包括農戶個人特征、農戶生產經營特征和非農就業狀況等，其中εi為隨機擾動項。

2.1.2農戶傾向得分匹配分析

Behrman J R指出，隨機分配實驗樣本可以消除樣本可觀察和不可觀察個人特征的差異所帶來的估計偏誤，從而測算出無偏的ATT值[15]。在隨機分配實驗樣本的情況下，Y1i，Y0i|Di，并且E[Y0i|Di=1]=E[Y0i|Di=0]，則ATT值可以表示為：

E[δ|Di=1]=E[Y1i|Di=1]-E[Y0i|Di=1]

=E[Y1i|Di=1]-E[Y0i|Di=0] （2）

但農戶采納病蟲害統防統治或選擇自防自治并不是由隨機分配確定，估計系數存在樣本自選擇所帶來的內生性偏差，影響農戶是否采納統防統治的個人特征差異會影響農戶農藥使用強度、使用農藥種類以及病蟲害防治費用。觀察值是非隨機實驗的，式（2）不成立，應用OLS得到的估計系數存在偏誤，因此只有在實驗結果與實驗分配方法無關的條件下，估計系數才是無偏的，但是觀察值E[Y0i|Di=1]是不存在的，因此解決這一問題常用一種方法就是傾向得分配對法（Propensity Score Matching，PSM），PSM首先假定X⊥T|P（X），是指如果兩個農戶采納統防統治概率如果相同，那么影響其采納統防統治的變量是平穩的，即影響農戶是否采納統防統治的因素在兩個農戶之間不存在顯著性異，可通過變量的平行性檢驗來確定；其次假定（Y0，Y1）⊥T|X，也就是說如果保持農戶采納統防統治的影響因素不變，因變量潛在結果與農戶是否采納病蟲害統防統治服務無關。在農戶是否采納統防統治服務可以由一組可觀察的變量完全解釋的情況下，將多維度的影響因素歸結為一維采納概率得分，從而把采納統防統治服務的農戶和與其特征類似的未采納統防統治服務的農戶在多個維度進行匹配，并假定按照這種匹配原則，因變量在采納統防統治服務組與非采納統防統治服務是隨機分配的。進而測算出統防統治對農作物病蟲害防治行為的凈影響。分析流程如下，首先對農戶是否采納統防統治服務的選擇行為方程進行回歸：

PS（ν）=Pr[ser_adop=1|ν]=E[ser_adop|ν] （3）

其中，ν表示影響農戶采納統防統治服務的變量，PS表示農戶采納統防統治服務的概率，即傾向分數（Propensity Score）。由式（3），可以擬合出每一個樣本農戶采納統防統治服務的傾向分數，用以在樣本農戶間進行匹配。

接著，根據Becker et al[16]提出的平均處理效應計算公式，求得病蟲害統防統治服務對參與農戶病蟲害防治的平均處理效應（average effect of treatment on the treated，ATT）：

ATT=E[G1i-G0i|Pi=1]

=E{E[G1i-G0i]|Pi=1，PS（νi）}

=E{E[G1i|Pi=1，PS（νi）]}-

E{E[G0i|Pi=0，PS（νi）|Pi=1]} （4）

其中，G1i、G0i分別表示農戶采納統防統治的病蟲害防治特征（施藥次數、無公害農藥應用比例）變量與未采納統防統治的農戶病蟲害防治特征變量。根據式（3）擬合得到的PS（νi）數據是連續的，連續變量的精確匹配是困難的，但可以通過其它非精確匹配方法進行匹配。常用的方法有：最近鄰匹配（knearest neighbors matching）、半徑匹配（radius matching）與核匹配（kernel matching）。最近鄰匹配的方法可以表示為：

c（i）={‖PSi（ν）-PSj（ν）‖} （5）

其中，下標i為采納統防統治服務的農戶，下標j為沒有采納統防統治服務的農戶，c（i）表示跟農戶i匹配成功的所有農戶j，他們是傾向分數與農戶j最為接近的一些農戶。

同理，半徑匹配的方法可以表示為：

c（i）={‖PSi（ν）-PSj（ν）‖

其中，r表示為處理組農戶在控制組進行匹配時，所設定的匹配得分距離，在此距離范圍之內的控制組農戶與處理組農戶匹配成功。然后，根據Becker et al[16]給出的計算ATT的公式有：

ATTM=1/NT∑i∈TGTi-1/NT∑i∈CωjGCi （7）

其中，M表示我們選擇的匹配原則，例如半徑匹配或核匹配原則等，T為采納統防統治服務的農戶（處理組），C為沒有采納統防統治服務的農戶（控制組）。NiC表示跟處理組農戶i成功匹配但是沒有采納統防統治的控制組中農戶的數量，那么式（7）中的權重可寫成ωj=∑iωi j，其中 ωi j=1/NCi。根據權重ωj不變且采納統防統治服務的農戶相對獨立，則可以通過以下公式估計ATT的方差。

Var（ATTM）=1/NTVar（GTi）+（1/NT）2∑j∈C（ωj）2Var（GCj） （8）

與最近鄰匹配和半徑匹配方法顯著不同，核匹配方法是一種基于非參數的估計。為了在控制組中尋找能跟農戶i匹配的對象，則把農戶i的PS值周圍沒有采納統防統治服務的農戶進行加權計算，根據農戶i和農戶j的PS值之差確定其權重。所以，用來跟農戶i進行匹配的農戶并不存在，而是通過相關因素創造出的虛擬農戶（hypothesized farmer），其所對應G的值就是反事實農戶病蟲害防治特征變量，也就是對各農戶j的病蟲害防治特征變量求加權平均。核匹配法ATT計算公式可以表示為：

ATTK=1/NT∑i∈TGTi-∑i∈CGCiB（（PSj-PSi）/hn）∑i∈CB（（PSk-PSi）/hn） （9）

其中，控制組中用來進行匹配的農戶j的數目由帶寬參數（bandwidth）hn 決定；匹配農戶病蟲害防治變量GCj的權重由高斯核函數（Gaussian kernel function） B來決定。

2.2數據來源及樣本概況

本文所用的數據來自2012年7月至2013年3月對我國水稻主產區病蟲害統防統治服務的調查。為了體現自然氣候和經濟發展水平不同地區推行統防統治的差異，使樣本具有較好的代表性，本文選擇經濟發展水平和大田作物種植狀況差異明顯的七個省份作為調研地區，調查地點具體包括位于東部沿海、經濟發達的糧食主產區的山東省、江蘇省；位于中部地區，經濟相對欠發達的河南省、湖北省、江西省和地處北方的水稻主產區的黑龍江省；西部地區選擇水稻大省四川省。以上樣本地區都是我國的水稻主產區，大田農作物種植面積廣泛，較適合開展病蟲害統防統治。調查人員共選取上述項目示范內51個村作為抽樣點。按分層抽樣方法在每個縣選取了1—2個村，每個村再隨機選擇20戶左右的農戶進行入戶調查。共獲得問卷1 151份，最終獲得有效問卷1 059份，問卷有效率92%。農戶調查問卷涵蓋農戶個體特征及家庭人口就業特征的基本情況、農戶生產特征和水稻生產成本收益情況3個方面。涵蓋信息廣泛，具有較高的真實性和代表性。

本文通過調研主要了解樣本留守務農和農戶外出務工人員的自身特征、農戶家庭特征、病蟲害防治的方式情況和效果、統防統治的推行情況和農戶對統防統治以及農業生產面源污染等問題。其中，農戶特征包括農戶的年齡、文化程度、健康狀況等；農戶家庭特征包括外出打工的類型、家庭收入構成、種植規模等；農業生產經營特征主要包括病蟲害防治社會化服務情況、統防統治服務價格和農技推廣情況等題項。具體變量及描述統計見表1。

3實證分析

3.1被解釋變量說明

本文所講的農藥施用強度以一季中晚稻種植過程中的施藥次數來衡量；農藥施用品種以一個樣本施用的符合無公害標準的農藥種類數占病蟲害防治過程中施用的農藥種類總數的比例來表示。從整體樣本來看，在2012年參加統防統治水稻施藥次數的均值為5.69，自防自治農戶水稻施藥次數為7.86，統防統治水稻施藥次數的均值比自防自治戶施藥次數的均值減少了27%。統防統治水稻施藥次數的中位數為6次，比自防自治戶水稻施藥次數的中位數8次少2次；在2012年統防統治水稻施用高效低毒無公害農藥的比例的均值為0.98，自防自治農戶水稻施用高效低毒無公害農藥的比例的均值為0.56，統防統治水稻施用高效低毒無公害農藥的比例的均值比自防自治戶水稻施用高效低毒無公害農藥的比例均值多達到75%，統防統治水稻無公害農藥應用比的中位數為0.99，比自防自治戶水稻無公害農藥應用比的中位數0.55次高80個百分點。具體數據見表2。

3.2參與統防統治的內生性

正如前文所述，很多不可觀測的因素（如農戶資源稟賦，農戶與農技部門的聯系等）可能同時影響農戶的病蟲害科學防治行為和農戶參與統防統治的傾向。例如，參加統防統治的農戶，可能本身的病蟲害防治水平就較高，而這種“錦上添花”式的選擇過程干擾了統防統治與農戶病蟲害規范防治行為之間的因果關系的評估，因為即使沒有采用統防統治，這些農戶在病蟲害規范防治方面的表現比其他農戶更好。所以，如果沒有充分考慮這一由農戶自選擇導致的內生性問題所產生的估計系數偏誤，將會對統防統治與農戶自防自治在農藥施用強度和施用品種方面的差異評價產生干擾，甚至得出錯誤的政策效果評價結論。

因此，我們首先通過表3在直觀上對樣本中采納病蟲害統防統治的農戶與自防自治農戶就各自影響其采納統防統治服務的變量特征進行比較，對比結果如表3所示：從農戶非農務工收入水平和打工地點看，我們發現處理組的大部分變量均值的數值均大于控制組變量均值，而且TT 均值檢驗結果表明，農戶年齡、受教育水平、外出打工地點和非農就業收入在統計學意義上存在高度顯著的組間差別。這說明采納病蟲害統防統治的農戶其非農就業特征和社會人口統計特征與自防自治農戶有明顯差別，從另一個角度證明可能存在農戶自選擇導致的內生性問題。

從上文分析可知，農戶是否采納統防統治與個人特征、社會人口特征、非農就業等變量有關，這些變量包括年齡、文化程度、農戶種植規模和種植品種等。另外，農戶是否參與統防統治還和農戶與農技部分之間的聯系有關：與農技部門聯系密切的農戶往往有較好的病蟲害防治信息渠道，因此其參加統防統治的可能性也越大。在我們的問卷中，設計了多個問題，包括每年與農技員聯系次數、與農技員的聯系方式、與農技員聯系難易程度和農技員的指導是否及時等作為衡量該農戶與政府農技部門聯系緊密程度的指標，并將其加入統防統治參與的選擇方程。選擇方程的估計結果由于篇幅所限未列出。本文根據統防統治參與方程估計中的確定的顯著影響因素作為傾向分數匹配的共同支持檢驗變量（common support），并運用PSM方法解決農戶采納統防統治服務的內生性問題。

3.3傾向得分匹配分析

根據傾向得分匹配分析的原理，我們基于上文中統防統治參與選擇模型擬合每一個樣本農戶采納統防統治服務的傾向分數，并依據該分數對采納統防統治的農戶和沒有采納統防統治的農戶進行匹配。筆者選擇文獻中經常使用的最近鄰匹配方法，并將半徑匹配方法和核密度匹配方法得到的ATT值作為檢驗穩健性的手段。圖1為匹配前的傾向得分分布圖，縱軸代表概率密度，橫軸代表傾向得分值。通過非參數K密度方法（kdensity）對處理組與控制組的傾向分數密度分布函數進行近似，虛線為控制組（自防自治戶）情況，實線為處理組（采納統防統治戶）情況。通過實線與虛線分布圖可以發現，自防自治組農戶傾向分數的頻率峰值在0.1左右，采納統防統治組農戶的傾向分數的頻率峰值在0.7左右。總體來看，采納統防統治組農戶的傾向得分顯著高于未采納統防統治服務組。

圖2是經過最近鄰匹配后的傾向分數分布圖，因為匹配后的控制組中只保留了與相應參與農戶傾向得分最為相近的樣本，而剔除了匹配失敗的樣本農戶，所以其傾向分數分布發生了明顯改變：表示控制組傾向分數概率分布的虛線向右移動明顯，且最高頻率值移動到0.5左右；處理組傾向分數概率分布的實線則與匹配前相比沒有發生改變。整體上看，兩組傾向分數分布特征基本近似，很直觀地表明：按照傾向得分最近鄰方法匹配后有效修正了兩組間傾向分數偏差，匹配效果理想。

經過最近鄰匹配后，依據前文介紹的方法計算農戶采納統防統治服務的平均效應值。除總體樣本外，我們還將對規模種植的農戶進行單獨分析，求得其ATT值。由于參與統防統治的回歸結果顯示，規模種植的農戶更可能參加統防統治，因此農戶是否采納統防統治的內生性問題在規模種植農戶中表現尤為明顯，對其單獨進行分析可以更加精確地測算統防統治服務對于農藥施用強度和施用品種的影響。同時，專業種植大戶也是農業部門在推行病蟲害統防統治時重點支持的對象，所以本文對其重點討論，力求通過對專業種植大戶采納統防統治服務的ATT值的分析來測算統防統治的實施效果。表4總結了對子樣本計算的ATT值。表4中，分別計算了匹配前、后處理組與控制組的平均農戶采納病蟲害統防統治的環境效應（用病蟲害防治次數或高效低毒無公害農藥的比例表示）以及它們的組間差異，匹配后的組間差（ATT值）就是本文所重點討論的實施統防統治的環境效應。我們首先分析用全體樣本估算的病蟲害防治次數的變化，其ATT值為-2.297，且顯著性水平為1%，說明采納統防統治使水稻種植過程中病蟲害防治次數減少2.3次。同時我們還可以看到，在沒匹配以前處理組與控制組的差別是 -2.177，表明如果忽視了農戶采納統防統治服務時所產生的內生性問題，統防統治的作用效果將會被低估。當因變量為高效低毒無公害農藥的比例時，ATT值為0.415，且在1%的水平上顯著，說明病蟲害統防統治使高效低毒無公害農藥的應用比例上升了41.5百分點。在沒匹配以前處理組與控制組的差別是0.01，可以看出如果忽視采納統防統治服務的內生性問題同樣會高估實施統防統治的效果。以上對比分析說明，病蟲害統防統治的推行確實能夠提升病蟲害防治的環境保護和食品質量安全水平，并且在促進農戶應用高效低毒無公害農藥的比例方面更為顯著。

接下來，我們重點關注規模種植（種植面積大于等于9畝）農戶。我們從統防統治戶樣本總體中提取規模種植的農戶作為子樣本，估算推行統防統治水稻病蟲害防治的影響采用的方法同上。表5的結果顯示，無論是農戶病蟲害防治次數還是高效低毒無公害農藥的比例，在忽視采納統防統治內生性的情況下，都會導致對病蟲害統防統治的實際效果的估計產生偏誤（匹配后農戶病蟲害防治次數ATT值提高了5個百分點；高效低毒無公害農藥的比例的 ATT 值降低了2.3個百分點）。并且，通過與總體樣本統防統治對施藥次數和品種影響的比較分析，可以看出子樣本匹配后的ATT值和顯著性均有不同程度下降：當因變量為農戶病蟲害防治次數時，與全樣本的ATT值-2.297相比，規模種植農戶組的ATT值是-1.206，顯著下降；當因變量為高效低毒無公害農藥的比例時，與全部樣本農戶的ATT值0.415相比，規模種植農戶組的ATT值降低到0.297。可以看出，采納統防統治服務對水稻病蟲害防治環境效應的提升效果在那些非規模種植農戶中表現得更顯著。原因可能是：首先規模種植農戶本身的資源稟賦較高，因而本身的病蟲害防治水平就高，其次，規模種植的農戶更關注生產成本，其參加統防統治前的施藥次數就不多。因此，相對于小農戶，其參加統防統治減少農藥應用次數的潛力小，這也正印證了在沒有統防統治的情況下，小農戶比大農戶傾向于應用更多的農藥；相對于小規模種植戶，即使沒有參加統防統治大戶也能更好多地應用無公害農藥，可能的原因是：一方面，大戶的經營規模比較大，比較利于政府監管部門控制其農藥來源渠道；另一方面，很多農技推廣都是以大戶為依托進行，其會受到更多的病蟲害防治方面的培訓，因而比小規模種植戶傾向于應用高效低毒的無公害農藥。

3.4穩健性檢驗

為驗證前文評價統防統治對農藥施用強度和施用品種影響的準確性，本文嘗試通過多種匹配方法檢驗研究結論的穩健性。筆者具體做法是對不同樣本情形的觀察值進行核匹配與半徑匹配，然后分別計算匹配后樣本的ATT值，將其分別與前文各估算結果比對，發現用本文所得到的研究結論不因匹配方法的改變而發生變化：統防統治的推行顯著改善了水稻病蟲害防治造成的農藥殘留超標和農業生態環境污染，尤其是在應用高效低毒無公害農藥的比例方面；同時，相對于統防統治總體樣本的效果，小規模種植農戶樣本的ATT值更高，說明采納統防統治對小規模種植的農戶作用效果更加顯著。

4結論與討論

統防統治的推行對于提升病蟲害防治的環境效應是明顯的。一方面，在統一施藥時間、統一防治技術和統一施用農藥種類的前提下，可以最大限度的發揮病蟲害防治的規模效應，病蟲害統防統治是農業生產規模經營在這一環節的具體體現。要準確評估統防統治對于減少過量施用農藥和濫用高毒農藥的作用，就必須剔除內生性對評價結果的干擾，因此，本文采用得分傾向匹配法（PSM）排除內生性的影響。分析結果表明：①統防統治在確保產量不減的情況下，減少了病蟲害防治的施藥次數，這也印證了小農戶在自防自治時，存在過量施用農藥現象。同時統防統治也促進了農戶應用無公害農藥的傾向。②與全體樣本相比，采納統防統治對于規模種植戶減少施藥次數效果不如全體樣本的效果明顯，這主要是規模種植戶自身的防治技術水平較高的原因，另一方面，統防統治顯著促進了規模種植戶無公害農藥的應用。

所以，無論在施用農藥次數，還是應用農藥的種類方面，統防統治都顯著促進了病蟲害規范防治。因此要加大對統防統治項目的財政支持力度，準確地防控病蟲害，提高防治效果，不斷提升病蟲害統防統治覆蓋率；重點鼓勵小規模種植戶采納病蟲害統防統治服務；促進農戶土地流轉，推動農業生產規模經營。加強對種植專業戶的病蟲害防治方面的培訓和引導，并著力發揮專業大戶的示范帶頭作用，提高統防統治的效率，同時還要依托專業種植大戶和農機專業合作社發展植保合作社等統防統治服務組織，促進食品質量安全水平的提升和農業生態環境的保護。

（編輯：劉呈慶）

參考文獻（References）

[1]徐偉.提升植保統防統治服務能力[J].江蘇農村經濟，2011（11）：31-32. [XU Wei. Improve regional pest control service abilities [J].Jiangsu rural economy，2011（11）：31-32.]

[2]黃季焜， 齊亮， 陳瑞劍. 技術信息知識、風險偏好與農民施用農藥[J]. 管理世界， 2008（5）：71-76. [HUANG Jikun， QI Liang， CHEN Ruijian. The knowledge about technology information， predilection for running risks and peasants application of pesticides[J]. Management world， 2008（5）：71-76.]

[3]王志剛， 呂冰. 蔬菜出口產地的農藥使用行為及其對農民健康的影響——來自山東省萊陽、萊州和安丘三市的調研證據[J]. 中國軟科學， 2009（11）：72-80. [WANG Zhigang， LV Bing. Farmers pesticide use and its impacts on farmers health in exportoriented vegetable region： evidences from an investigation in Laiyang， Laizhou and Anqiu Cities， Shandong Province [J]. China soft science， 2009（11）：72-80.]

[4]NADINE L F B J. Azinphosmethyl （AZM） phaseout： actions and attitudes of apple growers in Washington State[J]. Renewable agriculture and food systems， 2011， 26（4）：276-286.

[5]YILA O M， THAPA G B. Adoption of agricultural land management technologies by smallholder farmers in the Jos Plateau， Nigeria[J]. International journal of agricultural sustainability， 2008，6（4）：277-288.

[6]RAHMAN S. Farmlevel pesticide use in Bangladesh： determinants and awareness[J]. Agriculture， ecosystems & environment， 2003，95（1）：241-252.

[7]TAMINI L D A nonparametric analysis of the impact of agrienvironmental advisory activities on best management practice adoption： a case study of Quebec[J]. Ecological economics， 2011，70（7）： 1363-1374.

[8]HAMILTON J， SIDEBOTTOM J. Mountain pesticide education and safety outreach program： a model for community collaboration to enhance onfarm safety and health.[J]. North Carolina medical journal， 2011，72（6）：471.

[9]張云華，馬九杰，孔祥智，等.農戶采用無公害和綠色農藥行為的影響因素分析——對山西、陜西和山東15縣（市）的實證分析[J].中國農村經濟， 2004（1）：41-50. [ZHANG Yunhua， MA Jiujie， KONG Xiangzhi， et al. The econometric research on the behavior of adopting nonhazard and green pesticides based on the data from fifteen counties in Shanxi， Shaanxi and Shandong Province in China [J].China rural economy， 2004 （1）：41-50.]

[10]褚彩虹，馮淑怡，張蔚文.農戶采用環境友好型農業技術行為的實證分析——以有機肥與測土配方施肥技術為例[J]. 中國農村經濟， 2012（3）：68-77. [CHU Caihong， FENG Shuyi， ZHANG Weiwen. The econometric research on the behavior of adopting environmentfriendly agricultural techniques： taking organic fertilizer and soiltestformulationfertilization as an example[J].China rural economy， 2012（3）：68-77.]

[11]FRLICH M， LECHNER M. Exploiting regional treatment intensity for the evaluation of labor market policies[J]. Journal of the American Statistical Association， 2010，105：1014-1029.

[12]ASFAW S， MITHFER D， WAIBEL H. EU food safety standards， pesticide use and farmlevel productivity： the case of highvalue crops in Kenya[J]. Journal of agricultural economics， 2009，60（3）：645-667.

[13]ALI A， SHARIF M. Impact of farmer field schools on adoption of integrated pest management practices among cotton farmers in Pakistan[J]. Journal of the Asia Pacific economy， 2012，17（3）：498-513.

[14]鐘甫寧，寧滿秀，邢鸝，等.農業保險與農用化學品施用關系研究——對新疆瑪納斯河流域農戶的經驗分析[J]. 經濟學（季刊）， 2007（1）：291-308. [ZHONG Funing，NING Manxiu， XING Li， et al. Study on the relationship between agricultural insurance and agricultural chemical： taking farmers in Manas River Basin as the target[J]. China economic quarterly， 2007（1）： 291-308.]

[15]BEHRMAN J R， DEOLALIKAR A B. The intrahousehold demand for nutrients in rural south India： individual estimates， fixed effects， and permanent income[J]. Journal of human resources， 1990：665-696.

[16]BECKER S O， ICHINO A. Estimation of average treatment effects based on propensity scores[J]. The Stata journal， 2002，2（4）：358-377.

Effects of regional pest control adoption on pesticides application

YING Ruiyao1XU Bin2

（1.Jiangsu New Rural Science and Technology Innovation Research Base， Nanjing Agricultural University，

Nanjing Jiangsu 210095， China； 2.School of Business， Linyi University， Linyi Shandong 276000， China）

AbstractExcessive application of pesticides leads to nonpoint source pollution in agriculture. Regional Pest Control program is one of agricultural socialization services. whether it has decreased quantity of pesticides usage and increased the ratio of environmentalfriendly pesticides application is a worthwhile issue to probe into. It is a selfselection problem whether farmers adopt Regional Pest Control program， which brings about endogenous bias in results. Therefore， this paper utilizes seven provinces nationwide rice production data by propensity score matching method （Propensity Score Matching， PSM）， excluding the selfselection problem of endogenous influence， to evaluate the difference in species and quantity of pesticides application between Regional Pest Control and dispersed insect control by farmers， that is to say， whether this professional plant protection has achieved good environmental benefits and provided safe food for people. The research results show if endogenous influence is ignored， the effect of Regional Pest Control program on alleviation of abuse of pesticides will be overestimated. The quantity of the pesticide usage has been reduced significantly through implementation of Regional Pest Control program， and the ratio of environmentalfriendly pesticides application has been increased. Smallscale growers who adopt the program work better in alleviation of pesticides abuse than largescale growers. It is therefore recommended to raise the level of fiscal support， strengthen epidemic forecasting， improve the pesticides control and increase the programs coverage. Furthermore， the authority should strengthen the guidance and training of professional pest control， encourage land circulation and raise agricultural appropriate scale of operation.

Key wordsRegional Pest Control program； Propensity Score Matching； pesticides application； endogenous problem； food safety