城鎮化與居民乳制品消費
——基于省際動態面板的實證研究

鄭曉冬，方向明，賀子軒

(中國農業大學經濟管理學院，北京 100083)

改革開放以來，隨著大量農村勞動力涌入城市，我國的城鎮化率從1978年的17.9%上升至2014年的54.8%。在這一過程中，居民的食物消費模式發生了明顯變化，主要表現為糧食消費量逐步下降，動物性食品消費量不斷提高[1-2]。就乳制品而言，按照1992年不變價，2011年城鎮人均乳制品消費支出是1992年的5.6倍(附圖)，2015年全國人均乳制品消費量達每年36kg。乳制品營養豐富，是城鎮居民食物消費中的重要組成部分，然而目前居民乳制品消費仍處于較低水平，2011年城鎮居民乳制品消費率僅為25.6%[3]，2015年人均乳制品消費量僅為世界平均水平的1/31。因此，居民乳制品消費仍有較大潛力，探究乳制品消費的主要決定因素不僅有助于更好地促進居民的營養攝入與健康狀態，同時對于我國奶業的發展具有重要政策意義。

現有關于居民乳制品消費的研究主要集中在以下幾個方面，分別是居民乳制品消費特征討論[4-5]、乳制品消費相關影響因素研究[6-7]、食品安全事件對居民乳制品消費行為的影響分析[8-9]以及對未來乳制品消費的預測[10]等，然而直接論述與檢驗城鎮化與居民乳制品消費關系的研究較為少見。那么城鎮化對城鎮居民乳制品消費的影響如何？影響機制又是怎樣？這是本文所關注的重點問題，而回答這些問題能為更好地促進居民健康狀態和乳制品消費市場的發展提出相關對策建議。

1 文獻回顧與理論分析

梳理文獻發現，城鎮化主要通過收入水平、就業結構以及消費可及性等3方面對城鎮居民乳制品消費水平產生作用。第一，收入水平。城鎮化進程中，大量農村勞動力流入城市，在市場化條件下，城鎮居民將通過與農村移民進行分工合作產生的互補效應使城鎮居民獲得更大的工資收益[11]。郭克良[12]基于河北省數據的研究結果表明，城鎮化對城鎮居民收入有顯著正向影響。凱恩斯的消費理論認為人們的消費是收入水平的函數，收入的提高將直接增強居民食物消費的能力，這無疑將有助于提升城鎮居民乳制品消費水平；第二，就業結構。我國的城鎮化進程也是第三產業快速發展的過程，由于城鎮勞動力對城市的主要勞動力市場有壟斷性地位[13]，加之農村人口在教育程度與技能處于劣勢，農村勞動力在城市將補足原本由城鎮居民從事的第二產業(比如建筑業)崗位[14]，城鎮居民的就業結構將偏向非體力勞動的第三產業，從而在一定程度上轉變食物消費習慣，城鎮居民對主食類高能量食物需求將減小[15]，對較具營養價值的乳制品將有更大的消費需求。第三，消費的可及性。人口向城市聚集的過程也促使運輸零售終端不斷向城市集中，便利的消費渠道將促進城鎮居民乳制品的購買與消費[6]。從現有統計數據看，1998年全國連鎖零售門店總數為1.92萬個，2011年為19.58萬個，后者是前者的10倍左右。而20世紀90年代和21世紀初也是我國城鎮化快速發展的階段。胡定寰等[16]和Fuller 等[17]對昆明、北京、上海和廣州的消費者的調查發現，大賣場及連鎖便利商店等對居民乳制品的購買行為有顯著的正向影響。因此城鎮化進程中基礎設施和零售網絡系統的完善也將有利于促進城鎮居民乳制品消費狀況。附圖展示了我國城鎮化、城鎮居民乳制品消費以及上述三方面影響途徑的近20年的時間變化趨勢，可以發現，居民收入水平、就業第三產業GDP占比以及零售業相關指標與城鎮化和居民乳制品消費有較明顯的共同趨勢。

根據以上理論分析，本文基于1992—2011年省際面板數據，運用系統廣義矩估計克服可能存在的內生性問題，從宏觀角度討論城鎮化對居民乳制品消費的影響及影響機制。此外，對城鎮化水平進行分組探析城鎮化及其影響機制的異質性，從而為進一步促進居民乳制品消費水平以及奶業發展提供經驗證據和對策建議。

