急！！R语言面板分位数回归的系数图如何做？

急！！R语言面板分位数回归的系数图如何做？
类似于这种图

使用RStudio会方便很多：
1、安装R3.5.1
2、安装RStudio（最新版本1.2.1335 ：https://www.rstudio.com/products/rstudio/download/）
3、设置RStudio中的R路径以及CRAN镜像站点（推荐清华镜像站）
4、在RStudio里面下载Quantreg包
5、拟合/输出结果

详见下。Hope it helps

怎么看不到刚刚输入的信息呢？

#安装quantreg分位回归包
install.packages("quantreg")
#载入quantreg包
library("quantreg", lib.loc="~/R/win-library/3.5")
#看一下分位回归包里面的恩格尔演示数据
View(engel)
#画一下恩格尔演示数据
plot(engel)
#对数据的收入进行按平均值归一
engel<-within(engel,xx <- income - mean(income))
#分位拟合0.5分位值处的食品支出与收入的关系，显示分位拟合系数结果
rq(foodexp~income,tau=0.5,data = engel)
#以百分之一的间隔分位拟合0.02~0.98的食品花销~收入，显示分位拟合系数结果
rq(foodexp~income,tau=1:19/20,data = engel)
#以百分之五的间隔分位拟合0.05~0.95的食品花销~收入
rqFit=summary(rq(foodexp~xx,tau=1:19/20,data = engel))
#分位拟合系数画图
plot(rqFit,mfrow = c(1:2))

请用你自己的数据代替上述的engel数据以及分位间隔值

liusm99 发表于 2019-9-2 20:09
#安装quantreg分位回归包
install.packages("quantreg")
#载入quantreg包

你好，我当时用的是rqpd，代码如下：
a.form<-y~x
h.a<-rqpd(a.form,panel(lambda=0.5,taus= c(0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9), tauw=rep(1/10, 9)),data=data,control=list(tmpmax=10000))
a<-summary(h.a)
plot(a,mfrow = c(1:2))
显示为Error in xy.coords(x, y, xlabel, ylabel, log) :
  'x' is a list, but does not have components 'x' and 'y'

*画变量x1的0.05~0.95分位回归系数图(另一种作法)
forv i=5(5)95{
qreg y x*, q(`i')
predictnl b`i'=_b[x1], ci(l`i' u`i')
}
keep b* l* u*
keep in f
xpose,clear v
destring _v, g(q) i("b l u")
replace  _v=substr(_v,1,1)
reshape wide v1, i(q) j(_v) s
tw rarea v1l v1u q || line v1b q

邮箱多少。谢谢！

fjpsf 发表于 2019-9-3 20:30
邮箱多少。谢谢！

834596658@qq.com
我想问下R语言如何操作呢？

fjpsf 发表于 2019-9-3 20:30
邮箱多少。谢谢！

743491270@qq.com  求示例或相关代码

请问stata软件怎么画这个图呢？