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

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