附圖 城鎮化與城鎮居民乳制品消費的相關統計圖示注：數據來自歷年《中國統計年鑒》和《中國奶業年鑒》。城鎮乳制品消費指數和收入指數是以1992年為基期(100)，以1992年不變價的歷年城鎮乳制品消費和人均可支配收入與基期值之比乘以100%得出。由于就業結構與產業結構類似，這里以第三產業GDP占比變化代表就業結構的變化；消費可及性由全國連鎖零售門店總數、營業面積和從業人數表征，限于數據可得性，部分年份數據有所缺失

2數據來源：全球乳糖不耐受癥比例分布圖https：//en.wikipedia.org/wiki/Lactose_intolerance#cite_note-GHR2016～63

2 模型、變量與數據

2.1 模型設定

設定簡約式消費模型作為城鎮化與城鎮居民乳制品消費關系的基準估計模型，具體面板回歸模型如式(1)：

DCit=α0+α1URit+α2Xit+μi+εit

(1)

式(1)中，DC為本文的因變量，表示城鎮居民乳制品消費水平，下標i和t分別代表第i個地區和第t期。UR為關鍵自變量，表示城鎮化水平。Xit為一系列控制變量，ui表示不隨時間變化的個體效應，εit為隨機擾動項。由于個體的消費習慣存在著一定的時間慣性[18]，因而當期居民乳制品消費水平很可能受到上一期消費狀況的影響。據統計，中國超過80%的居民有不同程度的乳糖不耐受癥2，因此乳制品消費可能存在較強的棘輪效應。為此，在自變量中加入了因變量的滯后項反應乳制品消費的慣性，如式(2)所示：

DCit=β0+β1DCi，t-1+β2URit+β3Xit+μi+εit

(2)

式(2)中，DCi，t-1表示各地區滯后一期的城鎮居民乳制品消費狀況，其余變量含義與式(1)相同。由于因變量的滯后項與擾動項存在相關性，直接使用靜態面板固定效應模型所得的估計結果是有偏不一致的[19]。在這種情況下，需采用動態面板估計方法，即差分GMM或者系統GMM來克服這一內生性問題。由于差分GMM可能存在的弱工具變量問題，本文采用系統GMM法，將差分和水平變量的滯后項作為工具變量進行估計，使所得結果更有效率[20]。然而系統GMM的有效估計也需滿足擾動項差分無序列相關以及工具變量外生可靠等條件。因此本文在運用以上模型進行估計的同時，也進行了二階序列相關檢驗以及Sargan檢驗用以保證估計結果的可靠性。

而后按照中介效應模型的一般步驟[21]，通過加入中介變量后關鍵自變量的系數與顯著性變化檢驗中介變量的作用。加入中介變量的模型式(3)所示：

DCit=γ0+γ1DCi，t-1+γ2URit+γ3MEit+γ4Xit+μi+εit

(3)

式(3)中，ME為本文設定的中介變量，包括收入水平、就業結構與消費可及性。就本研究而言，如城鎮化對乳制品消費的影響顯著，而加入中介變量后，中介變量顯著影響城鎮居民乳制品消費的同時，城鎮化變量的系數減小且顯著性下降，則說明城鎮化通過中介變量影響城鎮居民乳制品消費。

2.2 變量說明與數據來源

本文所用數據為1992—2011年3中國內地29個省(市、區)20個年度的省際面板數據，由于西藏和重慶數據存在缺失，按照多數文獻的做法，將西藏的數據舍棄，同時將重慶的數據并入四川。

2.2.1 因變量 本文的因變量為城鎮居民乳制品消費水平，該指標由運用歷年各地區消費價格指數(CPI)調整為以1992年不變價計算的城鎮居民人均乳制品實際消費支出表示。

2.2.2 關鍵自變量 本文的關鍵自變量為城鎮化水平，由歷年各地區城鎮化率(城鎮常住人口占該地區常住總人口比例)表示。

2.2.3 中介變量 在探究城鎮化對城鎮居民乳制品消費的影響途徑時，運用了以下變量：(1)收入水平：該指標為以1992年不變價計算(經過歷年各地區CPI調整)的城鎮居民人均可支配收入表示。(2)就業結構：由于當地產業結構與就業結構直接聯系，因此就業結構指標由歷年各地區第三產業GDP占地區生產總值比例表示；(3)消費可及性：由歷年各地區的零售餐飲與交通運輸業GDP占地區生產總值的比例表示。從2004年開始，統計資料將批發零售與住宿餐飲業GDP分列，為保持與2004年前統計口徑一致，本文將2004—2011年批發零售業與住宿餐飲業GDP進行了合并；以上統計資料均來自歷年各地區的統計年鑒或年度統計公報。

