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

时间:2013-08-12 09:12:29

标签: r plot knitr quantreg

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

如何选择我需要显示的图?我使用knitr作为我的报告,但不想显示几十个变量(你很快就会使用假人)。有没有办法挑选樱桃?

1 个答案:

答案 0 :(得分:0)

默认情况下,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