分位數(shù)回歸模型及其應(yīng)用研究

1、2008年中國數(shù)量經(jīng)濟(jì)學(xué)會年會論文 第一組 計量經(jīng)濟(jì)學(xué)理論與方法分位數(shù)回歸模型及其應(yīng)用研究王桂勝 王桂勝：男，1970年生，首都經(jīng)濟(jì)貿(mào)易大學(xué)勞動經(jīng)濟(jì)學(xué)院副教授，清華大學(xué)經(jīng)管學(xué)院博士生。 （首都經(jīng)濟(jì)貿(mào)易大學(xué)，北京，100026）摘要：本文在對分位數(shù)回歸方法的含義和基本原理進(jìn)行全面分析說明的基礎(chǔ)上，對分位數(shù)回歸方法在panel data模型中的應(yīng)用作了深入分析，并對不同回歸估計方法在panel data模型中的估計效果進(jìn)行了比較分析。在此基礎(chǔ)上，通過分別采取一般最小平方法和分位數(shù)回歸法對中國15省區(qū)的人均消費和人均收入的回歸方程估計的統(tǒng)計結(jié)果比較，發(fā)現(xiàn)分位數(shù)回歸方法在進(jìn)行某些特殊的panel da

2、ta模型估計時具有一定的優(yōu)勢。關(guān)鍵詞：分位數(shù)回歸、面板數(shù)據(jù)模型、懲罰分位數(shù)回歸估計一、分位數(shù)回歸研究介紹 自koenker 和 bassett (1978)提出線性分位數(shù)回歸理論以來，分位數(shù)回歸（qr）即成為近幾十年來發(fā)展較快、應(yīng)用廣泛的回歸模型方法，它不僅深化了對傳統(tǒng)回歸模型的理解，而且也推廣了回歸模型的類型和應(yīng)用，使得回歸模型擬合有關(guān)統(tǒng)計數(shù)據(jù)更加準(zhǔn)確細(xì)致。分位數(shù)回歸模型是在穩(wěn)健估計模型基礎(chǔ)上發(fā)展形成。穩(wěn)健估計（robust estimation）理論包括基于一般凸損失函數(shù)的m 估計理論、基于樣本秩統(tǒng)計量的r估計理論和基于樣本次序統(tǒng)計量的l估計理論等。分位數(shù)回歸強調(diào)以解釋變量的分位數(shù)來估計推

3、斷因變量的分位數(shù)，通過建立分位數(shù)估計方程，并運用線性規(guī)劃方法或非參數(shù)估計等方法來估計相應(yīng)于不同分位數(shù)的解釋變量系數(shù)或未知參數(shù)。分位數(shù)回歸是中位數(shù)回歸和均值回歸的推廣。分位數(shù)回歸模型具體又分為四分位數(shù)回歸、十分位數(shù)回歸、百分位數(shù)回歸、logit分位數(shù)回歸、審查分位數(shù)回歸等模型。關(guān)于分位數(shù)回歸研究的最近發(fā)展，主要表現(xiàn)在分位數(shù)回歸技術(shù)方法和方法應(yīng)用等兩方面的研究上。具體包括koenker和zhijie xiao (2000) 解決分位數(shù)回歸過程中存在的特定推斷問題; kim 和muller (2000) 關(guān)于雙步分位數(shù)回歸的漸進(jìn)特性的研究; tasche (2001) 對最小分位數(shù)回歸的無偏性研究;

4、 chernozhukov 和han hong (2002) 提出對審查分位數(shù)回歸的三步評估法;吳建南、bret-schneider 等(2002) 用蒙特卡羅(monte carlo) 方法產(chǎn)生100 個隨機數(shù)據(jù)集合來比較顯著權(quán)重分析方法與分位數(shù)回歸的優(yōu)劣; kottas 和krnjajic 提出分位數(shù)回歸中的貝葉斯非參數(shù)模型;koenker(2004)將分位數(shù)回歸方法運用于panel data模型估計中，并提出了pqr估計技術(shù)及相關(guān)理論證明，在此基礎(chǔ)上carlos lamarche（2006）對pqr估計方法進(jìn)一步深入探討并結(jié)合實際數(shù)據(jù)進(jìn)行實證分析等等。此外是對分位數(shù)回歸方法的應(yīng)用研究。在

