使用 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'))

输出：

Subset of predictors using coefplot()