基于DID模型的福建省生態(tài)公益林保護對林農的經濟影響評估

陳 欽，黃巧龍，游玲娜，謝麗星，程秋旺，潘 輝

（1.福建農林大學 經濟學院，福州 350002；2福建農林大學 管理學院，福州 350002；3.閩江學院，福州 350108）

隨著福建省鐵路、公路建設的不斷擴大，水源保護區(qū)域的增加，生態(tài)保護的加強，有些商品林被政府劃為生態(tài)公益林，有些商品林被政府界定為重點生態(tài)區(qū)位商品林，禁止采伐。重點生態(tài)區(qū)位商品林雖然形式上屬于商品林，但是經濟實質上雷同于生態(tài)公益林。這會對林農產生經濟影響。本文從林農的微觀角度出發(fā)，分析福建省生態(tài)公益林保護對林農的經濟影響。從理論研究來看，針對生態(tài)公益林保護對林農經濟影響的研究較少，未見報道采用大樣本和DID模型進行研究；從實踐方面來看，本文在問卷調查的基礎上分析福建省生態(tài)公益林保護對林農的經濟影響以及影響金額，對政府制定福建省生態(tài)公益林生態(tài)補償價格，保障林農的切身利益，促進林業(yè)生態(tài)效益和經濟效益的提高，進而推進生態(tài)公益林保護，具有重要的實踐意義。

1 模型設定

本文使用DID方法來評估生態(tài)公益林保護對林農的經濟影響，用林農家庭林業(yè)收入作為經濟影響的衡量指標。DID模型的基本思想是把樣本分為控制組（Control）與處理組（Treated）(史桂芬,2012)[1]，其中，控制組為未受到生態(tài)公益林保護政策影響的樣本，處理組為受到生態(tài)公益林保護政策影響的樣本，將兩組進行比較從而評估該政策對林農林業(yè)收入的影響，首先收集兩個時期的數據，即2002年和2013年相關變量數據；其次將樣本分為兩組，2002年只有商品林，2013年擁有生態(tài)公益林的樣本為處理組；而2002年和2013年都只有商品林，沒有生態(tài)公益林的樣本為控制組。通過實證分析，分別得到處理組和控制組在2002年和2013年林業(yè)收入的變化值，這兩個變化值的差額就反映了生態(tài)公益林保護對林農林業(yè)收入的凈影響(古曉,2013)[2]。

在沒有考慮其他控制變量的情況下，DID模型的表達式如下：

其中，i代表林農，t代表時期，在t=0階段，所調查的林農的森林是商品林；在t=1階段，上述一部分林農的森林由商品林變成生態(tài)公益林，而另一部分林農的森林仍然是商品林。關鍵解釋變量FPCit是一組二值虛擬變量，代表t時期i林農情況，分別測量林農的森林是否由商品林變成生態(tài)公益林，1表示是，0表示否。Yit是t時期i林農的收入；t是時期虛擬變量，生態(tài)公益林保護實施前取值為0，實施后取值為1。ci是非觀測效應，控制隨時間不變的不可觀測因素；μit是特異性擾動項，代表因林農因時而變且影響Yit的那些非觀測擾動因素。α0、α1和δ是待估計參數，其中δ是研究最為關注的參數，衡量了生態(tài)公益林保護對變量Yit的經濟影響。

根據該模型，可以分別得到處理組和控制組在生態(tài)公益林保護實施前后林農林業(yè)收入變動的模型：當FPCit=0，t=0時，得到控制組和處理組在生態(tài)公益林保實施前林業(yè)收入的變動模型；當FPCit=0，t=1時，得到控制組在生態(tài)公益林保實施后林業(yè)收入的變動模型；當FPCit=1，t=1時，得到處理組在生態(tài)公益林保實施后林業(yè)收入的變動模型，具體如下所示：

根據式（2），可得到控制組在生態(tài)公益林保護實施前后林農林業(yè)收入的變化值為：

