Nextflow流程向Rscript注入变量失败，生成空Excel问题排查

1. Channel传递目录的方式错误 #

你用Channel.fromPath(params.bamstat_outdir)时，默认会把目录下的每个文件作为独立的Channel元素，导致collect流程为每个文件单独执行一次，而非接收整个目录。这会让R脚本里的文件查找逻辑失效。

修正方法：给fromPath添加type: 'dir'参数，将整个目录作为单个元素传入Channel：

stat_qc = Channel.fromPath(params.bamstat_outdir, type: 'dir', checkIfExists: true)

2. R脚本的核心问题 #

（1）正则表达式写法错误

R的dir()函数pattern参数是正则表达式，你写的"*.out"会被错误解析，需用"\\.out$"来匹配所有以.out结尾的文件。

（2）缺少读取文件与写入Excel的逻辑

你仅获取了文件名列表，但未读取文件内容，也没有调用write_xlsx()生成有效内容，这是Excel为空的直接原因。

（3）冗余的环境变量传递

没必要通过export传递参数，直接用Nextflow的!{}插值把目录路径传入R脚本更直接可靠。

修正后的collect.nf脚本：

process collect {
    publishDir "${params.qc_outdir}", mode: "copy"
    
    input:
    path stat_qc  # 现在接收的是整个目录
    
    output:
    path "mapping_stats.xlsx", emit: map_stat
    
    script:
    """
    Rscript -e '
    library(writexl)

    # 直接通过Nextflow插值传入目录路径
    stat_dir <- "!{stat_qc}"

    # 匹配所有.out文件，获取完整路径
    out_files <- list.files(stat_dir, pattern = "\\\\.out$", full.names = TRUE)

    # 无文件时提前终止，避免生成空Excel
    if(length(out_files) == 0) {
        stop("No .out files found in directory: ", stat_dir)
    }

    # 读取每个.out文件内容（示例按行读取，可根据你的文件格式调整，比如用read.table）
    stat_data <- lapply(out_files, function(file) {
        file_content <- readLines(file)
        data.frame(
            Filename = basename(file),
            Content = file_content,
            stringsAsFactors = FALSE
        )
    })

    # 合并所有数据
    all_stats <- do.call(rbind, stat_data)

    # 写入Excel文件
    write_xlsx(all_stats, "mapping_stats.xlsx")
    '
    """
}

3. 启动脚本的冗余代码 #

你的bash脚本中bamstat_file=ls $bamstat_outdir/*``这行未被使用，可直接删除，不影响流程运行。

验证提示 #

• 确认params.bamstat_outdir目录下确实存在.out文件
• 运行流程时可添加-with-trace参数查看日志，排查R脚本是否有报错
• 如果你的.out是特定格式（如制表符分隔表格），需将readLines替换为对应读取函数（如read.delim）

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