（R）用quantreg包绘制分位数回归系数分布

Question

  

可能重复：
  Multiple plots with high-level plotting functions, especially plot.rqs()

我正在使用一个依赖变量和一个自变量运行18个分位数回归。我希望在$ \ tau = 0.01,0.02，...，0.99 $的18个斜率估计值的分布中得到$ 6 * 3 $ tile图。

这是我到目前为止所能做到的：

require(quantreg)

op <- par(xlab="")
y=rnorm(200)
x=rnorm(200)
plot(summary(rq(y~x,tau=1:99/100)),parm=2,main="",mfrow=c(6,3),mar=c(0.5,2,2,2))
par(new=TRUE)
plot(summary(rq(y~x,tau=1:99/100)),parm=2,main="",mfrow=c(6,3),mar=c(0.5,2,2,2))

第二个情节一直到右下角，当我第三次尝试时，它甚至不在图形窗口中。

我试图将par（mfrow = c（6,3））放在plot（）之外，但它不起作用。

Answer 1

plot.summary.rqs的默认行为是构建一个瓦片图，每个参数都有一个图，这可能是你习惯看到的，但正如Roland所说，它正在以一种方式使用par（）打败你的努力。你可以通过几种方式克服这个问题。一个是重做该功能。但是如果你想保持这个功能没有改变，你可以制作一个辅助功能（带有新名字的原件的副本），用于注释par调用（有2个）并将其用于随后的框架内绘图。

Answer 2

首先，我会担心这些警告。

问题是plot.summary.rqs设置并重置了绘图参数，这不是很好的imho。一种解决方案是在没有此功能的情况下基于plot.summary.rqs定义您自己的绘图功能。