ggplot添加多条布朗运动曲线异常：仅显示两条曲线的问题排查

一、检查循环遍历列的逻辑是否正确 #

你的data.frame结构是奇数列=时间，偶数列=对应样本数值，比如5个样本对应10列：列1/2为样本1，列3/4为样本2，…，列9/10为样本5。如果循环时步长或索引范围错误，会导致只处理部分样本。

错误示例（逻辑缺陷） #

Plot_a_Sample_of_Brownian_Motions <- function(df) {
  # 错误：手动限制了遍历范围，仅处理前2个样本
  for (i in seq(1, 4, 2)) {
    time_col <- df[, i]
    value_col <- df[, i+1]
    Plot_a_brownian_motion(time_col, value_col)
  }
}

修正后的遍历逻辑 #

Plot_a_Sample_of_Brownian_Motions <- function(df) {
  # 获取所有奇数列的索引（对应各样本的时间列）
  time_indices <- seq(1, ncol(df), by = 2)
  # 遍历每个样本的时间和数值列
  for (i in time_indices) {
    time_col <- df[, i]
    value_col <- df[, i + 1]
    Plot_a_brownian_motion(time_col, value_col)
  }
}

二、确保绘图时不覆盖已有曲线 #

如果你的Plot_a_brownian_motion基于R基础绘图系统，默认每次调用都会新建画布，覆盖之前的曲线。需要添加叠加参数，并统一坐标轴范围避免曲线被裁剪。

改造单条曲线绘图函数 #

Plot_a_brownian_motion <- function(time, value, add = FALSE) {
  if (add) {
    # 叠加到已有画布上，用随机颜色区分样本
    lines(time, value, col = sample(colors(), 1))
  } else {
    # 初始化画布，设置覆盖所有样本的坐标轴范围
    plot(time, value, type = "l", xlab = "时间", ylab = "数值", 
         xlim = range(time), ylim = range(value))
  }
}

改造批量绘图函数（先初始化画布，再叠加所有曲线） #

Plot_a_Sample_of_Brownian_Motions <- function(df) {
  time_indices <- seq(1, ncol(df), by = 2)
  # 提取所有样本的时间和数值，计算全局坐标轴范围
  all_time <- unlist(df[, time_indices])
  all_values <- unlist(df[, time_indices + 1])
  
  # 绘制第一个样本作为基础画布
  first_time <- df[, time_indices[1]]
  first_value <- df[, time_indices[1] + 1]
  plot(first_time, first_value, type = "l", xlab = "时间", ylab = "数值",
       xlim = range(all_time), ylim = range(all_values))
  
  # 叠加剩余所有样本的曲线
  for (i in time_indices[-1]) {
    time_col <- df[, i]
    value_col <- df[, i + 1]
    lines(time_col, value_col, col = sample(colors(), 1))
  }
}

如果使用ggplot2，可以更简洁地实现批量绘图（推荐）：

library(ggplot2)
library(tidyr)

Plot_a_Sample_of_Brownian_Motions <- function(df) {
  # 将宽格式data.frame转换为长格式，方便ggplot批量绘图
  long_df <- data.frame(
    time = unlist(df[, seq(1, ncol(df), 2)]),
    value = unlist(df[, seq(2, ncol(df), 2)]),
    sample = rep(paste0("样本", 1:(ncol(df)/2)), each = nrow(df))
  )
  
  ggplot(long_df, aes(x = time, y = value, color = sample)) +
    geom_line() +
    labs(x = "时间", y = "数值", title = "布朗运动样本集")
}

三、验证输入数据的完整性 #

检查你的data.frame是否真的包含完整的5个样本数据：

# 检查列数是否为偶数（每个样本对应2列）
ncol(df) %% 2 == 0
# 查看前几行数据，确认所有列都有有效数值
head(df)
# 检查每个数值列是否存在全NA情况（全NA列无法绘制曲线）
sapply(df[, seq(2, ncol(df), 2)], function(x) sum(is.na(x)))

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