用R语言进行分位数回归：基础篇詹鹏（师大学经济管理学院）本文根据文献资料整理，以介绍方法为主要目的作者的主要贡献有：（1）整理了分位数回归的一些基本原理和方法；（2）归纳了用R语言处理分位 数回归的程序，其中写了两个函数整合估计结果；（3）写了一个分位数分解函数来处理MM2005的分解过程；（4）使用一个数据集进行案例分析，完整地展现 了分析过程第一节分位数回归介绍（一）为什么需要分位数回归？传统的线性回归模型描述了因变量的条件均值分布受自变量X的影响过程其中，最小二乘法是估计回归系数的最基本方法如果模型的随机误差项来 自均值为零、方差相同的分布，那么回归系数的最小二乘估计为最佳线性无偏估计（BLUE）；如果随机误差项是正态分布，那么回归系数的最小二乘估计与极 大似然估计一致，均为最小方差无偏估计（MVUL）此时它具有无偏性、有效性等优良性质但是在实际的经济生活中，这种假设通常不能够满足例如当数据中存在严重的异方差，或后尾、尖峰情况时，最小二乘法的估计将不再具有上述优良 性质为了弥补普通最小二乘法（OLS）在回归分析中的缺陷，1818年Laplace[2]提出了中位数回归（最小绝对偏差估计）。
在此基础上，1978年Koenker …[3] 和Bassett 把中位数回归推广到了一般的分位数回归（Quantile Regression）±0分位数回归相对于最小二乘回归，应用条件更加宽松，挖掘的信息更加丰富它依据因变量的条件分位数对自变量X进行回归，这样得到了所有分位数 下的回归模型因此分位数回归相比普通的最小二乘回归，能够更加精确第描述自变量X对因变量Y的变化围，以及条件分布形状的影响二）一个简单的分位数回归模型[4]假设随机变量的分布函数为（1）的最小 值，即Y的 分位数的定义为满足（3）样本中位数回归是使误差绝对值之和最小，即（4）样本分位数回归是使加权误差绝对值之和最小，即（5）上式可等价表示为：其中，为检查函数（check function），定义为：其中，；当z为假时，同线性方程y=kx比较，为指示函数(indicator function)，z是条件关系式，当z为真时，相当于直线的斜率k，可以看出，为分段函数，如下图所示通过求解(8)式得到参教估计值现假设因变量Y由k个自变量组成的矩阵X线性表示，对于条件均值函数对于条件分位数函数，通过求解(9)式得到参教估计值函数表示取函数最小值时的取值。
式中，(三) 分位数回归模型的参数估计算法1、主要算法(1)单纯形算法(Simplex Method)一 …[5" 「61. Koenker和Orey (1993)把分两步解决最优化问题的单纯形算法 扩展到所有回归分位数中该算法估计出来的参数具有很好的稳定性，但 是在处理大型数据时运算的速度会显著降低2)点算法(Interior Point Method)由于单纯形算法在处理大型数据时效率低下，Karmarker提出了点算法 [7] :Portnoy和Koenker把这种方法是用在分位数回归中，得出了处理大型数据时点算法的运算速度远快于单纯形算法的结论但点算法每计算一步都要进行因数分解，当自变量比较多的时候效率比较低其次，如果要达到和单纯形算 法一样的精度，就必须进行舍入步骤的计算，者也降低了算法的运行效率3)平滑算法(Smoothing Method)上述两种算法都有各自的优点和不足，而有限平滑算法则是一种同时兼顾运算效率以及运算速度的方法[Chen把这种算法扩展到计算回归分位数中「8]2、R语言quantreg包中的假设检验加载quantreg包以后，使用summary。
函数或summary.rq（）函数，可以得到参数系数的一些假设检验统计量其实，以上两个函数是一致的在使用summary的时候，如果sumamry（）加载的模型（对象）是分位数回归模型，则会自动调用summary.rq（）来处理这个对象summary.rq（）的调用格式为summary(object, se 二 NULL, covariance=FALSE, hs 二 TRUE,...)其中主要参数有：# object:分位数回归对象，根据rq（）函数等得到的结果 se:用于计算参数估计值标准差的方法，可以选取的值包括：rank:根据Koenker（1994）的秩检验得到标准差的估计值默认情况下假定残差是服从独立同分布如果补充另一个参数iid=FALSE，则采用Machado（1999）的方法计算标准差（参数的写法：se=” rank”，iid=FALSE）iid:（这个与上面提到的iid二FALSE不同，这里是参数se的一个取值，而上面的iid是一个逻辑参数）假定残差服从独立同分布，并按照KB（1978）的方法计算残差nid:用sparsity算法计算的参数估计值标准差ker:用Powell（1990）的核密度估计方法得到标准差。
boot:采用bootstrap自助抽样的方法计算标准差默认情况下，se=NUL L且convariance二FALSE，标准差的默认算法是se=” rank” ；其他情况下，se默认值为” nid” covariance:逻辑参数，是否返回参数估计量的协方差矩阵不同参数的结果，可参看下面的程序案例 、 「91（四）分位数分解（MM2005方法）我们可以进一步运用分位数分解法对各个影响因素进行分解分析[10]这里仅介绍MM2005方法为讲解方便，这里以各因素对城乡家庭收入的影响为例，观察各个影响因素在不同分位数上对城乡家庭收入差异的影响度的大小这里介绍Machado 和 Mata[11] （2005）提出的分位数分解法，将每个分位数上的城乡收入差异分解为两个部分：一部分是由于城乡家庭劳动力特征的不同回报率引起的（即分位数回归参数的不同引起的，The Return Effects），例如城乡家庭劳动力在相同的教育程度、工作年限以及所处当地的经济发展水平相同的特定因素下不同的回报率引起的家庭人均收入差异；另一部分是由于城乡家庭劳动力的特征变量分布不同引起的（即影响因素变量值的不同引起的，The Covariate Effect），城乡家庭人均收入这部分的差异会随着样本分布的不同而略有变化。
利用Machado和Mata分位数分解方法的关键是进行反事实分析（the counter-factual analysis），我们最关心的一种反事实分析就是，如果城市家心肥K庭劳动力按照农村家庭劳动力的分位数回归参数决定家庭人均收入的话，城市家庭的人均收入分布会如何？这里定义反事实分布为其中表示影响农村家庭人均收入的变量在每个分位数上的回归参数 表示如果城市家庭劳动力按照表示影响城市家庭人均收入的变量分布，农村家庭劳动力的分位数回归参数决定家庭人均收入的话，城市家庭的反事实人均收入的大小的具体计算步骤为：（1）确定不同的分位点，分别表示4 ）把（2）中得到的分位数回归参数和（3）中得到得城市家庭子样本变量分布相结合，得到一个新的样本，即反事实分布样本2）在农村家庭样本中，分别以做分位数回归，得到组分位数回归参数向量（3）将城市家庭样本数据表示为假定在t分位数下城市家庭人均收入、反事实家庭人均收入和农村家庭人均收入分别为则不同分位数下的城乡家庭人均收入分布差异可表示为：等式右边的第一项称为“回报影响（the return effect）”，它表示在不同的分位数下，由于城乡家庭劳动力的生产回报率不同所导致的城乡差异部分；等式右边第二项成为“变量影响（the covariate effect）”，它表示不同分位数下城乡家庭随机抽样的样本变量分布不同所导致的城乡差异部分。
五）非线性分位数回归和非参数分位数回归暂略第二节用R语言进行分位数回归(一)安装和加载包R语言的基本包中没有进行分位数回归的程序包，故需要在官网下载并安装相应的程序包quantreg在电脑上安装过quantreg包以后，下次不需要再次安 装了但每次使用分位数回归前，需要加载quantreg包install.package( "quantreg” ) #保持联网的情况下安装包library( "quantreg” ) # 加载包help.start()help(rq)example(rq)#进入R帮助首页#获取rq函数的帮助，也可以写成：?rq#显示分位数回归函数rq()的一个简单示例代码(二)一个简单的分位数回归模型及结果data(engel)#加载quantreg包自带的数据集，见说明①fit1 = rq(foodexp 〜income, tau = 0.5, data = engel) #进行分位数回归，见说明②fit1#直接显示分位数回归的模型和系数，见说明③summary(fit1)#得到更加详细的显示结果，见说明④r1 = resid(fit1)#得到残差序列，并赋值为变量r1c1 = coef(fit1)#得到模型的系数，并赋值给变量c1，见说明⑤summary(fit1, se =“nid” ) #通过设置参数se，可以得到系数的假设检验，说明⑥说明：①engel (1857)是考察食物支出与家庭收入之间关系的一个数据集，用函数head(engel)可以查看前六行的值：# income foodexp# 1 420.1577 255.8394# 2 541.4117 310.9587# 3 901.1575 485.6800# 4 639.0802 402.9974# 5 750.8756 495.5608# 6 945.7989 633.7978② 这里因变量为foodexp，即食物支出。
自变量为income，即家庭收入tau表示计算50%分位点的参数，这里可以同时计算多个分位点的分位数回归结果，如tau=c(0.1,0.5,0.9)是同时计算10%、50%、90%分位数下的回 归结果data二engel指明这里处理的数据集为engelmethod :进行拟合的方法，取值包括：A.默认值“br”，表示Barrodale & Roberts算法的修改版；B. “fn”，针对大数据可以采用的Frisch - Newton点算法；C. “pfn” ，针对特别大数据，使用经过预处理的Frisch - Newton逼 近方法；D. “fnc”，针对被拟合系数特殊的线性不等式约束情况；E. “lasso”和“ scad”， 基于特定惩罚函数的平滑算法进行拟合③ 直接运行fitl，会得到简单的计算结果，如：# Call:# rq(formula 二 foodexp 〜income, tau = 0.5, data 二 engel)## Coefficients:# (Intercept) income# 81.4822474 0.5601806## Degrees of freedom: 235 total; 233 residual④ 用summary。
函数可以得到回归模型的详细结果，包括系数和上下限 Call: rq(formula 二 foodexp 〜income, tau = 0.5, data 二 engel)## t。