R循环遍历22条染色体文件时read.table执行报错如何解决？

报错原因 #

• 路径拼接语法错误：R中无法直接通过逗号拆分多段字符串拼接文件路径，逗号为函数参数分隔符，拆分的路径片段会被识别为read.table的多个独立参数，导致参数逻辑混乱。
• 参数分隔符缺失：read.table中文件路径的最后一段字符串".QC.het"和后续参数header=T之间没有加逗号分隔，R无法识别两个独立参数，触发语法报错。

另外原代码中write.table也存在相同的路径拼接错误，执行到对应步骤也会触发同类报错。

修复后的完整代码 #

使用paste0()函数完成无分隔符的字符串拼接，修正参数分隔问题：

for(i in 1:22) {
  # 拼接输入文件路径
  input_file <- paste0("~/path/to/files/chromosome", i, ".QC.het")
  dat <- read.table(input_file, header = T) 
  # 计算F值均值和标准差
  m <- mean(dat$F)
  s <- sd(dat$F)
  # 筛选3倍标准差内的有效样本
  valid <- subset(dat, F <= m+3*s & F >= m-3*s)
  # 拼接输出文件路径
  output_file <- paste0("~/path/to/files/chromosome", i, ".valid.sample")
  write.table(valid[,c(1,2)], output_file, quote = F, row.names = F)
}

代码说明 #

• paste0()是R内置的无分隔符字符串拼接函数，多个待拼接内容直接放入函数参数中用逗号分隔即可
• 原代码中id <- i的定义无实际作用，直接使用循环变量i拼接路径即可
• 所有函数的不同参数之间必须用逗号分隔，避免语法识别错误

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