Snakefile中能否使用Dynamic Branching/Plumbing实现条件分支流程？

1. 第一种写法首先存在笔误，else分支也写了rules.B，如果修正为rules.C后仍同时运行两个规则，大概率是判断条件存在布尔值和字符串混用的问题：你在下游D规则的input里用config["deblur"] == "True"判断字符串，而分支判断用config["deblur"] == True判断布尔值，二者不匹配会导致判断失效，两个规则都被加入依赖队列。
2. 第二种写法触发语法错误是因为你将if-else分支语句直接写在了规则的属性块（比如input块）内部，规则的属性块仅支持键值对和表达式，不支持直接写分支语句。

以下三种方案都适用于config["deblur"]为运行前就确定的静态配置的场景：

方案1：目标规则直接指定条件依赖（最简洁） #

在你的最终目标规则（比如all规则）的input里直接用三元表达式指定上游分支即可，Snakemake会只解析返回的对应规则作为依赖：

rule all:
    input:
        # 其余全局依赖不变
        rules.D.output,
        rules.B if config["deblur"] else rules.C

# 其余规则定义不变，你已经写好的D规则input三元适配可以保留

方案2：全局定义分支变量（易维护） #

如果分支判断逻辑需要在多个规则里复用，可以在Snakefile全局提前判断分支，后续所有位置直接复用变量即可：

# 全局提前判断分支，统一判断逻辑避免不匹配
if config["deblur"]:
    branch_rule = rules.B
    branch_table = rules.qiime_deblur.output.qiime_deblur_table
else:
    branch_rule = rules.C
    branch_table = rules.qiime_denoise.output.qiime_denoise_table

rule all:
    input: rules.D.output, branch_rule

rule D:
    input:
        qiime_feature_table_input = branch_table
    # 其余属性不变

方案3：规则优先级匹配（适合分支输出结构一致的场景） #

如果B、C两个规则的输出结构完全相同，可以给两个规则设置相同的输出路径，再根据config值设置规则优先级，Snakefile会自动选择优先级高的规则生成中间文件，不需要在每个规则里写判断：

rule B:
    input: rules.A.output
    output: "middle/feature_table.qza"
    # 其余属性不变

rule C:
    input: rules.A.output
    output: "middle/feature_table.qza"
    # 其余属性不变

# 根据config设置规则优先级
if config["deblur"]:
    ruleorder: B > C
else:
    ruleorder: C > B

rule D:
    input:
        qiime_feature_table_input = "middle/feature_table.qza"
    # 其余属性不变

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