如何在Snakemake的shell代码块中访问配置变量？

嘿，这个问题我之前也碰到过！在Snakemake的shell指令里直接调用config确实有点小坑，但其实有几种干净的解决办法，完全不用妥协那些不必要的副作用：

方法一：直接在Shell指令中嵌入Python表达式引用Config #

Snakemake的shell:块支持直接嵌入Python代码（用{}包裹），你可以直接通过字典索引或者属性访问的方式调用config里的键，比如：

rule my_rule:
    input: "data/input.txt"
    output: "results/output.txt"
    shell:
        """
        # 直接引用顶级配置键
        echo {config["quux"]} > {output}
        
        # 如果是嵌套配置，比如config里有个section叫settings，里面有key叫quux
        echo {config.settings.quux} >> {output}
        """

这种方法最直接，完全不需要额外的中间变量，也不会产生任何元数据跟踪的副作用。

方法二：用Params但跳过跟踪 #

如果你习惯用params来传递变量，但又不想让参数被保存到.snakemake/params_tracking里，可以给rule加上params_tracking: False的设置，这样就能避免不必要的元数据存储：

rule my_rule:
    input: "data/input.txt"
    output: "results/output.txt"
    params:
        quux = config["quux"]
    params_tracking: False  # 关键：禁用这个rule的参数跟踪
    shell:
        """
        echo {params.quux} > {output}
        """

这种方法适合那些已经习惯了用params传递变量的场景，同时解决了副作用问题。

方法三：用Run块调用Shell（简洁版） #

你提到了用run:替代shell:，其实完全不用写复杂的逻辑，直接用Python的f-string调用shell就行，这样能自由访问所有Snakemake变量（包括config）：

rule my_rule:
    input: "data/input.txt"
    output: "results/output.txt"
    run:
        # 用f-string拼接shell命令，直接访问config
        shell(f"echo {config['quux']} > {output}")

这种方式灵活度很高，如果你需要在执行shell前加一些简单的Python逻辑，也能轻松扩展。

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