如何在R的比值比图中展示分类变量参考值及相关疑问

核心问题解答 #

1. 是否应排除截距，参考值怎么展示？ #

必须排除截距项，因为截距对应的不是“参考组相对于自身的OR”，而是参考组自身的患病优势（odds），这并非OR图要呈现的核心内容。OR图的本质是展示各分类水平相对于参考水平的比值比，参考水平的OR固定为1，但它没有对应的回归系数（建模时被设为基准），正确做法是：

• 从绘图数据中移除截距行
• 在图中为参考组保留位置，可标注“参考组”或仅留空，让读者明确其作为对比基准的角色

2. 指数化截距不等于1是否有问题？ #

完全没问题。原因如下：

• 截距的指数值是参考组的患病优势（odds = P(患病)/P(未患病)），而非比值比（OR）
• 参考组与自身对比的OR确实为1，但这个值不会出现在回归结果中（参考组的系数被设为0，指数化后为1，但回归输出无此行），你看到的截距指数值是参考组本身的优势，不是OR，所以不等于1是正常现象

代码修改示例 #

针对你的代码，只需在生成绘图数据B后，过滤截距项并调整绘图逻辑即可：

# 过滤截距项，仅保留分类变量的非参考水平
B_filtered <- B %>% filter(term != "(Intercept)")

# 修正varname的对应关系（更严谨的方式）
B_filtered$varname <- ifelse(grepl("agecat", B_filtered$term), "Age", "Gender")

# 绘制优化后的OR图
ggplot(B_filtered, aes(y = term, x = estimate)) +
  geom_point(shape = 18, size = 3) +
  geom_errorbarh(aes(xmin = conf.low, xmax = conf.high), height = 0.4, size=0.2) +
  geom_vline(xintercept = 1, color = "red", linetype = "dashed", alpha = 0.5) +
  xlab("Odds Ratio (95% CI)") +
  facet_col(~varname, scales = 'free_y', space='free') +
  # 可选：为参考组添加标注（以Age组的0-19为例）
  annotate("text", x = 1, y = "agecat0-19", label = "参考组", color = "gray50", size = 3)

额外说明 #

• 确保分类变量的水平顺序在因子设置时正确，这样y轴顺序会与你设定的参考组顺序一致
• 若想更清晰展示参考组，可在y轴标签中明确标注，或用灰色文本等样式区分

内容的提问来源于stack exchange，提问作者Lindsay