5、這一方面,barnes 和w. hughes (2002) 利用分位數(shù)回歸對跨部門公債市場的回收進(jìn)行了分析;buhai (2004) 在分析介紹了分位數(shù)回歸方法的基礎(chǔ)上,研究了它在持續(xù)時期模型和循環(huán)結(jié)構(gòu)等式模型中的應(yīng)用;leggett 和craighead 利用分位數(shù)回歸確定了時間分布和特定風(fēng)險驅(qū)動的影響。國內(nèi)也有很多學(xué)者將分位數(shù)回歸估計方法運用于醫(yī)學(xué)衛(wèi)生事業(yè)研究、公共管理事業(yè)等具有極值分布特點的統(tǒng)計數(shù)據(jù)研究中。具體如教育收益率估計分析、勞工歧視因素分析等等。本文將在系統(tǒng)介紹分位數(shù)回歸模型的含義和基本原理的基礎(chǔ)上，重點分析將分位數(shù)回歸方法與panel data模型結(jié)合分析的主要形式和估計原理以

6、及相關(guān)結(jié)果。第二節(jié)討論分位數(shù)回歸方法的一般原理以及pqr估計原理，在此基礎(chǔ)上提出關(guān)于panel data的分位數(shù)回歸模型的pqr估計方法；第三節(jié)分析分位數(shù)回歸方法在panel data模型中的應(yīng)用以及數(shù)據(jù)模擬；第四節(jié)對本文的分析結(jié)果進(jìn)行總結(jié)。二、分位數(shù)回歸模型原理（一）分位數(shù)回歸的基本原理一般線性回歸模型可設(shè)定如下：在滿足高斯-馬爾可夫假設(shè)前提下，可表示如下：其中，為待估解釋變量系數(shù)。上述模型即均值回歸模型表達(dá)式，是對等式兩邊取數(shù)學(xué)期望的結(jié)果。類似于均值回歸模型，也可設(shè)定中位數(shù)回歸模型如下：其中，為關(guān)于x的條件中位數(shù)，為隨機擾動項的中位數(shù)。以及分位數(shù)回歸模型如下：對于均值回歸模型，可采取最小平

7、方法（ols）估計未知參數(shù)；對于中位數(shù)回歸模型，可采取最小一乘法（或稱最小絕對偏差法lad）；而對于分位數(shù)回歸模型，則可采取線性規(guī)劃法（lp）估計其最小加權(quán)絕對偏差，從而得到解釋變量的回歸系數(shù)。分別可表示如下：ols法： 求解得：lad法： 求解得：qr法： 求解得：其中，。（二）panel data模型與分位數(shù)回歸方法考慮一般panel data模型，表達(dá)式如下： 其中，i代表不同的樣本個體，t代表不同的樣本觀察時點，u表示隨機誤差項，表示解釋變量的系數(shù)向量，i表示第i個樣本的不可觀察隨機效應(yīng)。關(guān)于的估計有兩種情形，即固定效應(yīng)情形和隨機效應(yīng)情形，相應(yīng)的有兩種估計方法。將上述panel dat

8、a模型寫成如下矩陣形式：固定效應(yīng)情形下，的估計量為 隨機效應(yīng)情形下，的估計量為 假設(shè)u n(0,r), n(0,w)，v=z+u。則有 ?？梢詫Υ饲樾蝡anel data模型采取gls估計法、也可采取反映個體影響的懲罰最小二乘法（pls）估計回歸系數(shù)，分別表述如下：gls法： pls法：二者的共同解為：。對于上述panel data模型也可采取分位數(shù)回歸法進(jìn)行參數(shù)估計。為此。建立以下條件分位數(shù)方程：上述分位數(shù)方程假設(shè)個體效應(yīng)是固定的。對此方程，koenker（2004）提出了懲罰分位數(shù)回歸（pqr）方法進(jìn)行估計。具體如下：其中，wj為相應(yīng)于各分位數(shù)的權(quán)數(shù)，為調(diào)節(jié)系數(shù)。如果=0，則為固定效應(yīng)分位

