禁/限用促生長(zhǎng)抗生素對(duì)生豬養(yǎng)殖戶收入的影響

楊 宇，徐惠敏，王秀敏，何軍東

(1.邢臺(tái)學(xué)院經(jīng)濟(jì)與貿(mào)易學(xué)院，河北 邢臺(tái) 054001；2.沈陽農(nóng)業(yè)大學(xué)經(jīng)濟(jì)管理學(xué)院，遼寧 沈陽 110861；3.邢臺(tái)學(xué)院農(nóng)業(yè)經(jīng)濟(jì)研究中心，河北 邢臺(tái) 054001；4.邢臺(tái)學(xué)院會(huì)計(jì)學(xué)院，河北 邢臺(tái) 054001)

作為世界畜禽養(yǎng)殖大國，我國對(duì)抗生素非常依賴，然而對(duì)食品動(dòng)物長(zhǎng)期廣泛甚至過量地使用抗生素，導(dǎo)致藥物殘留，環(huán)境污染和高耐藥性細(xì)菌病原體的暴發(fā)性增長(zhǎng)。近年來，世界各國紛紛出臺(tái)日趨嚴(yán)厲的管控措施，控管焦點(diǎn)直指養(yǎng)殖戶對(duì)促生長(zhǎng)用途抗生素(Growth-promoting antibiotics，GPAs)的使用。2006年，歐盟發(fā)布“禁抗令”，2013年美國食品藥品監(jiān)督管理局(FDA)發(fā)布“限抗令”[1]，2017年6月我國農(nóng)業(yè)部發(fā)布了《全國遏制動(dòng)物源細(xì)菌耐藥行動(dòng)計(jì)劃》[2]。我國政府“限抗”的政策趨勢(shì)已十分明確，不過對(duì)GPAs具體的管控思路還可慎重斟酌：是立即禁用風(fēng)險(xiǎn)品種，還是逐步推進(jìn)風(fēng)險(xiǎn)品種退出市場(chǎng)？是分畜禽生長(zhǎng)階段禁用還是全生長(zhǎng)階段禁用？不同的管控思路將對(duì)養(yǎng)殖戶的收入產(chǎn)生影響。本文以生豬養(yǎng)殖為例，重點(diǎn)考察使用GPAs對(duì)養(yǎng)殖戶的經(jīng)濟(jì)影響，試圖揭示養(yǎng)殖戶使用GPAs的動(dòng)因和內(nèi)在機(jī)理。

1 數(shù)據(jù)與方法

1.1 數(shù)據(jù)來源 本文的實(shí)證數(shù)據(jù)來源于2018年1～2月 由課題組組織的“河北省生豬養(yǎng)殖戶抗生素使用行為”的問卷調(diào)查。本次調(diào)查包括預(yù)調(diào)查和正式調(diào)查2個(gè)階段。首先，課題組在邢臺(tái)市開展了小規(guī)模預(yù)調(diào)查，并據(jù)此修正了養(yǎng)殖戶容易出現(xiàn)理解偏差的詞語，改善了某些讓他們感到敏感的問題的提問方式，形成了最終的調(diào)查問卷。接著對(duì)調(diào)研人員進(jìn)行了專門的培訓(xùn)，以保證獲得數(shù)據(jù)的準(zhǔn)確可靠。然后對(duì)河北省石家莊、唐山、保定和邢臺(tái)4個(gè)城市展開正式調(diào)查。具體的抽樣方法是先從每個(gè)樣本市分別選取若干個(gè)村(鎮(zhèn))，再從每個(gè)村(鎮(zhèn))隨機(jī)選取受訪的養(yǎng)殖戶。為了確保回收問卷的質(zhì)量，采用問卷為主、訪談為輔的調(diào)查方式，并且將調(diào)查對(duì)象限定為能夠?qū)ιi養(yǎng)殖發(fā)揮直接決策作用的家庭成員。本次調(diào)查共收集問卷 600 份，剔除關(guān)鍵數(shù)據(jù)缺失或回答不真實(shí)樣本后，共獲得有效問卷476份，有效率為79.3%。

