使用coefplot()的预测器子集

Question

我想用coefplot()绘制一个系数图，它只考虑到我正在使用的预测器的一个子集。例如，如果您有以下代码

 y1 <- rnorm(1000,50,23)
 x1 <- rnorm(1000,50,2) 
 x2 <- rbinom(1000,1,prob=0.63) 
 x3 <- rpois(1000, 2) 
 fit1 <- lm(y1 ~ x1 + x2 + x3)

然后跑了

 coefplot(fit1)

它会给出一个图表，显示拦截、x1、x2和x3的系数。我如何修改它，使我只得到系数，比如说，x1和x2？

Answer 1

您可以使用参数predictors，它只绘制所需的系数：

library(coefplot)
coefplot(fit1, predictors=c('x1','x2'))

输出：

​

​