Ask Your Question
0

How can wildcards be preserved while expanding a dictionary in Snakemake?

asked 2023-02-13 11:00:00 +0000

qstack gravatar image

edit retag flag offensive close merge delete

1 Answer

Sort by ยป oldest newest most voted
3

answered 2022-11-04 23:00:00 +0000

scrum gravatar image

To preserve wildcards while expanding a dictionary in Snakemake, you can use the expand function with the wildcards parameter. Here's an example:

rule example:
    input:
        "input/{sample}/file_{num}.txt"
    output:
        "output/{sample}/result_{num}.out"
    params:
        param_dict=lambda wildcards: {"sample": wildcards.sample, "num": wildcards.num}
    shell:
        "command --input {input} --output {output} --params {params}"

expanded_inputs = expand("input/{sample}/file_{num}.txt", sample=["A", "B"], num=[1, 2])
expanded_outputs = expand("output/{sample}/result_{num}.out", wildcards, sample=["A", "B"], num=[1, 2])
params_dict = {"sample": "{sample}", "num": "{num}"}

rule example:
    input:
        input=expanded_inputs
    output:
        output=expanded_outputs
    params:
        params_dict=params_dict
    shell:
        "command --input {input} --output {output} --params {params}"

In this example, we use the expand function to generate the input and output file paths with wildcards. We then pass the wildcards object into the expand function for the output paths, which preserves the values of the wildcards. Finally, we use a params_dict dictionary to pass the values of the sample and num wildcards as parameters to the shell command.

edit flag offensive delete link more

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account. This space is reserved only for answers. If you would like to engage in a discussion, please instead post a comment under the question or an answer that you would like to discuss

Add Answer


Question Tools

Stats

Asked: 2023-02-13 11:00:00 +0000

Seen: 7 times

Last updated: Nov 04 '22