1.2 變量設(shè)定與賦值 變量設(shè)定與賦值說明如表1所示，由于部分GPAs是由供應(yīng)商提供的飼料里面添加的，不少養(yǎng)殖戶無法確定促生長(zhǎng)用途抗生素的具體種類、數(shù)量以及支出成本，因此本文借鑒歐美學(xué)者的處理方法[3-5]，在問卷中，受訪養(yǎng)殖戶被問到：“是否向保育豬和育肥豬在其飼料或水中添加過抗生素？”回答選項(xiàng)為“是”或“否”?；卮稹笆恰钡酿B(yǎng)殖戶還被問及在保育和育肥2個(gè)不同生長(zhǎng)階段添加的抗生素被用于以下哪種目的？回答選項(xiàng)包括促進(jìn)生長(zhǎng)、預(yù)防疾病、治療疾病(可多選)。在本文中，GPAs的使用者指的是為達(dá)到促進(jìn)生豬生長(zhǎng)的目的使用抗生素的養(yǎng)殖戶。McBride(2008年)采用全要素生產(chǎn)率來衡量美國養(yǎng)殖戶的生豬生產(chǎn)力[5]，類似地本文也采用全要素生產(chǎn)率，即養(yǎng)殖戶出欄生豬重量與全部生產(chǎn)投入之比，單位是kg/元。2012-2016年，我國生豬養(yǎng)殖戶仔豬投入、飼料投入和勞動(dòng)力投入三者合計(jì)占生豬養(yǎng)殖戶總成本的比重平均為96.46%[6]，因此本文選擇這3個(gè)指標(biāo)作為成本變量來表示生產(chǎn)投入情況。模型中規(guī)定的外生變量包括個(gè)體特征和養(yǎng)殖場(chǎng)特征，以及一套生豬養(yǎng)殖操作。個(gè)體特征包括年齡、受教育程度和養(yǎng)殖經(jīng)驗(yàn)，以說明養(yǎng)殖戶知識(shí)與能力的差異。養(yǎng)殖場(chǎng)特征首先考慮了生豬養(yǎng)殖結(jié)構(gòu)，包括養(yǎng)殖規(guī)模和養(yǎng)豬專業(yè)化程度，其中養(yǎng)殖規(guī)模借鑒了孫世民等[7]與王建華等[8]的劃分標(biāo)準(zhǔn)?？赡苡绊慓PAs使用決策的其他養(yǎng)殖場(chǎng)特征還包括供應(yīng)合同、豬舍類型以及生豬買家是否要求禁用抗生素。通過供應(yīng)合同，供應(yīng)商可能提供含有GPAs的飼料給養(yǎng)殖戶；豬舍類型(是否封閉)表明生豬被限制活動(dòng)的程度，表示疾病傳播潛在風(fēng)險(xiǎn)大??；生豬購買者(或銷售商)是否要求生豬在任何時(shí)候都不喂養(yǎng)抗生素。為了隔離生產(chǎn)力與GPAs使用決策之間的關(guān)聯(lián)，產(chǎn)出方程中還增加了一些其他變量，包括蠕蟲病防控、全進(jìn)全出管理、鼠害防控、運(yùn)輸車輛消毒、分性別飼喂。

1.3 養(yǎng)殖戶組間差異性檢驗(yàn) 為了考察使用GPAs和沒使用GPAs養(yǎng)殖戶的差異性，我們將樣本分為2組，通過t檢驗(yàn)對(duì)組間均值進(jìn)行比較，當(dāng)P≤0.1時(shí)，否定原假設(shè)，表明2組均值存在顯著差異。t檢驗(yàn)結(jié)果表明，2組養(yǎng)殖戶在多個(gè)變量上存在顯著差異，即存在自選擇問題。

