Question

我对以下规则有疑问，看起来外卡正在改变一半。

> maker_flow={'boot1':'init', 'boot2':'boot1', 'final':'boot2'} #this defines the bootstraping workflow {now:last}
> rule maker:
>     input: snap= 'another_file_string_with_same_'.format(), gff=lambda wc: 'maker/{lrun}/{name}_{lrun}.maker.output/{name}_{lrun}.all.gff'.format(lrun=maker_flow[wc.run], name=config['proj']['name'])
>     output:   master='another_file_string_with_same'.format(), gff='maker/{{run}}/{name}_{{run}.maker.output/{name}_{{run}}.all.gff'.format(name=config['proj']['name'])

当我注释掉gff输入时，代码将会运行。否则我得到：

rule maker:
    input: maker/init/snap/hmm/test_10l_init.snap_hmm, maker/init/test_10l_init.maker.output/test_10l_init.all.gff
    output: maker/boot1/test_10l_boot1.maker.output/test_10l_boot1_master_datastore_index.log, maker/boot1/test_10l_boot1.maker.output/test_10l_boot1.all.gff
    log: /gpfs/projects/bgmp/shared/lizards/test_data/sandbox/snake/logs/test_10l_boot1_maker.log
    jobid: 4
    wildcards: run=boot1

InputFunctionException in line 95 of /gpfs/projects/bgmp/shared/lizards/test_data/sandbox/snake/Snakefile:
KeyError: 'init'
Wildcards:
run=init

这没有多大意义，因为gff文件的输入似乎已成功解决。而且通配符似乎已经改变了。

还有另一条规则可以创建＆＃39; init＆＃39;文件

rule maker_init:
    output:
        'maker/{run}/{name}_{run}.maker.output/{name}_{run}_master_datastore_index.log'.format(name=config['proj']['name'], run='init'), gff=protected('maker/init/{name}_{run}.maker.output/{name}_init.all.gff'.format(name=config['proj']['name'], run='init'))

Answer 1

在评估run = boot1的规则时，snakemake发现它需要文件“maker / init / test_10l_init.maker.output / test_10l_init.all.gff”。这个文件应该已经存在吗？如果是这样，似乎snakemake没有看到它，因为相反它看到它可以使用制造商规则模式来创建这个文件（所以现在它用run = init来评估规则，因此通配符更改）但是当它尝试时要在maker_flow中查找键'init'，它会得到一个KeyError。

snakemake：通配符正在改变输入函数

1 个答案: