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

输出：