在生豬保育階段，2組養(yǎng)殖戶大多數(shù)變量的均值存在顯著差異。在產(chǎn)出與投入方面，無論是單位收入還是仔豬、飼料和勞動(dòng)力投入，GPAs用戶均值都顯著高于非GPAs用戶。在個(gè)體特征方面，使用GPAs的養(yǎng)殖戶通常年齡更大，受教育程度更高，養(yǎng)殖經(jīng)驗(yàn)更豐富。在經(jīng)營(yíng)特征方面，變量供應(yīng)合同、養(yǎng)殖規(guī)模和豬舍類型的均值差異顯著，表明更大生產(chǎn)規(guī)模、簽訂供應(yīng)合同和封閉式豬舍的養(yǎng)殖場(chǎng)，使用GPAs更為普遍。在養(yǎng)殖操作特征方面，使用GPAs的養(yǎng)殖戶，更可能采用蠕蟲病防控，全進(jìn)全出管理，鼠害防控和運(yùn)輸車輛消毒等措施。

在生豬育肥階段，同保育階段相比，很多變量的均值不再存在顯著差異。在產(chǎn)出與投入方面，沒有變量存在顯著差異。在個(gè)體特征方面，變量養(yǎng)殖經(jīng)驗(yàn)差異變得不再顯著。在經(jīng)營(yíng)特征方面，變量供應(yīng)合同與養(yǎng)殖規(guī)模也不再有顯著差異。養(yǎng)殖操作特征方面，變量運(yùn)輸車輛消毒不再存在顯著差異。

表1 變量設(shè)定及賦值說明

1.4 模型構(gòu)建 本文采用Maddala(1983年)提出的處理效應(yīng)模型來解決由不可觀測(cè)變量導(dǎo)致的自選擇問題[9]。該模型估計(jì)2個(gè)方程，方程1是產(chǎn)出方程，解釋養(yǎng)殖戶養(yǎng)殖收入的決定因素，而變量使用GPAs包括在產(chǎn)出方程的解釋變量中；方程2是選擇方程，運(yùn)用Probit估計(jì)養(yǎng)殖戶使用GPAs的概率。通過引入干擾項(xiàng)的聯(lián)合分布來校正可能存在的樣本選擇偏差，允許無偏估計(jì)使用GPAs對(duì)養(yǎng)殖收入的影響。首先構(gòu)建1個(gè)簡(jiǎn)單的柯布 - 道格拉斯形式的生產(chǎn)函數(shù)：

yi=α0+β1x1i+β2x2i+β3x3i+δDi+εi

(1)

其中yi是因變量，表示產(chǎn)出水平；x1i表示第i個(gè)養(yǎng)殖戶生產(chǎn)性投入特征向量，包括仔豬投入、飼料投入和勞動(dòng)力投入3個(gè)變量，由實(shí)際投入數(shù)值取對(duì)數(shù)后獲得；x2i表示第i個(gè)養(yǎng)殖戶個(gè)體特征向量，包括養(yǎng)殖戶年齡、受教育程度和養(yǎng)殖經(jīng)驗(yàn)3個(gè)變量；x3i表示第i個(gè)養(yǎng)殖戶經(jīng)營(yíng)特征向量，包括專業(yè)化程度、是否簽有訂貨合同、豬舍類型和養(yǎng)殖規(guī)模4個(gè)變量；x4i表示第i個(gè)養(yǎng)殖戶養(yǎng)殖操作特征向量，包括Di代表第i個(gè)養(yǎng)殖戶是否使用GPAs，是二元分類虛擬變量。產(chǎn)出yi用養(yǎng)殖戶平均每頭生豬的收入來表示，產(chǎn)出yi是實(shí)際數(shù)值取對(duì)數(shù)后獲得的，ε表示隨機(jī)誤差項(xiàng)。

通過前面的t檢驗(yàn)，發(fā)現(xiàn)養(yǎng)殖戶使用GPAs存在自選擇問題。由于養(yǎng)殖環(huán)境條件差，管理水平低的養(yǎng)殖戶能夠從使用GPAs中獲得更大的生產(chǎn)力提升，因此這類養(yǎng)殖戶更傾向于使用GPAs，而且環(huán)境條件和管理水平等不可觀測(cè)變量也可能影響生豬生產(chǎn)力，所以方程(1)無法被直接估計(jì)，這種情況下我們引入潛變量Di*，構(gòu)建以下決定方程。

(2)

當(dāng)Di=1時(shí)，Di*≥0；當(dāng)Di=0時(shí)，Di*<0。其中向量Zi可以與zi有重疊變量，但至少有1個(gè)變量，比如Z1i不在zi中，且滿足Cov(Z1i，εi)，此時(shí)Z1i可以視為Di的工具變量。這里將買家要求禁止使用抗生素作為工具變量Z1i，而且假設(shè)Z1i直接影響Di養(yǎng)殖戶使用GPAs，但并不直接影響結(jié)果變量yi(只通過Di間接影響yi)。依據(jù)處理效應(yīng)模型，定義逆米爾斯函數(shù)為：

(3)

這樣我們就得到了對(duì)使用GPAs與未使用GPAs的養(yǎng)殖戶都適用的預(yù)期產(chǎn)出方程。

Ε[yi|xi,Zi]=βxi+γzi+δDi+ρσλi

(4)

這里采用極大似然估計(jì)法，同時(shí)估計(jì)所有模型參數(shù)。如果ρ=0，則不存在內(nèi)生性，可直接做方程(1)的最小二乘估計(jì)；如果ρ≠0，這正是模型內(nèi)生性選擇偏誤的來源。當(dāng)ρ<0時(shí)，最小二乘估計(jì)將低估使用GPAs對(duì)生豬生產(chǎn)力的貢獻(xiàn)；當(dāng)ρ>0時(shí)，使用GPAs對(duì)生豬生產(chǎn)力的貢獻(xiàn)將被高估。

1.5 統(tǒng)計(jì)分析 采用STATA15.0軟件對(duì)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，使用極大似然法對(duì)處理效應(yīng)模型進(jìn)行估計(jì)。用“P”表示使用GPAs及其對(duì)養(yǎng)殖收入的影響水平，其中P<0.1表示影響基本顯著；P<0.05表示影響比較顯著，P<0.01表示影響極顯著。

2 結(jié)果

為了探究在不同養(yǎng)殖階段使用GPAs對(duì)生豬生產(chǎn)的影響，我們分別估計(jì)養(yǎng)殖戶對(duì)保育豬和育肥豬是否使用GPAs及其對(duì)養(yǎng)殖收入的影響。表2分別顯示了對(duì)保育豬和育肥豬的選擇方程和產(chǎn)出方程的估計(jì)結(jié)果。

2.1 對(duì)選擇方程的估計(jì) 無論在生豬保育階段還是育肥階段，對(duì)選擇方程的估計(jì)有較多相似之處。養(yǎng)殖戶年齡越大，越傾向于使用GPAs。受教育程度高的養(yǎng)殖戶，更傾向于使用GPAs，擁有封閉式豬舍的養(yǎng)殖戶更可能使用GPAs。變量買方要求禁用抗生素對(duì)使用GPAs存在顯著的負(fù)向影響。變量養(yǎng)殖經(jīng)驗(yàn)、專業(yè)化程度、供應(yīng)合同以及養(yǎng)殖規(guī)模對(duì)使用GPAs沒有產(chǎn)生顯著影響。對(duì)選擇方程的估計(jì)結(jié)果的不同之處在于：在生豬育肥階段，變量養(yǎng)殖經(jīng)驗(yàn)對(duì)使用GPAs存在基本正向影響(只在10%統(tǒng)計(jì)水平上顯著)，而在保育階段則影響不顯著。

