R语言多区域δ¹⁸O深度散点图并排绘制无显示问题排查及优化

1. 数据预处理（合并+标准化） #

先合并两个数据源，统一列名并筛选所需字段，避免重复处理两个数据库的冗余代码：

# 加载所需工具包
library(tidyverse)

# 假设两个数据库的列名分别为：
# Database1: Region, MBSL, d18O_Core
# Database2: Region, MBSL, d18O_Plankton

# 合并并整理数据
combined_data <- bind_rows(
  Database1 %>% select(Region, MBSL, d18O = d18O_Core) %>% mutate(Source = "Core Tops"),
  Database2 %>% select(Region, MBSL, d18O = d18O_Plankton) %>% mutate(Source = "Plankton Tows")
) %>%
  # 筛选目标区域
  filter(Region %in% c("Labrador", "Greenland", "Iceland")) %>%
  # 去除深度或δ¹⁸O的缺失值（避免空图）
  drop_na(MBSL, d18O)

2. 绘制横向并排散点图（简化代码+解决无输出） #

用ggplot2的分面功能实现三图横向排列，无需重复编写三次绘图代码；同时确保代码能正确输出图形：

# 生成绘图对象
p <- ggplot(combined_data, aes(x = d18O, y = MBSL, color = Source)) +
  geom_point(alpha = 0.7) +
  # 反转y轴（深度数值越大代表越深，反转后更符合认知）
  scale_y_reverse() +
  # 设置分面为横向1行排列
  facet_wrap(~Region, nrow = 1) +
  # 添加标签和主题
  labs(
    x = expression(delta^{18}*"O (‰)"),
    y = "Depth [MBSL, m]",
    color = "Data Source"
  ) +
  theme_minimal() +
  theme(
    strip.text = element_text(size = 12, face = "bold"),
    legend.position = "bottom"
  )

# 强制输出图形（脚本模式下必须添加，否则可能无输出）
print(p)

3. 无图形输出的常见排查点 #

• 未加载ggplot2包：确保开头运行library(tidyverse)或library(ggplot2)
• 数据为空：检查combined_data是否有有效数据，可运行dim(combined_data)查看行数
• 脚本模式未调用print()：在R脚本中，ggplot对象需要用print()触发绘图；RStudio交互模式可直接运行对象名
• 坐标轴范围异常：若数据值范围极端，可手动设置xlim()/ylim()调整

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