Question

我只得到coefplot部分自变量。我的回归方程是固定效应回归如下：

aa1 <-glm(Eighty_Twenty ~ Market_Share_H+Market_Share_L+Purchase_Frequency_H+Purchase_Frequency_L+factor(product_group))
coefplot(aa1)

但是，我不想绘制factor(product_group)变量的系数，因为有产品组。相反，我会得到一个只有其他变量系数的系数图。我怎么能这样做？

Answer 1

在帮助页面（请参阅?coefplot.default）中，您可以选择绘图中所需的预测变量或系数。

# some example data
df <- data.frame(Eighty_Twenty = rbinom(100,1,0.5),
                 Market_Share_H = runif(100),
                 Market_Share_L = runif(100),
                 Purchase_Frequency_H = rpois(100, 40),
                 Purchase_Frequency_L = rpois(100, 40),
                 product_group = sample(letters[1:3], 100, TRUE))

# model
aa1 <- glm(Eighty_Twenty ~ Market_Share_H+Market_Share_L +
                           Purchase_Frequency_H + Purchase_Frequency_L + 
                           factor(product_group), df, family="binomial")


library(coefplot)

# coefficient plot with the intercept
coefplot(aa1, coefficients=c("(Intercept)","Market_Share_H","Market_Share_L",  
                              "Purchase_Frequency_H","Purchase_Frequency_L"))

# coefficient plot specifying predictors (no intercept)
coefplot(aa1, predictors=c("Market_Share_H","Market_Share_L"  ,      
                             "Purchase_Frequency_H","Purchase_Frequency_L"))

R中的coefplot与部分独立变量

1 个答案: