分位数回归—R语言实现

（备考中，暂停更新）4.14 于武汉

已于 2023-11-20 16:40:15 修改

阅读量4.9k

点赞数 16

文章标签： r语言回归

于 2023-03-03 15:24:55 首次发布

本文链接：https://blog.csdn.net/qq_45856698/article/details/129319946

版权

大家好，我是带我去滑雪，每天教你一个小技巧！

分位数回归—R语言实现

1、分位数基本介绍——什么是分位数？

2、分位数回归用途——为什么要使用分位数回归？

3、图形分析——如何进行分位数回归图像分析？

4、分位数回归原理——如何进行分位数回归？

5、分位数回归的假设检验——能不能用？好不好？

6、简单线性分位数回归实例——分位数回归R实操

1、分位数基本介绍——什么是分位数？

2、分位数回归用途——为什么要使用分位数回归？

传统的线性回归模型描述了因变量的条件均值分布受自变量X的影响过程。其中，最小二乘法是估计回归系数的最基本方法。如果模型的随机误差项来自均值为零、方差相同的分布，那么回归系数的最小二乘估计为最佳线性无偏估计；如果随机误差项进一步服从正态分布，那么回归系数的最小二乘估计与极大似然估计一致，均为最小方差无偏估计。此时它具有无偏性、有效性等优良性质。

但是在实际的经济生活中，这种假设通常不能够满足。例如当数据中存在严重的异方差，或后尾、尖峰情况时，最小二乘法的估计将不再具有上述优良性质。最小二次估计假定解释变量x只能影响被解释变量的条件分布的均值的位置，不能影响其分布的刻度或形状的任何其他方面。分位数回归的条件更加宽松。

3、图形分析——如何进行分位数回归图像分析？

假如，现在有一个如下散点图分布的数据，对其进行普通的回归分析，得到：

从拟合的曲线我们就可以看出问题了，原数据随着x增大，y的分布范围越来越大，可是因为普通的回归分析得到的是条件期望函数，也就是y的期望，所以平均即使y的分布变化了，平均来说y还是以同样的斜率稳定上升。

当我们使用0.9分位数回归，重新得出新函数图像：

这次，比起普通的回归分析，就能进一步显示出y的变化幅度其实是增大了。所谓的0.9分位数回归，就是希望回归曲线之下能够包含90%的数据点（y），这也是分位数的概念，分位数回归是把分位数的概念融入到普通的线性回归而已。

进一步的我们可以画出不同的分位数回归曲线，这样才能能更加明显地反映出，随着x的增大，y的不同范围的数据是不同程度地变化的，而这个结论通过以前的回归分析是无法得到的，这就是分位数回归的作用。

在实际研究中，例如研究社会的贫穷情况可能就会选择收入相对较低的群体即收入水平的低分位点,而研究制定税收政策会关注高收入人群即收入水平的高分位点。

4、分位数回归原理——如何进行分位数回归？

5、分位数回归的假设检验——能不能用？好不好？

分位数回归估计的检验包括两部分：一是与均值回归类似的检验，例如拟合优度检验、回归约束检验；二是分位数回归估计特色要求的检验，如斜率相等检验、斜率对称性检验。

6、简单线性分位数回归实例——分位数回归R实操

install.packages("quantreg")

library(quantreg)#导入分位数回归的包

data(engel)# 引入数据

mode(engel)#查看数据格式

names(engel)#查看变量名

head(engel)#查看数据的前五行

plot(engel$income,engel$foodexp)#画个散点图

可以看出，数据随着x的增大，y的分布范围越来越大。

#简单验证一下因变量foodexp是否服从正态分布

plot(density(engel$foodexp))

可以看出y是右偏分布。

qqnorm(engel$foodexp, main='QQ plot')

qqline(engel$foodexp, col='red', lwd=2)

结果呈一条直线才能说明大致服从正态分布，显然因变量foodexp明显不服从正态分布，但是，分位数回归不要求y服从正态分布，不仅如此，而且分位数回归还对异常值点不敏感。

下面我们继续，为了对比，我们做一个均值回归，再做一个分位数回归。

#设置0.05, 0.25, 0.5, 0.75, 0.95五个分位点，并且rq函数进行分位数回归，这样可以得到五条分位数回归线

rq_result <- rq(foodexp ~ income, tau=c(0.05, 0.25, 0.5, 0.75, 0.95))

summary(rq_result,se="nid")#系数的假设检验，se指定用于计算标准误差的方法

#上面就是每个分位点下回归线的回归系数，做个图看一下

plot(income, foodexp, cex=0.25, type='n', xlab='income', ylab='foodexp')

points(income, foodexp, cex=0.5, col='blue')

abline(rq(foodexp~income, tau=0.5), col='blue')#加中位数数回归的直线

abline(lm(foodexp~income), lty=2, col='red')#加均值回归线

taus <- c(0.05, 0.1, 0.25, 0.75, 0.9, 0.95)#将分位数回归线加上去

for (i in 1:length(taus)){

abline(rq(foodexp~income, tau=taus[i]), col='gray')