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

*******************
sysuse auto
preserve
loc s=0
mat b=J(19,4,.)
forv i=0.05(0.05)1{
loc s=`s'+1
qreg price weight length foreign , q(`i')
predictnl b=_b[weight], 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 lo u
}

clear
svmat b
tw rarea b3 b4 b1||scatter b2 b1,c(l)
restore

*********************

我感觉程序应该是这样的：
*******************
sysuse auto
preserve
loc s=0
mat b=J(19,4,.)
forv i=0.05(0.05)1{
loc s=`s'+1
qreg price weight length foreign , q(`i')
predictnl b=_b[weight], 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 lo up
}

clear
svmat b
tw rarea b3 b4 b1||scatter b2 b1,c(l)
restore

*********************

太有用了  谢谢各位~~！谢谢

非常感谢分享

lihoujian 发表于 2013-6-8 10:47
我感觉程序应该是这样的：
*******************
sysuse auto

如果变量中没有与up开头字母重复的，u可以自动识别成up

太有用了，非常感谢各位高手的命令分享

sungmoo 发表于 2012-11-12 13:51
*画变量x1的0.05~0.95分位回归系数图(另一种作法)
forv i=5(5)95{
qreg y x*, q(`i')

请教下，您的命令好像这一步有点问题
destring _v, g(q) i("b l u")
提示 _varname contains characters not specified in ignore(); no generate

thanks for sharing

h3327156 发表于 2011-8-13 00:18
这一题测好久…
我的感觉啦！
sungmoo版主用法都很精简。

drop b lo u的u应该为up,不然会出u ambiguous abbreviation的警告

voodoo 发表于 2010-4-22 12:29
有个现成的命令grqreg帮你得到你要的
ssc install grqreg
help grqreg

厉害，大神！