2.2.4 控制變量 借鑒已有研究[5-6，18]，本文納入的控制變量包括老齡化、價格水平、人力資本、奶業產出、經濟發展狀況以及三聚氰胺事件。選取這些變量的原因是老年人口的增加將提高乳制品的需求從而總體促進居民的乳制品消費；價格水平往往會改變居民的消費預期，價格的上升對人們的乳制品消費不利；人力資本的提高將改善人們對乳制品的認知和消費偏好，從而促進乳制品消費狀況。同時，對于某些地區(比如內蒙古)而言，食物的產量與其飲食習慣有較大的正向相關性，因而奶類的產量可能正向影響居民乳制品消費水平。此外，如三聚氰胺事件這樣的重大食品安全事件也將對影響消費者的乳制品消費有直接的抑制作用。在指標選取上，老齡化由老齡化率(65歲及以上人口占該地區常住人口比例)表示；價格水平由通脹率表示，采用CPI計算；人力資本運用人均受教育程度指標，由6歲及以上人口的人均受教育年限表示，本文將統計資料中顯示的“未上過學”、“小學”、“初中”、“高中”、“大專及以上”等教育程度分別按教育年限為0、6、9、12、16年進行賦值，在此基礎上根據人口抽樣數據測算歷年各地區的人均受教育年限；奶業產出指標為人均奶類產量，指標由包括牛奶、酸奶、奶粉等的產量除以當期該地區常住人口計算而得；經濟發展狀況指標由各地區GDP表示；三聚氰胺事件指標設為虛擬變量，如時期在2008年及以后，則賦值為1，否則賦值為0。以上指標中，通脹率和各地區GDP來自歷年《中國統計年鑒》，人均受教育程度與常住人口數指標來自歷年《中國人口統計年鑒》，奶類產量來自歷年《中國奶業年鑒》。

2.2.5 數據描述 為減小可能存在的異常值或異方差對估計結果的影響，本文對乳制品消費支出、收入水平、奶類產量以及地區GDP均取對數處理。表1可知，各變量均有較大的變異性，這為研究變量間的相互影響的數據變異性提供了基礎。

3 實證結果

3.1 城鎮化與城鎮居民乳制品消費

表2第(1)列為僅加入城鎮化指標與控制變量的回歸結果，第(2)～(5)列是在式(1)的基礎上分別加入各個中介變量和加入全部中介變量的估計結果。在進行具體的實證結果說明與分析前，首先對模型選擇的合理性與可靠性進行檢驗。在工具變量的選取上，參考已有研究，選取因變量的滯后兩期及其更高階作為工具變量進行系統GMM回歸估計，從表2可以看出，各模型的二階自相關檢驗AR(2)檢驗和Sargan檢驗均未通過5%的顯著性檢驗，故接受擾動項無序列相關的原假設以及難以拒絕所有工具變量均有效的原假設。同時，各模型的方差膨脹因子均不大于5，說明不存在明顯的多重共線性問題。以上檢驗結果表明，本研究的模型選擇滿足進行有效估計的前提。

3由于從2013年開始，《中國奶業年鑒》欄目和條目進行了調整，不再詳細匯報全國各地的居民乳制品消費支出數據，因此將研究時期設為1992—2011年。

表1 變量含義與描述性統計

資料來源：根據統計數據計算所得

表2 城鎮化與城鎮居民乳制品消費的系統GMM估計結果

注：括號內為標準誤；*、**、***分別表示在10%、5%和1%的統計水平上顯著

