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

输出：