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

Question

这个问题在这里已经有了答案:




9年前关闭。

Possible Duplicate:
Multiple plots with high-level plotting functions, especially 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 个)并将其用于框架内的后续绘图。