r - 在分位数回归中显示单系数图？

Question

我正在为我使用的分位数回归绘制回归摘要quantreg。显然plot.summary.rqs这里使用了这个方法。问题是使用了很多解释变量，每个解释变量都显示在图中。大多数系数的表现与 OLS 没有显着差异，所以我只想挑选并显示其中的一些。

如何选择需要显示的图？我正在使用knitr我的报告，但不想显示几十个变量（你可以使用假人快速到达那里）。有办法摘樱桃吗？

Answer 1

默认情况下，plot.summary.rqs绘制所有系数：

library(quantreg)
data(stackloss)

myrq <- rq(stack.loss ~ Air.Flow + Water.Temp + Acid.Conc., tau = seq(0.35, 0.65, 0.1), data=stackloss)
plot(summary(myrq)) # Plots all 4 coefficients

要挑选系数，parm可以使用该参数：

plot(summary(myrq), parm = 2) # Plot only second regressor (Air.Flow)
plot(summary(myrq), parm = "Water.Temp") # Plot only Water.Temp
plot(summary(myrq), parm = 3:4) # Plot third and fourth regressor