在for循环中使用ggplot函数绘制受试者轨迹图遇错求助

错误原因分析 #

1. 数据集拆分错误：原代码中gg_ID <- list(data_long$ID)仅将ID列的向量存入列表，而非按ID拆分的完整数据集，导致ggplot接收的是字符/数值向量而非数据框，触发报错。
2. 语法错误：geom_line()后多了一个冗余括号，且代码末尾的+符号破坏了ggplot的语法结构。
3. 变量名错误：geom_text中使用了不存在的score变量，实际数据中对应的变量是value。

完整解决方案代码 #

1. 宽格式转长格式（可选，替代手动创建T1-T104） #

如果你的宽格式数据已经存在，可通过tidyr快速转换为长格式：

library(tidyr)
library(dplyr)

# 转换宽格式到长格式，提取时间点为数值型
data_long <- df %>%
  pivot_longer(cols = starts_with("T"), 
               names_to = "time", 
               values_to = "value") %>%
  mutate(time = as.integer(stringr::str_remove(time, "T")))

2. 按ID拆分数据集并批量绘图+PDF保存 #

使用ggplot2绘图，patchwork排版，实现每10张图存为PDF的一页：

library(ggplot2)
library(patchwork)

# 按ID拆分数据集为列表
id_data_list <- split(data_long, data_long$ID)

# 配置PDF参数与每页图数量
plots_per_page <- 10
total_pages <- ceiling(length(id_data_list) / plots_per_page)

# 打开PDF输出设备
pdf("subject_trajectory_plots.pdf", width = 16, height = 10)

# 循环生成每页的图
for(page_num in 1:total_pages) {
  # 计算当前页的ID索引范围
  start_pos <- (page_num - 1) * plots_per_page + 1
  end_pos <- min(page_num * plots_per_page, length(id_data_list))
  current_data <- id_data_list[start_pos:end_pos]
  
  # 为当前页的每个ID生成轨迹图
  page_plots <- lapply(current_data, function(subject_data) {
    ggplot(subject_data, aes(x = time, y = value, group = ID)) +
      geom_line(na.rm = TRUE, size = 1.5, color = "#2c3e50") +
      geom_text(aes(label = value), vjust = -0.3, color = "black", size = 3.5) +
      xlab("Time Point") +
      ylab("Score") +
      ggtitle(paste("Subject ID:", unique(subject_data$ID))) +
      theme_bw() +
      theme(plot.title = element_text(size = 10),
            axis.title = element_text(size = 9),
            axis.text = element_text(size = 7))
  })
  
  # 组合当前页的图为2行5列布局
  combined_page <- wrap_plots(page_plots, ncol = 5, nrow = 2)
  print(combined_page)
}

# 关闭PDF设备，完成保存
dev.off()

关键优化说明 #

• 单个受试者的图无需按ID着色（每个子集仅含一个ID），去掉col = ID可避免多余图例，简化视觉效果。
• 通过wrap_plots实现自动排版，可灵活调整ncol/nrow参数改变每页图的排列方式。
• 调整theme中的字体大小，确保多图排版时文字不会重叠。

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