根據式（3），可得到處理組在生態(tài)公益林保護實施前后林農林業(yè)收入的變化值為：

根據式（4）和（5），生態(tài)公益林保護對林農林業(yè)收入的凈影響為：

為了剔除其它因素對林業(yè)收入的影響，本文進一步引入控制變量，采用非觀測效應綜列數據模型如下：

式中，Zit是一組隨時間變化的可觀測的影響因變量Yit的控制變量；Xit是一組不隨時間變化，或隨時間同等變化的可觀測的影響因變量Yit的控制變量；β和γ是待估參數矩陣，衡量了控制變量對因變量Yit的作用。

2 理論分析

經濟學認為在經濟活動中，人們都是理性的，人們參加經濟活動時都持有預期收益大于預期成本的觀點，這是人們參與活動的基本出發(fā)點和動力。人們投入的主要目的是獲得收益，只有當其收益大于成本時，才是經濟合理的。一般而言，農戶作為理性的經濟主體，其生態(tài)公益林建設和保護是為了追求潛在的收益，即通過比較其成本和收益水平，在可接受的范圍內，產生行為。由此可見，在生態(tài)公益林建設和保護過程中，成本和收益水平是農戶決策的核心與關鍵。

2.1 生態(tài)公益林保護對林農林業(yè)收入的影響

生態(tài)公益林保護對林農林業(yè)收入的影響主要體現在商品林和生態(tài)公益林的收益差距上。為了保護森林資源，提高生態(tài)環(huán)境質量，政府采用自上而下的辦法將部分商品林強制劃定為生態(tài)公益林，禁止對其進行采伐，導致林業(yè)收入減少(沈潔,2014)[3]。為此，政府制定了生態(tài)補償標準，生態(tài)公益林的經營者可以獲得一定的經濟補償，但是現行的生態(tài)補償標準低，對于福建省絕大部分生態(tài)公益林而言明顯不夠，使得很多林農獲得的補償價不能彌補營林成本。

2.2 收入影響因素的研究綜述

農業(yè)收入會受到農業(yè)生產條件、生產經營費用、勞動力等因素的影響；同時，根據公共財政理論，財政支農也會影響農業(yè)收入(楊靜,2014)[4]。高陽、趙正等(2014)運用實證方法分析林農林業(yè)收入的影響因素，結果表明：林農林業(yè)收入的影響因素包括受教育程度、林下經營項目、林種等[5]。林和平(2009)分別分析了林地資源數量、林地資源經營方式與林業(yè)收入之間的關系，得出：耕地資源、林地資源、人力資源以及樹種等都對林業(yè)收入有顯著的影響[6]。劉偉平(2009)認為，經濟環(huán)境和價格等經濟因素、勞動力數量、林業(yè)政策等都會影響林農的林業(yè)收入[7]。石康橋、蘇時鵬等(2014)通過比較集體林權制度改革前后林農林業(yè)收入的影響因素，得出：林改前林農林業(yè)收入的主要影響因素包括林地投入、經營類型等；而林改后林農林業(yè)收入的主要影響因素為受教育程度、林地面積等[8]。姚林香、舒成(2010)通過對江西7縣農戶樣本數據的分析，發(fā)現：農業(yè)人口數、耕地面積等對農民收入具有顯著影響[9]。

2.3 變量的選擇及說明

