我试图在一个变量中同时将所有路径提供给python
中的snakemake
脚本,如下所示:
rule neo4j:
input:
script = 'python/neo4j.py',
path_to_cl = 'results/clusters/umap/{sample}_umap_clusters.csv',
path_to_umap = 'results/umap/{sample}_umap.csv',
path_to_mtx = 'data_files/normalized/{sample}.csv'
output: 'results/neo4j/{sample}/cells.csv', 'results/neo4j/{sample}/genes.csv',
'results/neo4j/{sample}/cl_nodes.csv', 'results/neo4j/{sample}/cl_contains.csv',
'results/neo4j/{sample}/cl_isin.csv', 'results/neo4j/{sample}/expr_by.csv',
'results/neo4j/{sample}/expr_ess.csv'
shell:
"python {input.script} -path_to_cl {input.path_to_cl} -path_to_umap {input.path_to_umap} -path_to_mtx {input.path_to_mtx} -output {output}"
当我在output
脚本中访问python
参数时,它只看到第一条路径:'results/neo4j/{sample}/cells.csv'
。我也试过命名每个路径,但它没有解决问题。如何将规则的output
中的所有路径作为数组或字典传递,以便以后在python
中访问它们?
如果我正确理解你的问题,你的问题是neo4j.py
脚本不接受-output
参数的多个文件:shell命令可能以完整的文件列表结束(请查看-p
的snakemake
选项),但仅限于脚本会考虑第一个问题。
如果确实如此,可能更简洁的方法是修改neo4j.py
脚本的接口,使其为每个输出文件使用一个参数。
然后,您将修改您的规则,如下所示:
rule neo4j:
input:
script = 'python/neo4j.py',
path_to_cl = 'results/clusters/umap/{sample}_umap_clusters.csv',
path_to_umap = 'results/umap/{sample}_umap.csv',
path_to_mtx = 'data_files/normalized/{sample}.csv'
output:
cells = 'results/neo4j/{sample}/cells.csv',
genes = 'results/neo4j/{sample}/genes.csv',
nodes = 'results/neo4j/{sample}/cl_nodes.csv',
contains = 'results/neo4j/{sample}/cl_contains.csv',
isin = 'results/neo4j/{sample}/cl_isin.csv',
by = 'results/neo4j/{sample}/expr_by.csv',
ess = 'results/neo4j/{sample}/expr_ess.csv'
shell:
"""
python {input.script} \\
--path_to_cl {input.path_to_cl} \\
--path_to_umap {input.path_to_umap} \\
--path_to_mtx {input.path_to_mtx} \\
--cells {output.cells} \\
--genes {output.genes} \\
--nodes {output.nodes} \\
--contains {output.contains} \\
--isin {output.isin} \\
--by {output.by} \\
--ess {output.ess}
"""
一些可能有用的python模块来设置脚本的界面:
如果您不想将每个输入文件作为单个参数传递,则可以简单地传递输出目录,并让脚本从这个单个参数构建输出路径。鉴于您想要的文件名,这似乎是可能的:
rule neo4j:
input:
script = 'python/neo4j.py',
path_to_cl = 'results/clusters/umap/{sample}_umap_clusters.csv',
path_to_umap = 'results/umap/{sample}_umap.csv',
path_to_mtx = 'data_files/normalized/{sample}.csv'
output:
'results/neo4j/{sample}/cells.csv',
'results/neo4j/{sample}/genes.csv',
'results/neo4j/{sample}/cl_nodes.csv',
'results/neo4j/{sample}/cl_contains.csv',
'results/neo4j/{sample}/cl_isin.csv',
'results/neo4j/{sample}/expr_by.csv',
'results/neo4j/{sample}/expr_ess.csv'
shell:
"""
python {input.script} \\
--path_to_cl {input.path_to_cl} \\
--path_to_umap {input.path_to_umap} \\
--path_to_mtx {input.path_to_mtx} \\
--out_dir results/neo4j/{wildcards.sample}
"""
rule hello:
output:
"woot", "hoot"
run:
for f in output:
print(f)
print(output[1])
打印“woot”,“hoot”,“hoot”。