如何用ggplot2修改森林图中特定疾病对应的点形状？

核心思路 #

给数据框添加一个标记列，区分需要特殊形状的行，然后在ggplot中通过shape映射实现不同点样式。

步骤1：处理数据 #

先给你的数据框新增一列，标记每个行对应的点形状：

# 假设你已经读取了Excel数据到df中
# 先安装并加载readxl包（如果还没装）
# install.packages("readxl")
library(readxl)
df <- read_excel("你的Excel文件路径.xlsx")

# 添加形状标记列：默认用实心圆点（形状代码16）
df$point_shape <- 16
# 指定第6、8、14、16行改为实心方形（形状代码15）
df$point_shape[c(6, 8, 14, 16)] <- 15

形状代码参考：15=实心方形，17=实心三角形，21=空心方形，22=空心菱形，可根据需求替换。

步骤2：修改绘图代码 #

在ggplot的aes中加入shape映射，再手动设置样式并隐藏不必要的图例：

library(ggplot2)

fp <- ggplot(data=df, aes(x=Disease, y=OR, ymin=Lower, ymax=Upper, shape=factor(point_shape))) +
  geom_pointrange() + 
  geom_hline(yintercept=1, lty=2) +  # 添加y=1的参考虚线
  coord_flip() +  # 翻转坐标轴，将疾病名称放在y轴
  xlab("Disease") + ylab("OR (95% CI)") +
  theme_bw() +
  # 手动指定形状，隐藏图例（不需要展示形状分组）
  scale_shape_manual(values=c(16, 15), guide="none")

print(fp)

注意事项 #

如果你的目标行是按疾病名称而非原始行号确定的，可改用以下方式标记：

# 先获取目标疾病名称
target_diseases <- df$Disease[c(6, 8, 14, 16)]
# 根据疾病名称匹配修改形状
df$point_shape[df$Disease %in% target_diseases] <- 15

这种方式更稳妥，避免数据行顺序变动导致错误。

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