R遍历受试者子目录处理AD矩阵、宽转长并保存的for循环编写问题

原有代码核心问题 #

• 路径拼接错误：R不能直接用/拼接字符串路径，需要用file.path()函数处理跨系统兼容的路径
• 目录读取逻辑错误：list.dirs()的输入路径错误，没有基于project1的根目录定位受试者文件夹
• 循环变量逻辑错误：循环中把路径变量subj.AD当成了数字索引遍历，无法匹配实际的受试者路径
• 读写文件名错误：硬编码读取matrix.csv，没有指向每个受试者下AD/AD.csv；输出文件名没有动态提取受试者ID，路径也完全错误
• 语法错误：write.csv行末尾引号未闭合，会直接报语法错误

完整可运行代码 #

# 依赖包安装（已安装可注释）
# install.packages("reshape2")
library(reshape2)

# 定义project1的根路径，根据你本地的实际路径修改
proj_root <- "group/S/project1/"

# 读取所有受试者文件夹路径，过滤只保留Sub-开头的文件夹避免读入无关文件
subject_dirs <- list.dirs(path = proj_root, full.names = TRUE, recursive = FALSE)
subject_dirs <- subject_dirs[grepl("Sub-\\d+", basename(subject_dirs))]

# 遍历每个受试者文件夹处理
for (subj_dir in subject_dirs) {
  # 提取受试者ID（如Sub-001）
  subj_id <- basename(subj_dir)
  # 拼接AD.csv的完整读取路径
  ad_csv_path <- file.path(subj_dir, "AD", "AD.csv")
  # 读取矩阵文件，分隔符、header参数根据你的AD.csv实际格式调整
  df <- read.csv(ad_csv_path, sep = ",", header = FALSE, stringsAsFactors = FALSE)
  # 矩阵转换、保留上三角、宽转长处理
  df_matrix <- as.matrix(df)
  df_matrix[lower.tri(df_matrix)] <- NA
  df_melt <- melt(df_matrix, na.rm = TRUE)
  # 拼接输出路径，存到受试者根目录，命名符合要求
  output_path <- file.path(subj_dir, paste0(subj_id, "_AD_reshaped.csv"))
  # 输出文件，关闭多余行号输出
  write.csv(df_melt, file = output_path, row.names = FALSE)
}

注意事项 #

• 如果你的AD.csv第一行是列名，把read.csv的header = FALSE改成header = TRUE
• 如果AD.csv的分隔符不是逗号是制表符，把sep = ","改成sep = "\t"
• 运行前先确认proj_root的路径是你本地project1文件夹的完整绝对路径，避免相对路径定位错误

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