r - (R) 使用 quantreg 包绘制分位数回归系数分布

Question

可能重复：
具有高级绘图功能的多个绘图，尤其是 plot.rqs()

我正在运行 18 个分位数回归，其中包含一个因变量和一个自变量。我想要 18 个斜率估计值在 $\tau = 0.01,0.02,...,0.99$ 上的分布的 $6 * 3$ 平铺图。

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

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设置和重置绘图参数，这不是很好的恕我直言。一种解决方案是在plot.summary.rqs没有此功能的情况下定义您自己的绘图功能。