Snakemake调用R脚本时write.table报错：invalid connection

你遇到的Error in isOpen(file, "w") : invalid connection错误，核心原因及修复方案如下：

核心错误原因 #

1. 输出路径类型错误：snakemake@output返回的是列表对象，而非单个字符串路径。直接将列表传给write.table的file参数，R无法识别为有效文件路径，引发无效连接错误。
2. 路径格式潜在问题：规则中输出路径的双斜杠//虽部分系统兼容，但可能导致路径解析异常。
3. merge逻辑隐藏坑：读取文件时用row.names=1将gene_id设为行名，后续merge(..., by = "gene_id")会找不到该列，导致合并失败。

修复方案 #

1. 修正输出路径获取方式 #

在R脚本中，将列表转为单个路径字符串：

out_file_path <- snakemake@output[[1]]

2. 统一路径格式（可选） #

将Snakemake规则中的输出路径双斜杠改为单斜杠：

rule prepare_files:
    input: 
        x = "path_to_x/quant.sf", 
        y = "path_to_Y/quant.sf"
    output: 
        "out/Bigquantif.txt"   
    script:
        "R/manips.R"

3. 确保输出目录存在（可选） #

添加代码自动创建输出目录，避免手动创建遗漏：

dir.create(dirname(out_file_path), recursive = TRUE, showWarnings = FALSE)

4. 修正merge逻辑 #

基于行名合并数据，适配row.names=1的读取方式：

merged_data_final <- merge(x_fi, y_fi, by = "row.names", all = TRUE)
rownames(merged_data_final) <- merged_data_final$Row.names
merged_data_final$Row.names <- NULL

完整修正后的manips.R #

out_file_path <- snakemake@output[[1]]
x_path <- snakemake@input$x
y_path <- snakemake@input$y

# 确保输出目录存在
dir.create(dirname(out_file_path), recursive = TRUE, showWarnings = FALSE)

x_fi <- read.table(x_path, header = TRUE, sep = "\t", row.names = 1)
y_fi <- read.table(y_path, header = TRUE, sep = "\t", row.names = 1)

# 基于行名合并数据
merged_data_final <- merge(x_fi, y_fi, by = "row.names", all = TRUE)
rownames(merged_data_final) <- merged_data_final$Row.names
merged_data_final$Row.names <- NULL

write.table(x = merged_data_final, file = out_file_path, sep = "\t", quote = FALSE, row.names = TRUE)

内容的提问来源于stack exchange，提问作者Dima's