《张晓峒分位数回归讲义_经济学_高等教育_教育专区.docx》由会员分享，可在线阅读，更多相关《张晓峒分位数回归讲义_经济学_高等教育_教育专区.docx（27页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、1第第 15 章章 分位数回归模型分位数回归模型15.1 总体分位数和总体中位数15.2 总体中位数的估计15.3 分位数回归15.4 分位数回归模型的估计15.5 分位数回归模型的检验15.6 分位数的计算与分位数回归的 EViews 操作15.7 分位数回归的案例分析以往介绍的回归模型实际上是研究被解释变量的条件期望。人们当然也关心解释变量与被解释变量分布的中位数，分位数呈何种关系。这就是分位数回归，它最早由 Koenker 和 Bassett(1978)提出，是估计一组回归变量 X 与被解释变量 Y 的分位数之间线性关系的建模方法。正如普通最小二乘 OLS 回归估计量的计算是基于最小化残
2、差平方和一样，分位数回归估计量的计算也是基于一种非对称形式的绝对值残差最小化，其中，中位数回归运用的是最小绝对值离差估计(LAD，least absolute deviations estimator)。它和 OLS 主要区别在于回归系数的估计方法和其渐近分布的估计。在残差检验、回归系数检验、模型设定、预测等方面则基本相同。分位数回归的优点是，（1）能够更加全面的描述被解释变量条件分布的全貌，而不是仅仅分析被解释变量的条件期望（均值），也可以分析解释变量如何影响被解释变量的中位数、分位数等。不同分位数下的回归系数估计量常常不同，即解释变量对不同水平被解释变量的影响不同。另外，中位数回归的估计方
3、法与最小二乘法相比，估计结果对离群值则表现的更加稳健，而且，分位数回归对误差项并不要求很强的假设条件，因此对于非正态分布而言，分位数回归系数估计量则更加稳健。15.1 总体总体分位数分位数和和总体总体中中位数位数在介绍分位数回归之前先介绍分位数和中位数概念。对于一个连续随机变量 y，其总体第分位数是 y()的定义是：y 小于等于 y()的概率是，即=P(y y()=F(y()其中 P()表示概率，F(y()表示 y 的累积（概率）分布函数(cdf)。比如 y(0.25)=3，则意味着 y 3 的概率是 0.25。且有y()=F-1(y()即 F(y()的反函数是 y()。当=0.5 时，y()
4、是 y 的中位数。=0.75 时，y()是 y 的第 3/4 分位数，=0.25 时，y()是 y 的第 1/4 分位数。若 y 服从标准正态分布，y(0.5)=0，y(0.95)=1.645，y(0.975)=1.960。另外，如果随机变量 y 的分布是对称的，那么其均值与中位数是相同的。当其中位数小于均值时，分布是右偏的。反之，分布是左偏的。对于回归模型，被解释变量 yt对以 X 为条件的第分位数用函数 y()tX 表示，其含义是：以 X为条件的 yt小于等于 y()tX 的概率是。这里的概率是用 yt对 X 的条件分布计算的。且有y()tX=F-1(y()tX)其中 F(y()tX)是
5、yt在给定 X 条件下的累积概率分布函数(cdf)。则 y()tX 称作被解释变量 yt对 X的条件分位数函数。而 F(y()tX)=f(y()tX)则称作分位数概率密度函数。其中F(y()tX)表示 F(y()tX)2对 y()tX 求导。15.2 总体总体中中位数的估计位数的估计在介绍分位数回归之前，先来看中位数的估计和中位数回归。下面以连续变量为例介绍定理15.1。定理 15.1连续变量用 y 表示，其概率密度函数用 f(y)表示，累计概率密度函数用 F(y)表示，y 的中位数用 y(0.5)表示，则 y 与任一值的离差绝对值的期望)(yE以=y(0.5)时为最小。证明：)(yE=dyy
6、fydyyfy)()()()(-=)()()()(-ydFyydFy(15.1)根据莱布尼兹公式，若dyyfFba),()(，则有dyyfFba),()(。令-),(yyf，则有babadydyyF-)-()(。运用于式（15.1），得)(tyE=dyyfydyyfy)()()()(-=)(-)(-ydFydF=1-)(2)(-(1-)()(-1-)(-FFFydFF式（15.1）求极小的一阶条件是)(tyE=0，即1-)(2F=0，0.5)(F。这意味着等于中位数 y(0.5)。=y(0.5)与定理 15.1 等价的表述是y以=y(0.5)（中位数）时为最小。因此，中位数回归估计量可以通过最
7、小绝对离差法（least absolute deviation,LAD）估计。其中 X 和 分别为(k1)阶列向量。同理，对于线性回归模型 yt=X +ut，通过求(0.5)Xty最小，估计 的中位数回归系数估计量(0.5)，从而得到 yt的中位数回归估计量(0.5)5.0()(XX ty。15.3 分位数回归分位数回归Koenker 和 Bassett(1978)证明，若用ty)(表示 yt的分位数回归估计量，则对于以检查函数（check function）w为权数，yt对任意值的加权离差绝对值和tyw只有在=ty)(时取得最小值。其中tyw=)()(1(:TyttTyitiiyy(15.2
8、)3(0,1)。据此，分位数回归可以通过加权的最小绝对离差和法（weighted least absolute deviation,WLAD）进行估计。根据式(15.2)，对于线性回归模型 yt=X +ut，求第分位数回归方程系数的估计量)(的方法是求下式（目标函数）最小，TutTutttuuQ0)(0)()()()1(TXyttTXyttttyy)()(:)(:)()()(1(XX(15.3)其中tu)(表示第分位数回归方程对应的残差。(0,1)。第分位数的回归方程表达式是ty)(=)(X其中 X，都是 k1 阶列向量。)(称作分位数回归系数估计量，或最小绝对离差和估计量，估计方法称作最小绝
9、对离差和估计法。当=0.5 时，式（15.3）变为TttTXyttTXyttyyyQtt1)0.5(:)0.5(:)0.5(0.5)(0.5)(0.5)0.5()0.5(XXXty)0.5(=)0.5(X称作中位数回归方程，)0.5(称作中位数回归系数估计量。一旦得到估计的分位数回归方程，就可以计算分位数回归的残差tu)(。ttttyyyu)()(-)(X对一个样本，估计的分位数回归式越多，对被解释变量 yt条件分布的理解就越充分。以一元回归为例，如果用 LAD 法估计的中位数回归直线与用 OLS 法估计的均值回归直线有显著差别，则表明被解释变量 yt的分布是非对称的。如果散点图上侧分位数回归
10、直线之间与下侧分位数回归直线之间相比，相互比较接近，则说明被解释变量 yt的分布是左偏倚的。反之是右偏倚的。对于不同分位数回归函数如果回归系数的差异很大，说明在不同分位数上解释变量对被解释变量的影响是不同的。15.4 分位数回归模型的估计分位数回归模型的估计由于目标函数（15.3）不可微，因此传统的对目标函数求导的方法不再适用。估计分位数回归方程参数)(的一种较好的方法是线性规划方法。基于 Barrodale 和 Roberts(1973，以下简写为 BR)提出的单纯形法(simplex algorithm)，Koenker和 DOrey(1987)提出一种估计分位数回归系数的方法。EView
11、s 中应用的是上述算法的改进形式。BR算法由于其非有效性和大样本下的一些非优良特性曾备受批评。Koenker和Hallock(2001)以及 Portnoy 和 Koenker(1997)通过模拟证实，与内点法(interior point method)等替代方法相比，BR 算法的估计次数往往较多，大约是样本容量的平方次数。然而，改进的 BR 算法的估计次数在一定程度上是可以接受的，大约是样本容量的线性倍次数，在实际中是可以使用的。4分位数回归方程的 BR 算法原理略。下面讨论分位数回归系数估计量的渐近分布。在弱条件下，分位数回归系数渐近服从正态分布(Koenker,2005)。回归系数的方
12、差协方差矩阵的计算在分位数回归的系数估计中占有重要位置。其方差协方差矩阵的估计方法根据分位数密度函数是否与解释变量相关分为三种方法：误差项独立同分布(i.i.d.)假设下的直接估计方法。由 Koenker 和 Bassett(1978)提出。误差项独立但不同分布(i.n.i.d.)条件下的直接估计方法。误差项独立同分布(i.i.d.)和独立但不同分布(i.n.i.d.)条件下都可使用的自举法。（1）独立同分布假设下的参数渐近分布独立同分布假设下的参数渐近分布Koenker 和 Bassett(1978)在独立同分布假设下得出分位数回归系数渐近服从正态分布，可以表述为在弱条件下：)()()(n)
13、1(,0(12)(JsN(15.5)其中)(lim)(limTXXTXXJniiin(15.6)(/1)(11)(FfFs(15.7)其中 s()称为稀疏函数(Sparsity function)或分位数密度函数(quantile density function)。s()是分位数函数的导数，或在第分位数条件下概率密度函数的倒数(见 Welsh,1988)。另外，模型误差项独立同分布假设意味着 s()与解释变量 X 无关，因此，分位数方程只和 X 在局部期间相关，即所有的条件分位数平面互相平行。事实上，式(15.5)中的)1(12)(Js就是误差项独立同分布假设下解释变量的回归系数估计量的渐近
14、方差协方差矩阵表达式，而2)()1(s代表的是一般回归方程中随机误差项的方差。误差项独立同分布假设下，分位数回归参数估计量的渐近方差协方差矩阵表达式中含有s()，但 s()是未知分布的函数，而且必须要估计。EViews 提供了三种估计 s()的方法。两种是基于 Siddiqui(1960)的方法分别提出的差分商方法(Siddiqui Difference Quotient)（Koenker(1994)以及Bassett和Koenker(1982)），一种是核密度（KernelDensity）估计法。简述如下：Siddiqui 差分商法：差分商方法是用实际的分位数函数构造一个简单的差分商，从而求
15、得 s()的估计量，表达式如下：nnnhhFhFs2)()()(11(15.8)其中带宽 hn随着样本容量 n而趋向于 0。要计算()需要做两件事，一是得到分位数函数)(1F在两个点上的值，二是确定带宽。EViews 中提供了两种 Siddiqui 差分商法。计算分位数密度函数的第一种方法由 Bassett 和 Koenker(1982)提出，EViews 将其称之为Siddiqui(mean fitted)方法。这种方法需要重新估计两个分位数回归模型在-hn和+hn上的拟和值，进而用不同的估计参数计算分位数函数的拟和值。最终 s()的估计量的数学表达式如下，对5任意 X*有：nnnhhhXs
16、2)()(*)(15.9)独立同分布假设意味着 X*可以取任何值，Bassett 和 Koenker 建议取 X 的均值，其优点是：估计的精度在该点达到最大；且估计的分位数函数对是单调的，因此对一个恰当的 hn，()的值总是正的。另一种 Siddiqui 差分商法由 Koenker(1994)提出。其计算量相对较小，只需计算原分位数回归方程中残差的第-hn和+hn实际分位数，计算时排除在估计中设为零的 k 个残差，并插入新值以获得分位数的分段线性形式。EViews 中把这种方法叫做 Siddiqui(residual)方法。上述两种 Siddiqui 方法都需要估计带宽 hn。EViews 提
17、供了三种估计带宽的方法：Bofinger(1975)法，Hall-Sheather(1988)法和 Chamberlain(1994)方法。Bofinger(1975)提出的估计带宽的表达式为：5/1221415/1 1)(2)(5.4Thn(15.10)可以近似最小化()的均方误差(MSE)。另外两个带宽的表达式中含有显著性水平，因此常常用来进行假设检验。其中Hall和Sheather(1988)的表达式为：3/121213/23/11)(2)(5.1ZThn(15.11)其中T表示样本容量，表示正态分布的积累分布函数，表示正态分布的密度函数，Z=-1(1-/2)为选择的显著性水平对应的Z值
18、。Chamberlain(1994)的表达式为：TZhn)1(15.12)图1是样本容量1300时Hall和Sheather(1988)方法在第0.1、0.3、0.5、0.7、0.9分位数下得到的带宽。图2是样本容量11000时三种方法在第0.5分位数下的带宽比较图(=0.05，MATLAB计算）。图1图2从图2可以看出随着样本的增加，三种带宽都减小，并且在小样本时，减小的速度较大，在大样本情况下减小的速度较小。并且在大样本情况下，带宽的大小顺序为：Bofinger的最大，Hall和Sheather的次之，Chamberlain的最小。核密度法(Kernel Density)：6根据(15.7
19、)式有 s()=)(1F=1/f()(1F)，Falk(1988)和 Welsh(1988)提出了用核密度法估计)(1F进而得到 s()的方法。而 Powell(1986)、Jones(1992)以及 Buchinsky(1995)则通过估计1/f()(1F)来得到 s()。EViews 中使用的方法属于后者，沿用了 Powell(1984,1989)中的计算方法，其选项名称为 Kernel(residual)：TiTiTcuKcTs1)(1)/()/1/(1)(15.13)其中()表示分位数回归的残差；cT为带宽；K 表示核密度函数。EViews 中可以选择的核密度函数有 Epanechni
20、kov 核函数、均匀(Uniform)核函数、三角(Triangular)核函数、二权(Biweight)核函数、三权(Triweight)核函数、正态(Normal)核函数、余弦(Cosinus)核函数。EViews 中使用了 Koenker(2005)提出的带宽，表达式为：)()(11TTThhkc(15.14)其中 k 表示 Silverman(1986)的一个稳健估计量；hn是 Siddiqui 带宽。(2)独立但不同分布假设下的参数渐近分布当分位数密度函数独立但不同分布即与解释变量 X 相关时)()(T的渐近分布服从Huber sandwich 形式：)()()(T)()()1(,0
21、(11JHHN(15.15)其中 J 同(15.6)式，H 的表达式如下：)/)(lim)(TqfXXHiiiiiT(15.16)其中)(iiqf是个体 i 在第分位数上的条件密度函数。如果条件密度函数不依赖于观测值，式（15.15）中的方差就退化为(15.5)式中的方差。对于 H，EViews 提供了两种计算方法。第一种是 Hendricks 和 Koenker(1992)提出的 Siddiqui差分法；另一种是 Powell(1984,1989)提出的核密度法。这两种方法与在独立同分布假设时计算s()的算法相同，因此在 EViews 选单中的名称相同，分别为 Siddiqui(mean f
22、itted)和 Kernel(residual)。Siddiqui 差分商法这种方法需要对每个个体估计-hn和+hn两个分位数回归模型，将拟和值代入下式：)()(/(2 )()(/(2)(11TTiTTiiTiiTiihhXhhqFhqFhqf(15.17)由于分位数密度函数非同分布，因此，我们需要为每一个个体估计)(iiqf，这时当取XXi时，不能保证(15.17)式为正，因此，Hendricks 和 Koenker 对其进行了修正：)()(/(2,0max()(TTiTiihhXhqf(15.18)其中是一个很小的正数，避免上式中分母为零。将(15.18)式代入(15.16)式，得到 H
23、的估计量为TXXqfHiiiii/)()(15.19)7核密度法Powell(1984,1989)提出的用核密度法估计 H 的表达式为：iiTiTiTXXcuKcTH1)(1)/()/1()(15.20)其中()表示分位数回归的残差；cn为带宽；K 表示核密度函数；各参数含义与(15.13)式相同。（3）参数渐近分布的自举法前面的方法都是先求出分位数密度函数，然后再得到参数的渐近分布。自举法则可以省略这一步，直接得到参数的方差协方差阵。EViews 中给出了四种自举方法，分别为：残差自举，XY对自举，以及两种马尔可夫链边际自举法 MCMB 和 MBMB-A。其中前两种方法见 Buchinsky
24、(1995)。残差自举法(residual bootstrap)这种方法要求解释变量与随机误差项不相关。它是对残差和解释变量分别进行有放回的再抽样，构造样本容量为 m 的新序列 u*和 X*(其中 m 可以小于原样本容量 T)，然后运用初始参数估计量构造被解释变量，即*)(*uXy，最后用 X*和 Y*估计新的参数().如此重复 K 次，则参数方差协方差阵的估计量为：)()()()(1)()(1jBjjBTmTV(15.21)其中)(是自举参数估计量序列的均值。EViews 选单中称这种方法为 Residual。XY 对自举法(XY-pair or design bootstrap)这是最常用
25、的一种自举方法，它不要求随机误差项与解释变量相互独立。使用这种方法时，我们从原始数据中有放回的抽取 K 次样本容量为 m 的子序列(y*,X*)，然后用每个子序列计算()，最后运用(15.21)式计算参数方差协方差阵的估计量。EViews 选单中这种方法称为 XY-Pair。马尔可夫链边际自举法(Markov Chain Marginal Bootstrap)以上两种自举法往往计算量过大，当方程中含有 p 个参数时，每次自举都需要解一个 p 维的线性规划问题。He 和 Hu(2002)提出了一种新的自举法，将一个 p 维的最优问题简化为求解一个含 p 个元素的序列的一维问题。这个序列的一维解就
26、构成了一个马尔可夫链，其样本方差协方差阵可由(15.21)式计算，且当原序列样本容量 T 和自举次数 K 较大时具有一致性。EViews 选单中把这种方法称为 MCMB 法。然而，给定链长 B(即自举次数)，上述方法计算的参数序列之间往往存在较强的自相关从而导致参数方差协方差阵估计量的统计特性较差，有可能对任何链长 B，估计量都不能收敛。Kocherginsky、He 和 Mu(KHM,2005)提出了一种修正的方法消除可能存在的自相关。即通过先对参数空间进行某种转换，运用 MCMB 算法进行估计，然后再转换回原来的空间，这种方法叫做MCMB-A。它要求独立同分布的假设条件，但它对异方差的情况
27、表现的比较稳健。Kocherginsky、He 和 Mu 还建议对于满足 Tmin(,1-)5p 的情况，当 T 1000，p 10 时，B 应取在 100 至 200之间。对于 Tp 在 10,000 到 2,000,000 之间的情况，建议 B 取在 50 至 200 之间，当然，还取决于用户的耐心。15.5 分位数回归模型的分位数回归模型的检验检验评价分位数回归函数好坏的统计量主要有 3 个，拟合优度、拟似然比检验和 Wald 检验。（1）拟拟合合优度优度(Goodness-of-Fit)Koenker 和 Machado(1999)提出了分位数回归的拟合优度的概念。它与一般回归分析中的
28、 R28很类似。假设分位数回归直线为)()(Xy将解释变量矩阵和参数向量都分为两部分，即),1(ZX和),()(1)(0)(，且有)(1)(0)(Zy定义：)()(1(min)()(:)(1)(0:)(1)(0)(TXyttTXyttttZyZyQ(15.22)()(1(min)(:)(0)(:)(0)(TXyttTXyttttyyQ(15.23)式(15.22)和(15.23)分别表示无约束分位数回归目标函数（最小绝对离差和）和约束的分位数回归目标函数（最小绝对离差和）的极小值。无约束目标函数中的减项既包含常数项也包含所有回归因子。约束目标函数中的减项仅包含常数项，其他参数都约束为零。则 K
29、oenker 和 Machado 拟和优度准则表达式如下：)()()(*1QQR(15.24)很明显，上述统计量与传统的 R2非常相似。因为)()(QQ，所以 R*()的值在 0 和 1 之间，解释变量的作用越强，)(Q越远远小于)(Q，)(*R越接近 1。反之越接近 0。所以)(*R可用来考察解释变量对被解释变量第分位数回归拟和的好坏。（2）拟似然比检验拟似然比检验(Quasi-Likelihood Ratio Tests)Koenker 和 Machado(1999)根据目标函数在施加约束条件前后得到的两个极小值构造了两个拟似然比检验统计量（QLR）。这两个拟似然比检验也称作分位数检验(q
30、uantile-tests)。两统计量的表达式如下：)()1()(2)()()(sQQLT(15.25)log()()1(2)()()()(QQsQT(15.26)两个统计量都渐近服从自由度为 q 的2分布，其中 q 是原假设目标函数中约束条件的个数。)(Q和)(Q分别代表约束的和无约束目标方程的极小值。另外，两个统计量的分母都含有稀疏项 s()，上面给出的稀疏项 s()的 3 种计算方法都可在式（15.25）和（15.26）中使用。EViews 估计的是其在备择假设下的估计量。使用上述两统计量的前提是必须满足分位数密度函数 s()与解释变量 X 不相关。然而，尽管有时并不满足独立同分布的假设
31、，EViews 在进行分位数回归的时候，不管选择何种估计参数渐9近分布的方法，总会估计稀疏函数 s()，从而构造拟似然比（QLR）检验统计量。因此，这种检验方法与下面的 Wald 统计量相比稳健性较差。（3）Wald 检验检验给定分位数回归参数估计量的渐近方差协方差矩阵，我们就可以构造 Wald 形式的统计量进行各种约束形式的参数检验。31.3.5 系列系列分位数分位数回归回归检验检验前面的分析主要集中在单个分位数回归模型的假设检验上，而有些时候也需要对一系列分位数回归的回归系数进行联合检验。比如，需要通过检验不同分位数模型的斜率是否相等来判断一个模型是否具有位移特征。同时考虑多个分位数回归式
32、称作系列分位数回归分析（quantile processtesting）。EViews 在做单方程分位数回归的同时，有专门命令执行系列分位数回归分析。操作路径是在一个分位数回归估计结果窗口，点击 View 键，选 Quantile Process/ProcessCoefficients 功能。定义系列分位数回归系数列向量为，),()()()(21m(15.27)则有)(n),0(N(15.28)其中由如下形式的块矩阵ij(kmkm)组成：)()(),min(11jijijiijJHH(15.29)i,j=1,2,m.k 为方程待估参数个数。其中 J 的表达式见(4)式。H 的表达式见(15.1
33、9)或(15.20)式，取决于选择的估计方法。特别的，当误差项独立同分布的假设成立时，简化为：J0(15.30)其中0中的元素如下：)()(),min(11jijijiijFfFf(15.31)i,j=1,2,k.除了以上的方法以外，的估计量还可以由任何一种自举方法得到。（1）斜率相等检验Koenker 和 Bassett(1982a)提出了一种对异方差很稳健的判断不同分位数回归方程斜率是否相等的检验。零假设如下：)(1)(1)(1021:mH 其中1指常数项以外的解释变量所对应的(k-1)维参数列向量。因此，零假设共含有(k-1)(m-1)个约束条件。接下来构造 Wald 形式的统计量检验零
34、假设是否成立，它渐近服从自由度为(k-1)(m-1)的2分布。（2）对称性检验将 Newey 和 Powell(1987)检验最小二乘估计量对称性的方法扩展到分位数回归中。假设我们要检验的分位数回归模型有 m 个，m 是奇数，且中间值(m+1)/2是 0.5，其他都关于 0.5 对称，即j=1 m-j+1,j=1,(m-1)/2。参数估计量按照k的大小排序。则对称性检验的零假设为：)0.5()()(02:1 jmjH(15.32)10其中 j=1,(m1)/2。m 是奇数，代表分位数回归个数。即关于 0.5 对称的分位数回归参数估计量的两两平均值等于中位数回归参数估计量。我们可以构造 Wald
35、 形式的统计量检验上述 k(m-1)/2 个约束条件是否成立。该统计量服从自由度为 k(m1)/2 的2分布。另外，Newey 和 Powell 指出，如果我们已知随机误差项服从独立同分布，但不一定对称，则我们只需检验常数项的对称性。即)0.5(0)(0)(002:1 jmjH(15.33)这时约束条件减少为(m-1)/2 个。15.6 分位数的计算分位数的计算与与分位数回归的分位数回归的 EViews 操作操作（1）分位数的计算对一个离散的随机变量 yt，取其容量为 T 的样本序列，计算第分位数的方法如下：首先将数据从小到大排序，标号为 i，i=1,2,T。然后利用下表所列的方法计算随机变量
36、yt的第分位数的排列序号的 i；如果 i 为整数，则随机变量 yt的第分位数即为 yi，如果 i 不是整数，则随机变量 yt的第分位数为：y()=yi+(i i)(yi+1 yi)其中i表示不大于 i 的最大整数。给定一个具体的随机变量 yt，对于一个容量为 T 的样本，则 yt的第分位数的序号 i 的计算方法如下。在大样本情况下，各方法收敛到同一值。RankitOrdinaryVander WaerdenBlomTukeyGumbel(1/2)/T/T/(T+1)(3/8)/(T+1/4)(1/3)/(T+1/3)(1)/(T1)计算分位数的 EViews 6.0 的命令为：scalar q
37、=quantile(y,s)，其中 y 表示求分位数的序列；表示要取的分位数；s 取 16 依次表示上表中 6 种计算方法，计算所得结果存入标量 q 中。例：打开 6garch-03 文件，在空白处键入命令：scalar q=quantile(DASH,.5,1)scalar q=quantile(DASH,.25,1)意即对序列 DASH 求中位数。得结果 DASH(0.5)=-0.78，DASH 序列的中位数是-0.78。DASH(0.25)=-13.33，DASH 序列的第 0.25 分位数是-13.33。用 DASH 画分位数图如下。打开 DASHt序列窗口，点击 View 键选 Gr
38、aph 功能。在打开的Graph Option 窗口，Type 选择页的 Specifi 选择框选 Distribution，在 Details 的 Distribution 选择框中选 Emprical Quantile 如图。点击“确定”键，得分位数图如图。11-120-80-40040801201600.00.20.40.60.81.0QuantileDASH（2）分位数回归主要包括 3 部分内容。（1）介绍怎样进行分位数回归。（2）对输出结果的分析。（3）对分位数回归相关功能键的介绍。在 EViews 中进行分位数回归的路径有两个，分别是（1）点击主选单中的 Quick 键，选 Equ
39、ation Estimation，弹出 Equation Estimation 窗口。或者（2）点击主选单中的 Object 键，选 New Object，Equation，弹出 Equation Estimation 窗口。在该窗口的 Method 下拉选单中，选择如图所示的选项 QREG-Quantile Regression(includingLAD)，EViews 将打开如图所示的分位数回归对话框（Equation Estimation）。图 1Equation Estimation（方程估计）窗口包括两个选项模块，一个是 Specification（设定方程），一个是 Options
40、（选项）。可以在 Equation specification（方程设定）框中输入要估计的表达式。同一般线性回归模型一样，它可以是一行用空格隔开的被解释变量和解释变量（如图 1 所示），也可以是一个明确的参数为线性的表达式。Equation Estimation（方程估计）窗口与 OLS 估计的 Equation Estimation（方程估计）窗口相比，只多了对话框 quantile to estimate 的选项。在该处填入要估计的分位数。系统默认为 0.5，即做中位数回归（LAD）。用户可以选择任意一个 0 和 1 之间的数（当数值接近 0 和 1 时估计会变得困难）。激活 Option
41、s（选项）模块（点击对话框上的 Options（选项）。得到如图 2 的 quantile regressionOptions（分位数回归选择）选择框、Iteration control（迭代控制）选择框和 Bootstrap settings（自举设定）选择框。12quantile regression Options 对话框中的选择主要包括三部分。图2（1）Coefficient covariance（系数估计量方差协方差矩阵）选项框其下拉选单中包括三个选项：Ordinary(IID)，Huber-Sandwich 和Bootstrap，代表了可选的估计回归系数估计量方差协方差矩阵的方法（
42、具体介绍见15.4节）。EViews 默认的是Huber-Sandwich方法。（2）Weight（权数）选项框可以输入作为权重的序列或者一个序列的表达式，从而对估计式加权。（用于 WLS 估计）（3）Sparsity Estimation（稀疏函数估计）选项区其中包括 5 种选择框。稀疏函数的介绍见 15.4 节。Method（方法）选项框。当第一个选项框 Coefficient covariance 中选项为 Ordinary(IID)或 Bootstrap 时，Method（方法）选项框中包括三个选项：Siddiqui(mean fitted),Kernel(residual)和 Sid
43、diqui(residual)。当 Coefficient covariance 选项框中选项为 Huber-Sandwich 时，这里的 Method 选项框中只包括 Siddiqui(mean fitted)和 Kernel(residual)两种选择。Bandwidth Method（带宽）选项框。其下拉选单中包括四个选择，即 Bofinger(1975)，Hall-Sheather(1988)和 Chamberlain(1994)计算带宽方法，或者你自己给出一个特定的带宽。Size（置信尺度）选项框。当选择 Hall-Sheather 和 Chamberlain 方法时，置信度的选择默
44、认为 0.05。Quantile Method（分位数方法）选项框。EViews 提供了六种求解经验分位数的方法。Kernel（核函数）选项框。表示核函数的选用种类。EViews 中可以选择的核密度函数有 Epanechnikov 核函数、均匀核函数(Uniform)、三角核函数(Triangular)、二权核函数(Biweight)、三权核函数(Triweight)、正态核函数(Normal)和余弦核函数(Cosinus)。注意，不管系数方差协方差矩阵（Coefficient covariance）是否会用到，每次进行分位数回归时，系统都会自动给出一个稀疏函数估计值。Iteration co
45、ntrol（迭代控制）选项框包括3个选项。（1）Max（最大）。迭代的最大次数，默认为 500。（2）Starting（初始值）。表示迭代的初始值，默认为 0，也可以选择其他选项，如下拉选单中的 OLS，即用 OLS 估计量作为初始值进行迭代。（3）Display settings（设定显示）。选择是否需要在输出结果中给出这些设置。13Bootstrap settings（自举设定）（1）Method（方法）。代表不同的自举方法。EViews提供了四种方法，分别是Residual,XY-pair,MCMB,MCMB-A。默认方法为XY-pair方法。（2）Replications（循环次数）。
46、EViews 默认为100次。用户可以自己设定次数。（3）No.of obs（自举样本容量）。空白表示与原样本容量一致。Koenker(2005)的研究表明，选择自举样本容量小于数据样本容量时，能够获得更加准确的结果，特别是当数据样本容量较大时。（4）output（输出）。在这里键入一个名称可以得到自举的样本矩阵。（5）Random generator（生成随机数）和seed（种子）。本选项用于控制产生随机数。其中前者用于选择随机数产生方法，seed用于选择随机数种子，Clear（清除）按钮用于清空以往选定的随机数种子。估计结果。按照 EViews 默认设置得到的一个分位数回归估计结果如下：输
47、出结果上部给出的是估计设定，其中包括(按顺序)被解释变量（DASH）、方法：分位数回归（中位数）、操作日期、样本范围、样本容量（421）、标准误差和方差协方差矩阵估计方法(Huber-Sandwich 方法)、稀疏函数的估计方法（Kernel 方法）、带宽方法（Hall-Sheather 方法，带宽=0.12963）以及对估计结果的评价。输出结果中部给出的是回归系数估计量、标准差、t 统计量及其相应 p 值，这与 OLS 估计完全一样。可以看出，上述回归系数估计量都具有显著性。在中位数回归关系条件下，B 股收益DBSH 每增加一个单位，A 股收益 DASH 平均增加 3.38 个单位。输出结果
48、下部给出的是对分位数回归估计式的评价统计量。分别为Pseudo R-squared：伪拟合优度(伪 R2)，Adjusted R-squared：调整的伪拟合优度，S.E.of regression：分位数回归式的标准误差，Quantile dependent var：分位数回归式中只有常数项存在的系数估计值（也即被解释变量的分位数估计值）。Objective：目标函数极小值，Objective(const.only)：分位数回归式中只有常数存在的目标函数极小值，Sparsity：分位数密度函数（稀疏函数）估计值（本例是用核估计法计算的）。Quasi-LR statistic：准似然比估计量的
49、值14Prob(Quasi-LR stat)：准似然比估计量的值所对应的概率值。此外，由于这里使用的是 Huber-Sandwich 方法，因此稀疏函数值(Sparsity)并没有用来计算参数估计量标准差。与上述结果类似，我们也可以通过改变估计设定，运用自举方法获得参数估计量的方差协方差矩阵。例如选择 MCMB-A 方法进行自举，并且将自举次数增加至 500。对于稀疏函数的计算方法，选择 Siddiqui(mean fitted)，点击 OK 键，得到新设定所对应的估计结果。分位数回归中的 Views 和 Procs 功能键。分位数回归方程窗口中的大部分 Views 和 Procs 功能都与
50、OLS 回归相同，下面对一些计算细节其进行必要的补充说明。使用上述功能时需要注意以下计算细节：（1）这里的残差是指某一特定分位数回归函数条件下的残差，计算公式为ttyu)(-)(X；标准化残差指用自由度调整过的残差的标准误差。而在计算 QLR 统计量时则使用的是 Koenker和 Machado(1999)给出的目标函数极小值的平均值，即)(1)(QT。（2）构造 Wald 检验和置信椭圆时使用的是参数估计量方差协方差矩阵的稳健估计量。（3）进行遗漏和多余变量检验(omitted and redundant tests)以及 Ramsey RESET 检验时，报告的都是特定约束下的 QLR 统