如何在Snakefile中指定工作目录外的输入路径？

这里有几种实用方法帮你摆脱繁琐的相对路径：

• 直接使用绝对路径
  直接写出数据目录的完整绝对路径，比如：

  rule process_data:
      input: "/home/yourname/datasets/sample.fastq"
      output: "processed/sample_clean.fastq"
      shell: "trimmomatic SE {input} {output}"
  

  优点是一目了然，不用嵌套..；缺点是换机器或路径变动时需要修改代码。

• 定义全局路径变量
  在Snakefile开头定义一个全局变量存储数据目录路径，后续直接引用：

  # 开头定义路径
  DATA_DIR = "/home/yourname/datasets"
  
  rule process_data:
      input: f"{DATA_DIR}/sample.fastq"
      output: "processed/sample_clean.fastq"
      shell: "trimmomatic SE {input} {output}"
  

  后续要改路径只需要修改一处变量，维护更方便。

• 用配置文件管理路径
  创建config.yaml配置文件，写入数据目录路径：

  data_dir: "/home/yourname/datasets"
  

  在Snakefile中加载配置并引用：

  configfile: "config.yaml"
  
  rule process_data:
      input: f"{config['data_dir']}/sample.fastq"
      output: "processed/sample_clean.fastq"
      shell: "trimmomatic SE {input} {output}"
  

  这种方式适合多环境切换或协作场景，修改配置文件即可，无需改动Snakefile核心逻辑。

• 命令行临时传递路径
  运行Snakemake时通过--config参数临时指定数据目录，比如：

  snakemake --config data_dir=/home/yourname/datasets
  

  Snakefile中同样通过config引用：

  rule process_data:
      input: f"{config['data_dir']}/sample.fastq"
      output: "processed/sample_clean.fastq"
      shell: "trimmomatic SE {input} {output}"
  

  适合临时更换数据目录的场景，不用修改任何文件。

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