snakemake:访问字典值的正确语法

时间:2019-05-08 13:39:21

标签: python snakemake

这是我要做的事的一个例子:

mydictionary={
'apple': 'crunchy fruit',
'banana': 'mushy and yellow'
}

rule all:
    input:
        expand('{key}.txt', key=mydictionary.keys())

rule test:
    output: temp('{f}.txt')
    shell:
        """
        echo {mydictionary[wildcards.f]} > {output}
        cat {output}
        """

由于某些原因,我无法访问词典内容。我尝试使用双曲括号,但是文本文件的内容变为文字{mydictionary[wildcards.f]}(而我希望字典中相应条目的内容)。

1 个答案:

答案 0 :(得分:2)

我很确定方括号标记只能用变量值的字符串表示形式替换变量,但不支持方括号内的任何代码求值。也就是说,{mydictionary[wildcards.f]}将尝试查找一个字面名为"mydictionary[wildcards.f]"的变量。同样,{mydictionary}[{wildcards.f}]只会将字符串值粘贴在一起。因此,我认为您无法仅在shell部分中完成您想做的事情。相反,您可以在params部分中完成所需的操作:

rule test:
    output: temp('{f}.txt')
    params:
        value=lambda wcs: mydictionary[wcs.f]
    shell:
        """
        echo '{params.value}' > {output}
        cat {output}
        """