如何画分位数回归后的这种图

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

比如，一组数据x和y。进行分位数回归（0.1，0.2，0.3，0.4，0.5，0.6，0.7，0.75，0.8，0.9，0.95），将得到一系列的截距和斜率的值，以及它们的置信区间。请问，如何将相关的数据提取出来，并用以绘制如下图形？谢谢先。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏22 回帖

关键词：分位数回归 分位数 数据提取 置信区间 位数

这个真的有点难

有个现成的命令grqreg帮你得到你要的
ssc install grqreg
help grqreg

*手工操作（0.05~0.95，19个分位点，x的估计系数的图像）
preserve
loc s=0
mat b=J(19,4,.)
forv i=0.05(0.05)1{
loc s=`s'+1
qreg y x, q(`i')
predictnl b=_b[x], ci(lo up)
mat b[`s',1]=`i'
mat b[`s',2]=b[1]
mat b[`s',3]=lo[1]
mat b[`s',4]=up[1]
drop b l u
}
clear
svmat b
tw rarea b3 b4 b1||scatter b2 b1,c(l)
restore

sungmoo 发表于 2010-4-22 12:53
*手工操作（0.05~0.95，19个分位点，x的估计系数的图像）
preserve
loc s=0
mat b=J(19,4,.)
forv i=0.05(0.05)1{
loc s=`s'+1
qreg y x, q(`i')
predictnl b=_b[x], ci(lo up)
mat b[`s',1]=`i'
mat b[`s',2]=b[1]
mat b[`s',3]=lo[1]
mat b[`s',4]=up[1]
drop b l u
}
clear
svmat b
tw rarea b3 b4 b1||scatter b2 b1,c(l)
restore

predictnl b=_b[x], ci(lo up)
这句真是很妙！怎么想到的？！

感谢3楼voodoo和荣誉版主！我赶紧试试去。

非常感谢，感谢指导

看后很有启发，特别是4楼那段手工程序呀

谢谢楼主。。好厉害啊~~初学者很受用~~

是噢，真的很厉害