复现顶刊图，换数据即可用！

pythonic生物人2026年06月03日约 2 分钟

本次使用R语言复现一下Nature上的图！

✅读入测试数据

# 读入测试数据
plot_data <- read.csv("r_data/volcano_y.csv", header = TRUE)

head(plot_data)

✅绘制中间灰色点图

geom_point(data = filter(plot_data, group == "Neutral"),
             colour = "black", alpha = 0.4, size = 2)

✅绘制上侧点图

geom_point(data = filter(plot_data, group == "Up"),
             colour = "black", shape = 21, size = 3, stroke = 0.6, fill = "red")

✅绘制下侧点图

geom_point(data = filter(plot_data, group == "Down"),
             colour = "black", shape = 21, size = 3, stroke = 0.6, fill = "#2d2d9c")

✅绘制中间辅助线

geom_hline(yintercept = 0, linewidth = 0.4, colour = "black")

✅绘制文字标签

geom_text_repel(data = filter(plot_data, label != ""), aes(label = label),
                  family = "Arial", fontface = "bold.italic",
                  nudge_x = ifelse(fold_change > 0, -2, 5),
                  nudge_y = ifelse(fold_change > 0, -11, 11),
                  size = 3, colour = "black", box.padding = 1, point.padding = 0.8,
                  force = 4, force_pull = 0.5, min.segment.length = 0,
                  segment.size = 0.3, segment.color = "black", max.overlaps = Inf)