從第(1)列的估計結果來看，城鎮化對城鎮居民乳制品消費的影響系數為正，且通過了1%的顯著性檢驗，說明城鎮化能夠顯著提高城鎮居民的乳制品消費水平。從城鎮化變量系數值看，城鎮化每提高1個百分點，城鎮居民人均乳制品消費將上升1.1個百分點。在控制變量中，老齡化、價格水平、人力資本、奶業產出、經濟發展狀況以及三聚氰胺事件對城鎮居民乳制品消費的影響系數符合均與預期相符。同時，滯后一期的乳制品消費的影響系數顯著為正，說明城鎮居民的乳制品消費有較強的棘輪效應，這也與本文預期一致且證明了模型選擇的合理性。從第(2)～(4)列估計結果可以看出，分別在第(1)列基礎上單獨加入收入水平、就業結構和消費可及性等3個中介變量后，中介變量對乳制品消費的影響系數正向顯著，觀察系數值可知，收入水平、就業結構和消費可及性每增加1個百分點，城鎮居民乳制品消費水平將分別提高0.36、1.52、3.29個百分點。同時城鎮化變量系數值和顯著性均下降，其中在加入收入水平中介變量的模型中城鎮化變量系數降幅最大。說明城鎮化通過收入水平、就業結構和消費可及性影響城鎮居民乳制品消費狀況，其中收入水平發揮更大作用。第(5)列估計結果表明，在自變量中加入多個中介變量后，中介變量顯著正向影響乳制品消費，而城鎮化變量的系數不再顯著，且系數值接近于0，說明城鎮化中介變量幾乎可以解釋城鎮化的影響機制。

3.2 城鎮化影響的異質性

城鎮化對城鎮居民乳制品消費的影響及其機制在不同的城鎮化階段很有可能有所差異，得出這一差異有利于提高促進乳制品消費對策的針對性。本文按照各地區城鎮化率的中位數為界，將城鎮化階段分為高水平城鎮化組和低水平城鎮化組，而后進行城鎮化的影響及其機制的異質性探析。如表4所示，各個系統GMM模型均通過二階序列相關和Sargan檢驗，說明模型設定良好，同時滯后一期乳制品消費估計系數顯著為正，再次表明城鎮居民乳制品消費的時間慣性。觀察城鎮化系數可知，與高水平城鎮化組相比，低水平城鎮化組中的城鎮化影響系數更大，且在統計上更加顯著，說明在城鎮化對城鎮居民乳制品消費的影響在城鎮化初期更大，而后城鎮化率較高時發揮的作用較小，這表明城鎮化的影響存在一定程度的邊際遞減傾向。從各個中介變量結果看，在低水平城鎮化組中，收入水平和就業結構的影響系數均大于高水平城鎮化組，而高水平城鎮化組中的消費可及性影響系數比低水平城鎮化組稍大，這在一定程度上說明低水平城鎮化階段乳制品消費的促進主要依靠收入增長和就業結構改善的作用。根據邊際理論，處于城鎮化率較低階段時，居民收入水平較低，第三產業的就業占比仍較小，提高收入并改變就業結構能夠更加顯著促進居民的消費能力和改變消費習慣，隨著城鎮化進程的推進，收入水平與就業結構的作用就將下降。因此，為促進城鎮居民的乳制品消費，在低水平城鎮化階段應更加側重居民收入水平的提高和就業結構的轉型，零售網絡等消費可及性作用的發揮應以保證居民的消費能力和消費習慣為基礎。

4 結論

本文根據1992—2011年全國29個省(市、區)的面板數據，運用動態面板系統GMM方法實證檢驗了城鎮化對城鎮居民乳制品消費的影響及其機制，同時進一步討論了不同城鎮化水平下城鎮化的影響及其機制的異質性。結果表明：(1)城鎮化對居民乳制品消費有顯著的促進作用，平均而言，城鎮化率每上升1個百分點，城鎮居民人均乳制品消費水平提高1.1個百分點。(2)城鎮化主要通過收入水平、就業結構與消費可及性對城鎮居民乳制品消費產生影響，其中收入水平的作用最為明顯。(3)將城鎮化水平進行分組估計的結果顯示，城鎮化的乳制品消費促進作用存在邊際遞減傾向，即城鎮化的影響在低水平城鎮化階段更加顯著。其原因主要在于收入水平和就業結構的變化對居民消費能力的促進作用和對消費習慣的影響在低水平城鎮化階段更加明顯。

