【发布时间】:2013-08-13 13:08:39
【问题描述】:
我正在为使用quantreg 进行的分位数回归绘制回归摘要。
显然这里使用了plot.summary.rqs 方法。问题是使用了很多解释变量,每个解释变量都显示在图中。大多数系数的表现与 OLS 没有显着差异,所以我只想挑选并显示其中的一些。
如何选择需要显示的图?我在报告中使用knitr,但不想显示数十个变量(并且您可以使用假人快速到达那里)。有办法摘樱桃吗?
【问题讨论】:
我正在为使用quantreg 进行的分位数回归绘制回归摘要。
显然这里使用了plot.summary.rqs 方法。问题是使用了很多解释变量,每个解释变量都显示在图中。大多数系数的表现与 OLS 没有显着差异,所以我只想挑选并显示其中的一些。
如何选择需要显示的图?我在报告中使用knitr,但不想显示数十个变量(并且您可以使用假人快速到达那里)。有办法摘樱桃吗?
【问题讨论】:
默认情况下,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
【讨论】: