如何在ggplot中让散点图折线按深度递减顺序正确连接点

问题根源 #

当前代码存在3个核心错误导致折线连接逻辑混乱：

• 错误地将Depth转为因子类型，ggplot会按因子默认水平顺序而非数值逻辑连接点
• 无效列赋值：df$Phase <- as.factor(df$Var)和df$Temp <- df$val，原数据集无Var/val列，直接覆盖了正确的原始数据
• 数据未按分组排序，折线连接顺序完全随机

修正方案与完整代码 #

步骤1：修正数据预处理 #

删除无效赋值，对数据按Phase分组后，按Depth递增（对应深度从浅到深）排序，确保折线连接顺序符合深度变化逻辑。

步骤2：调整绘图逻辑 #

保留Depth的数值类型，用scale_y_reverse()实现深度从浅到深的y轴展示，同时匹配图例参数与实际分组数量。

# 加载依赖包
library(ggplot2)
library(dplyr)

# 原始数据集
df <-  structure(list(X = 1:43, Phase = c("Phase1a", "Phase1a", "Phase1a", 
                                   "Phase1a", "Phase1a", "Phase1a", "Phase1a", "Phase1a", "Phase1a", 
                                   "Phase1a", "Phase1a", "Phase1a", "Phase1a", "Phase1a", "Phase1a", 
                                   "Phase1b", "Phase1b", "Phase1b", "Phase1b", "Phase1b", "Phase1b", 
                                   "Phase1b", "Phase1b", "Phase1b", "Phase1b", "Phase1b", "Phase1b", 
                                   "Phase1b", "Phase1b", "Phase2", "Phase2", "Phase2", "Phase2", 
                                   "Phase2", "Phase2", "Phase2", "Phase2", "Phase2", "Phase2", "Phase2", 
                                   "Phase2", "Phase2", "Phase2"), Depth = c(0.3, 2, 1, 3, 4, 5, 
                                                                            6, 7, 8, 9, 10, 11, 12, 13, 14, 0.3, 2, 1, 3, 4, 5, 6, 7, 8, 
                                                                            9, 10, 11, 12, 13, 0.3, 2, 1, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 
                                                                            13), Temp = c(14.12, 14.07105, 14.1025, 13.8125, 13.56, 13.47, 
                                                                                          13.3325, 13.0775, 12.9225, 12.425, 11.025, 9.955, 10.6333333333333, 
                                                                                          9.355, 9.23, 17.7183333333333, 17.84, 17.93, 17.39, 17.065, 16.92, 
                                                                                          16.795, 16.585, 14.92, 13.53, 12.515, 11.685, 11.32, 12.11, 11.75, 
                                                                                          11.984, 12.044, 11.97, 11.954, 11.948, 11.944, 11.858, 11.772, 
                                                                                          11.786, 11.774, 15.1, 15.1, 15.12), DO = c(9.8, 9.4475, 9.5425, 
                                                                                                                                     9.505, 9.76, 9.39, 9.345, 9.25, 9.615, 9.315, 9.2475, 9.43, 7.76, 
                                                                                                                                     7.34, 6.61, 8.975, 8.12, 8.085, 8.21, 8.235, 8.14, 8.235, 8.145, 
                                                                                                                                     8.035, 7.99, 7.595, 6.455, 5.78, 6.74, 10.6846153846154, 10.6, 
                                                                                                                                     10.57, 10.582, 10.486, 10.502, 10.54, 10.656, 10.7, 10.52, 10.616, 
                                                                                                                                     8.41, 8.11, 7.88), Cond = c(36.25, 36.25, 36.25, 36, 35.6666666666667, 
                                                                                                                                                                 35.725, 35.725, 35.225, 34.95, 34.675, 36.1, 32, 45.3333333333333, 
                                                                                                                                                                 32.5, 31, 42.4333333333333, 43.5, 44, 43.5, 42.5, 42, 42.5, 42, 
                                                                                                                                                                 40.5, 38.5, 38.5, 49.5, 41.5, 41, 37.9, 38.24, 38.24, 38.24, 
                                                                                                                                                                 38.24, 38.22, 38.22, 37.98, 38.16, 38.56, 38.92, 42, 42, 42), 
               SPCond = c(45.9, 45.65, 45.875, 45.4, 45.2, 45.4, 45.4, 45.4, 
                          45.425, 45.35, 49.1, 44.5, 61, 46, 44, 49.3333333333333, 
                          51, 51, 51, 50, 50, 50.5, 50.5, 50.5, 50.5, 51, 65.5, 56.5, 
                          54, 50.7692307692308, 50.9, 50.9, 50.7, 50.9, 50.9, 50.7, 
                          50.7, 50.9, 51.7, 52.1, 52, 52, 52)), class = "data.frame", row.names = c(NA, 
                                                                                                    -43L))

# 按分组排序数据：每个Phase内按Depth从小到大排列
df_sorted <- df %>%
  group_by(Phase) %>%
  arrange(Depth, .by_group = TRUE) %>%
  ungroup()

# 设定Phase因子顺序（可选，保证图例顺序符合预期）
df_sorted$Phase <- factor(df_sorted$Phase, levels = c("Phase1a", "Phase1b", "Phase2"))

# 绘图
ggplot(data = df_sorted, aes(x = Temp, y = Depth, group = Phase)) +
  geom_line(colour = "black") +
  geom_point(aes(shape = Phase, color = Phase, fill = Phase), size = 2) +
  # 匹配3个Phase分组的图例参数
  scale_shape_manual(values = c(24, 21, 22)) +
  scale_color_manual(values = c('black', 'black', 'black')) +
  scale_fill_manual(values = c('blue', 'green', 'red')) +
  theme_bw() +
  theme(
    legend.title = element_blank(),
    panel.border = element_rect(colour = "black", fill = NA, size = 1),
    legend.position = 'right',
    legend.background = element_blank()
  ) +
  geom_vline(xintercept = 4, linetype = "dotted") +
  labs(x = "Temperature (°C)", y = "Depth (m)") +
  scale_x_continuous(position = "top") +
  # 反转y轴，实现深度浅的在上、深的在下的展示
  scale_y_reverse()

核心修正说明 #

• 数据排序：通过分组排序确保每个Phase的点按深度从浅到深排列，折线连接顺序完全符合深度变化逻辑
• 变量类型：保留Depth的数值属性，避免因子类型导致的非逻辑连接
• y轴优化：用scale_y_reverse()替代离散型反转，更适配数值型深度变量，同时实现温度随深度增加（y轴向下）而下降的效果
• 图例适配：调整图例参数数量与实际分组数一致，消除不必要的警告

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