技术求助：如何使用R ggplot绘制指定样式的柱状图？

Hey Peter, 别着急！这种带置信区间的分组OR值柱状图在R里完全可以实现，最灵活的工具就是ggplot2包。先帮你梳理下步骤和代码：

首先，先把你的数据整理成整洁数据格式（tidy data），因为你现在的A、B组数据是分开的，合并后更方便绘图。先补全你没写完的OR_B相关数据（你可以根据实际值替换）：

# 加载所需包
library(tidyverse)

# 原始数据
groupe2 <- rep(c(rep("P",4), rep("I",4)), 2)
groupe <- rep(c("PPP","PPI","PIP","PII","IPP","IPI","IIP","III"), 2)
OR_A <- c(1.00,0.86,0.88,0.90,0.77,0.68,0.77,0.70)
ICinf_A <- c(NA,0.70,0.72,0.76,0.60,0.50,0.61,0.61)
ICsup_A <- c(NA,1.06,1.07,1.06,1.00,0.92,0.96,0.81)
# 补全OR_B及置信区间（请替换为你的实际数据）
OR_B <- c(1.00,0.97,1.01,0.81,0.75,0.69,0.72,0.67)
ICinf_B <- c(NA,0.82,0.85,0.68,0.60,0.55,0.59,0.58)
ICsup_B <- c(NA,1.14,1.20,0.97,0.93,0.86,0.88,0.78)

# 合并成整洁数据框
df <- tibble(
  groupe = groupe,
  groupe2 = groupe2,
  OR = c(OR_A, OR_B),
  ICinf = c(ICinf_A, ICinf_B),
  ICsup = c(ICsup_A, ICsup_B),
  # 添加A/B组标识
  dataset = rep(c("A", "B"), each = 8)
)

接下来就可以用ggplot2绘图了，这里提供两种贴合你数据结构的常见样式思路：

思路1：按A/B组分面展示 #

如果想把A、B组的图分开对比，用分面（facet）会很清晰：

ggplot(df, aes(x = groupe, y = OR, fill = groupe2)) +
  # 绘制柱状图，调整宽度避免重叠
  geom_col(position = position_dodge(width = 0.8), width = 0.7) +
  # 添加置信区间误差棒
  geom_errorbar(aes(ymin = ICinf, ymax = ICsup), 
                position = position_dodge(width = 0.8), width = 0.2) +
  # 分面展示A/B组
  facet_wrap(~dataset) +
  # 添加OR=1的参考线（常用于对比基准）
  geom_hline(yintercept = 1, linetype = "dashed", color = "red") +
  # 调整标签和样式
  labs(title = "OR值及置信区间", x = "分组", y = "OR值") +
  theme_bw() +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))

思路2：同一图内对比A/B组，按groupe2分面 #

如果想在同一视图里直接对比A、B组，同时保留groupe2的分组逻辑，可以这样调整：

ggplot(df, aes(x = groupe, y = OR, fill = dataset)) +
  geom_col(position = position_dodge(width = 0.8), width = 0.7) +
  geom_errorbar(aes(ymin = ICinf, ymax = ICsup), 
                position = position_dodge(width = 0.8), width = 0.2) +
  # 按groupe2分面展示
  facet_wrap(~groupe2) +
  geom_hline(yintercept = 1, linetype = "dashed", color = "red") +
  labs(title = "OR值及置信区间（按groupe2分组）", x = "分组", y = "OR值", fill = "数据集") +
  theme_bw() +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))

关键细节说明： #

• 你数据里ICinf_A第一个值是NA（对应OR=1.00的参考组），geom_errorbar会自动忽略这个NA值，不会显示误差棒，这完全符合医学/统计领域展示参考组的习惯。
• 如果你的目标样式有特殊要求（比如柱状图排列方式、颜色搭配、添加数值标注等），可以补充细节，我再帮你调整代码。

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