本文關注的是生態(tài)公益林保護對林農的經濟影響，選取林業(yè)收入作為模型的因變量。模型選取的自變量包括虛擬變量、隨時間變化的可觀測的影響林業(yè)收入的控制變量以及不隨時間變化，或隨時間同等變化的可觀測的影響林業(yè)收入的控制變量。綜合有關學者的研究成果以及理論分析，首先，本文選取時期、商品林是否變成生態(tài)公益林作為虛擬變量。如果由于高速公路、高鐵、水庫等建設，林農的商品林被劃為生態(tài)公益林，禁止采伐，林農就沒有木材采伐收入，為此林農可以獲得森林生態(tài)效益補償收入、生態(tài)公益林林下經營補貼收入等，可見，商品林是否變成生態(tài)公益林對林農林業(yè)收入產生影響，本文試圖研究這個關鍵變量的影響大小。其次，選取家庭勞動力、人均GDP、職務、受教育程度、耕地面積、林地面積和樹種作為控制變量，具體說明如下：①家庭勞動力。家庭勞動力越多，林業(yè)生產要素中的人力資源數量越多，所以林業(yè)收入越多，假設家庭勞動力與林業(yè)收入正相關；②人均GDP。人均GDP越高，表明當地經濟發(fā)展水平越高，木材和林產品的銷售量越大，所以林農的林業(yè)收入越多，假設人均GDP與林業(yè)收入正相關；③林農職務。擔任職務（如村干部）的林農，一般社會資本更多，更容易獲得木材采伐指標和更快了解到林業(yè)信息，林業(yè)收入也會更高，假設林農擔任的職務與林業(yè)收入正相關；④受教育程度。受教育程度越高，有利于林農及時地對政策和市場作出反應，林農林業(yè)經營能力越強，林業(yè)收入越高，假設受教育程度與林業(yè)收入正相關；⑤耕地面積。耕地面積越大，對耕地投入越多，農業(yè)收入越多，對林業(yè)的經濟依賴越小，對林業(yè)的投入隨之減少，導致林業(yè)收入減少，假設耕地面積與林業(yè)收入負相關；⑥家庭林地面積，一般情況下林地面積越多的林農，林業(yè)收入越高，假設林地面積與林業(yè)收入正相關；⑦樹種。不同樹種的經濟價值不一樣，一個樹種經濟價值越大，其經濟收益越高，假設樹種的經濟價值與林業(yè)收入正相關。變量的選擇及其賦值說明具體如表1所示。

3 計量結果及其分析

考慮到可能存在影響林業(yè)收入的其他因素，為了剔除這些因素對林業(yè)收入的影響，本文選擇人均GDP、家庭勞動力、林農職務、受教育年限、家庭耕地面積、林地面積、樹種作為控制變量，分析福建省生態(tài)公益林保護對林農林業(yè)收入的凈影響。運用Stata12.0軟件，將上述控制變量放入模型中進行估計，得出引入控制變量后DID模型的計量結果，如表2和表3所示，模型的Prob＞F=0.000 0，說明模型整體顯著；R方為0.631 7，調整后R方為0.624 0，說明模型擬合度較好。

表2 DID模型計量結果一

從表2可以看出，DID值的系數為-8 700，且在10%的統(tǒng)計水平上顯著，說明在控制其他因素的情況下，福建省生態(tài)公益林保護使樣本林農戶均林業(yè)收入降低了8 700元，也就是說福建省生態(tài)公益林保護對林農的經濟影響為-8 700元。政府確定生態(tài)公益林補償價格可以參考這個數字，并且據此計算每年每畝生態(tài)公益林補償標準，從而為政府修訂生態(tài)公益林補償政策提供重要依據。可見，政府應重視生態(tài)公益林補償政策的完善和實施，彌補生態(tài)公益林保護給林農帶來的經濟損失，從而促進生態(tài)公益林保護工作的有效進行。

從表3可知，時期P值為0.000＜0.01，說明2013年與2002年相比，樣本林農林業(yè)收入發(fā)生了顯著變化；人均GDP對林業(yè)收入影響的回歸系數為正，且P值為0.002＜0.01，說明人均GDP對林農林業(yè)收入在1%水平上具有正向的顯著影響；家庭耕地面積的回歸系數為負，且P值為0.071＜0.1，說明耕地面積對林業(yè)收入在10%水平上具有負向的顯著影響；家庭林地面積的回歸系數為正，且P值為0.000＜0.01，說明林地面積對林業(yè)收入在1%水平上具有正向的顯著影響。

表3 DID模型計量結果二