如何用for循环生成Snakemake规则？执行异常及方案咨询

这是Python闭包延迟绑定导致的典型问题：循环里定义的run代码块并没有立刻捕获thiskey和thisperson的当前值，而是在规则实际执行时才读取这两个变量的状态——此时循环已经结束，变量停留在最后一次迭代的person2和Bruce，所以所有规则都复用了这组值。

解决思路1：用函数封装固定变量值 #

通过辅助函数把循环变量作为参数传入，利用函数参数的即时绑定特性，强制捕获当前迭代的变量值：

config={'person1': 'Amy', 'person2': 'Bruce'}

def make_person_rule(key, name):
    rule:
        name: f'{name}s_rule'
        output: f'{key}.txt'
        run:
            print(f'{key} is called {name}!')

for thiskey, thisperson in config.items():
    make_person_rule(thiskey, thisperson)

rule people:
    input: [f'{k}.txt' for k in config.keys()]

解决思路2：使用Snakemake参数化规则（更符合工具设计风格） #

Snakemake本身推荐用通配符+参数化规则替代手动循环生成规则，既简洁又能从根源避免变量绑定问题：

config={'person1': 'Amy', 'person2': 'Bruce'}

rule person_rule:
    output: '{key}.txt'
    params:
        name=lambda wildcards: config[wildcards.key]
    run:
        print(f'{wildcards.key} is called {params.name}!')

rule people:
    input: expand('{key}.txt', key=config.keys())

原方式的合理性 #

你用循环生成规则的思路本身可行，但忽略了Python闭包的延迟绑定特性，导致变量值未被正确捕获。上面两种方案都能解决问题，其中参数化规则的方式是Snakemake的标准用法，更推荐采用。

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