我正在尝试使用snakemake制作一个简单的管道,从网上下载两个文件,然后将它们合并到一个输出中。
我认为可行的是以下代码:
dwn_lnks = {
'1': 'https://molb7621.github.io/workshop/_downloads/sample.fa',
'2': 'https://molb7621.github.io/workshop/_downloads/sample.fa'
}
import os
# association between chromosomes and their links
def chromo2link(wildcards):
return dwn_lnks[wildcards.chromo]
rule all:
input:
os.path.join('genome_dir', 'human_en37_sm.fa')
rule download:
output:
expand(os.path.join('chr_dir', '{chromo}')),
params:
link=chromo2link,
shell:
"wget {params.link} -O {output}"
rule merger:
input:
expand(os.path.join('chr_dir', "{chromo}"), chromo=dwn_lnks.keys())
output:
os.path.join('genome_dir', 'human_en37_sm.fa')
run:
txt = open({output}, 'a+')
with open (os.path.join('chr_dir', "{chromo}") as file:
line = file.readline()
while line:
txt.write(line)
line = file.readline()
txt.close()
此代码返回错误:
No values given for wildcard 'chromo'. in line 20
另外,在合并规则中,运行中的python代码不起作用。
snakemake 包中的教程没有涵盖足够的示例来了解非计算机科学家的详细信息。如果有人知道一个很好的资源来学习如何使用snakemake,如果他们能分享我将不胜感激:)。