在Nextflow脚本中运行R命令时未生成输出文件求助

问题背景 #

运行nextflow run main.nf命令后，预期的output.rds文件未生成，相关文件内容如下：

main.nf文件 #

#!/usr/bin/env nextflow
nextflow.enable.dsl=2 

process my_script {

    """
    Rscript script.R
    """
}

workflow {
  my_script 
}

nextflow.config文件 #

process {
   executor = 'k8s'
   container = 'rocker/r-ver:4.1.3'
}

script.R文件 #

FUN <- readRDS("function.rds");
input = readRDS("input.rds");
output = FUN(
singleCell_data_input = input[[1]], savePath = input[[2]], tmpDirGC = input[[3]]
);
saveRDS(output, "output.rds")

排查与解决步骤 #

• 1. 确保输入文件被正确传入容器
  R脚本依赖function.rds和input.rds，但Kubernetes执行器不会自动挂载本地文件到容器。需在进程中声明输入：

  process my_script {
      input:
      path 'function.rds'
      path 'input.rds'
  
      """
      Rscript script.R
      """
  }
  

  同时保证这两个文件位于Nextflow运行的工作目录下，或使用绝对路径指定。

• 2. 明确声明输出文件
  Nextflow不会自动收集进程内生成的文件，需显式声明输出，否则文件会随Pod销毁丢失：

  process my_script {
      input:
      path 'function.rds'
      path 'input.rds'
  
      output:
      path 'output.rds'
  
      """
      Rscript script.R
      """
  }
  

  声明后，Nextflow会将文件复制到本地工作目录（默认./work下的对应进程子目录）。若需指定输出路径，可添加publishDir配置：

  process my_script {
      input:
      path 'function.rds'
      path 'input.rds'
  
      output:
      path 'output.rds'
  
      publishDir './results', mode: 'copy'
  
      """
      Rscript script.R
      """
  }
  

• 3. 添加错误捕获排查脚本问题
  修改R脚本添加错误日志，确认执行是否成功：

  tryCatch({
      FUN <- readRDS("function.rds")
      input <- readRDS("input.rds")
      output <- FUN(
          singleCell_data_input = input[[1]], 
          savePath = input[[2]], 
          tmpDirGC = input[[3]]
      )
      saveRDS(output, "output.rds")
      cat("脚本执行成功，已生成output.rds\n")
  }, error = function(e) {
      cat("脚本执行失败：", e$message, "\n")
      quit(status = 1)
  })
  

  运行流程时启用调试日志：nextflow run main.nf -log debug，或查看Kubernetes Pod日志，定位具体报错。

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