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

輸出：