2.2 對(duì)產(chǎn)出方程的估計(jì) 如表2所示，在產(chǎn)出方程中，生豬保育階段和育肥階段各個(gè)變量的估計(jì)有很多相似之處。其中變量飼料投入、勞動(dòng)力投入、養(yǎng)殖經(jīng)驗(yàn)、養(yǎng)殖規(guī)模、蠕蟲病防控、全進(jìn)全出管理和運(yùn)輸車輛消毒對(duì)養(yǎng)殖收入都具有顯著的正向影響。結(jié)果顯示，飼料投入與勞動(dòng)力投入對(duì)生產(chǎn)力存在顯著的正向影響，但是仔豬投入對(duì)養(yǎng)殖收入影響不顯著。養(yǎng)殖經(jīng)驗(yàn)與生產(chǎn)力顯著正相關(guān)，養(yǎng)殖規(guī)模對(duì)養(yǎng)殖收入的正向影響，蠕蟲病防控、全進(jìn)全出管理和運(yùn)輸車輛消毒對(duì)生豬養(yǎng)殖收入的正向影響。

在產(chǎn)出方程中，對(duì)保育階段和育肥階段各個(gè)變量的估計(jì)也有不同之處。一是在保育階段，變量養(yǎng)殖戶年齡對(duì)生豬養(yǎng)殖收入產(chǎn)生了顯著的負(fù)向影響，在育肥階段影響則不顯著；二是在育肥階段，供應(yīng)合同對(duì)生豬養(yǎng)殖收入存在顯著正向影響(只在10%統(tǒng)計(jì)水平上顯著)，在保育階段影響則不顯著。還有一處不同也正是本文考察的重點(diǎn)，即保育豬和育肥豬產(chǎn)出方程中變量添加GPAs的系數(shù)差異。在生豬保育階段使用GPAs與生豬養(yǎng)殖收入存在顯著的正相關(guān)關(guān)系。相反，在生豬育肥階段使用GPAs，對(duì)產(chǎn)出的影響不顯著。

表2 生豬養(yǎng)殖戶處理效應(yīng)模型的極大似然估計(jì)

生豬保育階段的選擇方程與產(chǎn)出方程誤差的相關(guān)性估計(jì)(Rho系數(shù))表明，兩者具有顯著的負(fù)相關(guān)性。這一結(jié)果意味著負(fù)選擇偏差，表明如果沒有考慮選擇偏差，那么GPAs選擇對(duì)保育階段生豬生產(chǎn)力的影響將會(huì)被低估。相反，對(duì)于育肥豬來說，2個(gè)方程之間誤差的相關(guān)性并不顯著，表明兩者不存在選擇偏差。

2.3 穩(wěn)健性檢驗(yàn) 為了評(píng)估結(jié)果的穩(wěn)健性，我們使用亞治療劑量抗生素(Sub-therapeutic antibiotics，STAs)作為因變量，重新估計(jì)模型。這是因?yàn)樵蜃兞縂PAs被定義為用于促進(jìn)生長(zhǎng)的抗生素，而促進(jìn)生長(zhǎng)和預(yù)防疾病用途的抗生素均通過添加到飼料里飼喂動(dòng)物，部分受訪者可能混淆添加在飼料里抗生素的用途，但是對(duì)是否添加抗生素記憶更加準(zhǔn)確。

表3顯示了穩(wěn)健性檢驗(yàn)的估計(jì)結(jié)果，僅包括在生豬保育階段和育肥階段，STAs變量對(duì)生豬養(yǎng)殖收入影響的估計(jì)。在保育階段STAs變量的系數(shù)是0.024 7，并且在1%統(tǒng)計(jì)水平上顯著，表明使用STAs對(duì)生豬養(yǎng)殖收入產(chǎn)生顯著的正向影響。在育肥階段使用STAs對(duì)生豬養(yǎng)殖收入依然不存在顯著影響。表明用GPAs作為因變量，結(jié)果是穩(wěn)健的。

表3 穩(wěn)健性檢驗(yàn)

3 討論

國內(nèi)學(xué)者主要對(duì)獸藥使用行為做了相關(guān)研究，大多從個(gè)體特征[10]、認(rèn)知特征[11-12]、渠道協(xié)作[13]、政府規(guī)制[14]等角度探討了獸藥使用行為影響因素以及各影響因素之間層級(jí)關(guān)系[7]。除了王建華等[8]基于獸藥使用的不同情境，將獸藥分為治療類獸藥和保健類獸藥，國內(nèi)學(xué)者一般在研究中沒有按照用途對(duì)獸用抗生素使用行為進(jìn)行區(qū)分，而針對(duì)養(yǎng)殖戶GPAs使用行為方面的研究目前沒有發(fā)現(xiàn)。