本文的研究結論對理解城鎮化與居民食物消費的關系，以及進一步促進居民乳制品消費和中國奶業發展有一定的現實意義。首先，城鎮化進程對居民的乳制品消費有積極影響，這間接驗證了城鎮化將改變居民的食物消費模式，甚至將改變人們的健康狀態。單從乳制品消費看，城鎮化能在一定程度上改善人們的營養攝入，從而可能有利于居民的健康水平。然而良好的健康狀態需要合理的飲食結構和營養攝取，因此，本文僅對城鎮化與居民營養與健康之間的關系進行了初步探索，后續需要更多研究對這一問題進行探討。其次，從城鎮化對城鎮居民乳制品消費的影響機制來看，要進一步提高乳制品消費水平，應從收入水平、就業結構和消費可及性三方面著手。在城鎮化水平較低地區，更應側重提高居民收入水平以及改善就業結構，提高居民的消費能力和轉變消費習慣，而在城鎮化水平較高的地區更應注重增加乳制品的消費可及性，完善零售網絡，從而更好地促進居民乳制品消費。◇

表4 不同城鎮化水平下城鎮化的影響及其機制差異的系統GMM估計結果

注：括號內為標準誤；*、**、***分別表示在10%、5%和1%的統計水平上顯著。表中的控制變量包括老齡化、價格水平、人力資本、奶業產出、地區經濟發展狀況以及三聚氰胺事件等

[1]鄭志浩, 高穎, 趙殷鈺. 收入增長對城鎮居民食物消費模式的影響[J]. 經濟學 (季刊), 2015,15(1)：263-288.

[2]Kearney J.Food consumption trends and drivers[J].Philosophical Transactions of the Royal Society B Biological Sciences，2010，365(1554)：2793-2807.

[3]蘇暢，張兵，王惠君，等.中國15省(區、直轄市)1991—2015年中老年居民奶類攝入狀況及變化趨勢[J].中國食物與營養，2017，23(12)：76-80.

[4]王莉，田國強.我國奶制品消費行為與態度的調查研究[J].中國乳業，2013(9)：26-30.

[5]馬冠生，崔朝輝，周琴，等.中國居民奶及奶制品消費現狀分析[J].中國食物與營養，2007(11)：36-39.

[6]聶迎利.中國城鎮居民奶類消費特點及影響因素研究[J].農業經濟問題，2008(S1)：165-170.

[7]王貴榮，陳彤，王建軍.關于奶制品消費影響因素的分析[J].中國統計，2009(10)：54-55.

[8]劉媛媛，曾寅初.風險溝通、政府信息信任與消費者購買恢復——基于三聚氰胺事件消費者調查[J].北京社會科學，2014(3)：52-62.

[9]周應恒，卓佳.消費者食品安全風險認知研究——基于三聚氰胺事件下南京消費者的調查[J].農業技術經濟，2010(2)：89-96.

[10]陳天霞.國際視野下的中國奶業：消費特征與未來需求量預測[J].世界農業，2011(4)：1-6.

[11]沈坤榮，余吉祥.農村勞動力流動對中國城鎮居民收入的影響——基于市場化進程中城鄉勞動力分工視角的研究[J].管理世界，2011(3)：58-65.

[12]郭克良，張子麟.城鎮化對于居民收入的提升效應——以河北省為例[J].西北農林科技大學學報(社會科學版)，2015，15(3)：82-87.

[13]喬明睿，錢雪亞，姚先國.勞動力市場分割、戶口與城鄉就業差異[J].中國人口科學，2009，2009(1)：32-41.

[14]葛曉巍，葉俊濤.劉易斯拐點下農民工就業結構及產業結構變化——基于蘇、浙、粵的調查[J].經濟學家，2014(2)：67-72.

[15]黃季焜.社會發展、城市化和食物消費[J].中國社會科學，1999(4)：102-116.

[16]胡定寰，Fuller F，Reardon T.超市的迅速發展對中國奶業的影響[J].中國農村經濟，2004(7)：11-17.

[17]Fuller F H，Beghin J C，Rozelle S.Urban demand for dairy products in China：evidence from new survey data[R].Working Paper，Iowa State University，2004.

[18]Duesenberry J S.Income，saving and the theory of consumer behavior[J].Review of Economics & Statistics，1949，33(3)：111.

[19]Nickell S.Biases in dynamic panel models with fixed effects[J].Econometrica，1981，49(6)：1417-1426.

[20]Che Y, Lu Y, Tao Z, et al. The impact of income on democracy revisited[J]. Journal of Comparative Economics，2013, 41(1):159-169.

[21]溫忠麟，張雷，侯杰泰，等. 中介效應檢驗程序及其應用[J]. 心理學報， 2004,36(5):614-620.