如何在R语言的glm模型输出中显示参考变量？

方案1：用parameters包快速实现 #

这是最省心的方法，不需要手动修改输出，直接就能把参考水平的行加上：

先安装并加载包：

install.packages("parameters")
library(parameters)

然后对你的模型执行参数提取：

model <- glm(y ~ x1_1 + x1_2, data = data, family = "binomial")
parameters(model, include_reference = TRUE)

这个函数会自动把参考组x1_0的行插入到输出里，明确标记出参考水平，完全不会改动原模型的截距和变量解释逻辑。

方案2：手动修改summary输出（无额外包） #

如果不想装新包，直接修改summary的系数表就行：

model <- glm(y ~ x1_1 + x1_2, data = data, family = "binomial")
sum_out <- summary(model)
# 生成参考水平的空行
ref_row <- data.frame(Estimate = NA, `Std. Error` = NA, `z value` = NA, `Pr(>|z|)` = NA, row.names = "x1_0")
# 把参考行插入到截距和其他自变量之间
sum_out$coefficients <- rbind(
  sum_out$coefficients[1,,drop=FALSE],
  ref_row,
  sum_out$coefficients[-1,,drop=FALSE]
)
# 打印修改后的结果
sum_out

这样输出的系数表就和你期望的样式完全一致，保留了原模型的所有设定，只是多了参考水平的行来明确标注。

内容的提问来源于stack exchange，提问作者j-indigo