9、數(shù)回歸估計量（feqr）；如果>0，則為懲罰分位數(shù)回歸估計量（pqr）。此外，koenker還對分位數(shù)回歸和懲罰分位數(shù)回歸估計量的漸進(jìn)特性進(jìn)行了詳細(xì)的考察。在此基礎(chǔ)上，運用蒙特卡羅模擬法在小樣本設(shè)定下比較分析了不同回歸估計方法的效果，如下表所示： 表 1 最小平方法與分位數(shù)回歸法的估計效果比較誤差項分布 ls pls lsfe qr pqr qrfe正態(tài)分布偏差 0.0031 0.0048 0.0056 0.0048 0.0067 0.0047均方誤差平方根 0.0847 0.0604 0.0668 0.0977 0.0781 0.0815t分布偏差 -0.0062 -0.0054 -0

10、.0051 -0.0063 -0.0101 -0.0082均方誤差平方根 0.1377 0.1031 0.1143 0.1274 0.0881 0.0921卡方分布偏差 -0.0068 0.0002 0.0032 -0.0052 0.0063 0.0072均方誤差平方根 0.2155 0.1503 0.1650 0.2362 0.1506 0.1513由上表可見，對于某些非正態(tài)分布，分位數(shù)回歸估計效果相對較好一些（如表1中的t分布）。三、分位數(shù)回歸方法的應(yīng)用分位數(shù)回歸方法對于具有某些非高斯分布的隨機誤差項和隨機效應(yīng)的回歸方程的系數(shù)估計具有較好的效果，這已在理論上得到了論證。本文將運用實際經(jīng)濟(jì)數(shù)

11、據(jù)對比最小平方法與分位數(shù)回歸方法在panel data模型中的估計效果。所用數(shù)據(jù)來自中國1996-2002年東北、華北和華東15個省、市、自治區(qū)的居民人均消費（cp，不變價格）和人均收入（ip，不變價格） 注：本數(shù)據(jù)資料來自張曉峒著eviews使用指南與案例，機械工業(yè)出版社，2007年，第261頁。建立關(guān)于人均消費與人均收入關(guān)系的panel data模型，分別采取最小平方法和分位數(shù)回歸法進(jìn)行估計，由于樣本資料原因，本文僅考慮合并數(shù)據(jù)下和固定效應(yīng)下的兩種估計方法的估計結(jié)果，并進(jìn)行比較分析。1、 合并數(shù)據(jù)下最小平方法估計運用eviews軟件對樣本數(shù)據(jù)進(jìn)行處理，結(jié)果如下：其中，判定系數(shù)為0.98,常

12、數(shù)項t統(tǒng)計量值為2.0,回歸系數(shù)t統(tǒng)計量值為79.7。2、固定效應(yīng)下最小平方法估計對15個省區(qū)設(shè)計15個虛擬變量d1，d2，d15。其定義為： 運用eviews軟件，對樣本數(shù)據(jù)進(jìn)行處理，得到如下結(jié)果： 其中，判定系數(shù)為0.99,常數(shù)項t統(tǒng)計量值為1.6,回歸系數(shù)t統(tǒng)計量值為55。3、合并數(shù)據(jù)下分位數(shù)回歸法估計將15個省區(qū)的數(shù)據(jù)當(dāng)作一個個體的數(shù)據(jù)進(jìn)行分位數(shù)回歸估計，運用stata軟件運行后，得到如下結(jié)果：表 2 各分位數(shù)下合并數(shù)據(jù)分位數(shù)回歸估計結(jié)果分位數(shù)ip常數(shù)項判定系數(shù)（r2）系數(shù)t統(tǒng)計量系數(shù)t統(tǒng)計量0.20.724648.63187.25272.730.84440.50.742258.852

13、05.25913.090.87430.750.779229.60139.57041.140.89300.80.798931.3275.44620.620.89684、固定效應(yīng)下分位數(shù)回歸法估計與2中固定效應(yīng)回歸一樣，可設(shè)計15個虛擬變量d1，d2，d15，并運用stata軟件進(jìn)行分析可得如下結(jié)果：表 3 各分位數(shù)下固定效應(yīng)分位數(shù)回歸估計結(jié)果分位數(shù)ip常數(shù)項判定系數(shù)（r2）系數(shù)t統(tǒng)計量系數(shù)t統(tǒng)計量0.20.679936.28490.54115.560.90960.50.671519.91597.55093.270.92120.750.732825.23351.14422.220.93290.80

14、.740824.76309.19151.900.9359比較分析上述兩種回歸方法的統(tǒng)計結(jié)果，發(fā)現(xiàn)在固定效應(yīng)情形下兩種方法回歸效果均比合并數(shù)據(jù)情形下更好；在同一情形下做回歸分析，顯然分位數(shù)回歸分析結(jié)果更加穩(wěn)定，各系數(shù)估計顯著程度更高。因而，分位數(shù)回歸估計在panel data模型中可以發(fā)揮重要作用。四、結(jié)論本文在對分位數(shù)回歸方法的含義和基本原理進(jìn)行全面分析說明的基礎(chǔ)上，對分位數(shù)回歸方法在panel data模型中的應(yīng)用作了深入分析，并對不同回歸估計方法在panel data模型中的估計效果進(jìn)行了比較分析。一般而言，分位數(shù)估計方法在估計具有非正態(tài)分布的誤差項或不可觀察的隨機效應(yīng)時具有一定優(yōu)勢。本文

15、在理論分析之后，提供了一個應(yīng)用案例分析，通過對我國人均收入和人均消費的各種回歸分析，充分證明了分位數(shù)回歸的較好效果。當(dāng)然，由于樣本數(shù)據(jù)的不足，缺乏對在隨機效應(yīng)情形下兩種回歸方法估計的效果比較，以及在不同調(diào)節(jié)系數(shù)下對懲罰分位數(shù)回歸估計效果的考察。這將在分位數(shù)回歸方法的應(yīng)用分析中進(jìn)一步開展深入探討。參考文獻(xiàn)1emanuel parzen, quantile probability and statistical data modeling, j.statistical science, vol. 19, no. 4, (nov., 2004), pp. 652-662.2roger koenker

16、 and zhijie xiao, inference on the quantile regression process, j.econometrica, vol. 70, no. 4, (jul., 2002), pp. 1583-1612.3keming yu, zudi lu, julian stander, quantile regression: applications and current research areas,j. the statistician, vol. 52, no. 3, (2003), pp. 331-350.4. roger koenker and

17、jose a. f. machado, goodness of fit and related inference processes for quantile regression,j.journal of the american statistical association, vol. 94, no. 448, (dec., 1999), pp. 1296 1310.5. keming yu and m. c. jones, local linear quantile regression, j.journal of the american statistical associati

18、on, vol. 93, no. 441, (mar., 1998), pp. 228 237.6. james w. taylor and derek w. bunn, a quantile regression approach to generating prediction intervals,j .management science, vol. 45, no. 2, (feb., 1999), pp. 225-237.7. rene gelinas, alain martel, pierre lefrancois, sos: a quantile estimation proced

19、ure for dynamic lot-sizing problems,j. the journal of the operational research society, vol. 46, no. 11, (nov., 1995), pp. 1337-1351.8. p. k. bhattacharya and ashis k. gangopadhyay, kernel and nearest-neighbor estimation of a conditional quantile,j. the annals of statistics, vol. 18, no. 3, (sep., 1

20、990), pp. 1400-1415.9. portnoy s. koenker r., the gaussian hare and the laplacian tortoise :computability of squared - error versus absolute- error estimatorsj . statistical science ,1997(12) :279 - 300.10. koenker , r. and bassett , g. 1978.“regression quantiles. ”j . econometrica. vol. 46 , pp. 33

21、 - 50.11. koenker, r. (2004): “quantile regression for longitudinal data,” j,journal of multivariate analysis, 91,7489.12.吳建南、馬偉，估計極端行為模型:分位數(shù)回歸方法及其實現(xiàn)與應(yīng)用，j數(shù)理統(tǒng)計與管理，2006、9月，第25卷第5期，536-543。 research on quantile regression modeling and its applicationwang guisheng(capital university of economics and business,beijing,100026)abstract: this article principally introd