國外學(xué)者專門針對(duì)GPAs使用行為做了細(xì)致研究，主要包括GPAs對(duì)生豬產(chǎn)出的影響[15]，GPAs對(duì)生豬生產(chǎn)性能其他方面的貢獻(xiàn)[16]，禁/限抗令對(duì)生豬養(yǎng)殖業(yè)產(chǎn)業(yè)結(jié)構(gòu)的影響[17]，減量使用GPAs的激勵(lì)與懲罰機(jī)制[18]等方面。但大都忽略了自選擇性問題,即是否使用GPAs是養(yǎng)殖戶自我選擇的結(jié)果，那些年齡大的養(yǎng)殖戶，衛(wèi)生條件和管理水平比較差的養(yǎng)殖場(chǎng)往往更傾向于使用GPAs，因?yàn)樗麄兺ㄟ^使用GPAs能獲得更大的生產(chǎn)力提升[4]。由于某些因素(比如環(huán)境條件，管理水平等因素)是不可觀測(cè)的，并且對(duì)使用GPAs以及生產(chǎn)力都存在影響。在這種情況下，采用最小二乘法所得到的估計(jì)量經(jīng)常是有偏且不一致[19]。

本文以我國養(yǎng)殖戶GPAs使用行為作為研究對(duì)象，運(yùn)用處理效應(yīng)模型解決由不可觀測(cè)變量造成的自選擇性問題，并且采用現(xiàn)實(shí)中的調(diào)查數(shù)據(jù)來驗(yàn)證Dritz等[14]通過試驗(yàn)得到的結(jié)論。因此，本試驗(yàn)在內(nèi)容和數(shù)據(jù)處理方法上有一定新意。

試驗(yàn)結(jié)果表明，在保育階段使用GPAs可顯著提高生豬養(yǎng)殖收入，這一結(jié)果獲得了穩(wěn)健性檢驗(yàn)的支持。而且極大似然估計(jì)揭示出在保育階段存在負(fù)選擇性偏誤，如果不考慮使用與沒使用GPAs養(yǎng)殖戶之間的差異，那么使用GPAs對(duì)生豬養(yǎng)殖收入的影響將會(huì)被低估。在保育階段使用GPAs可以有效彌補(bǔ)養(yǎng)殖戶在管理水平、衛(wèi)生條件以及其他不可觀測(cè)因素方面存在的差異，那些因?yàn)榄h(huán)境條件和管理水平較差而生產(chǎn)效率低下的養(yǎng)殖戶通過使用GPAs能獲得最大的生產(chǎn)力提升。反之，如果在生豬保育階段禁止使用GPAs會(huì)降低生豬養(yǎng)殖收入，增加生豬養(yǎng)殖成本。短期來看這些成本可能是由于豬的死亡率增加和產(chǎn)量下降造成的，長(zhǎng)期來看，可能來自生豬養(yǎng)殖戶生物安全管理的改變和其他投入的必要調(diào)整。另外，在育肥階段使用GPAs對(duì)生豬養(yǎng)殖戶收入未發(fā)現(xiàn)顯著影響。需要指出的是，本文的局限在于僅探究了使用GPAs和生豬養(yǎng)殖收入之間的關(guān)系，并未考察使用GPAs對(duì)生豬生產(chǎn)可能造成的其他影響，例如GPAs的使用可能會(huì)減少生豬養(yǎng)殖收入變化，進(jìn)而達(dá)到降低養(yǎng)殖風(fēng)險(xiǎn)的目的。關(guān)于GPAs對(duì)生產(chǎn)力波動(dòng)影響的進(jìn)一步研究將有助于解釋為什么GPAs在生豬育肥階段被廣泛使用。