Snakemake 中的 MissingRule 异常答案

【问题标题】：MissingRule Exception in SnakemakeSnakemake 中的 MissingRule 异常
【发布时间】：2021-09-27 12:04:39
【问题描述】：

我是snakemake 工作流管理的新手，我正在努力掌握通配符输入的工作原理。我尝试对一些 SRR 数据进行 QC，但 snakemake 给出“MissingRuleException 错误”。

我的配置文件(config.yaml)包含内容：

样本：sample.csv

路径：/Users/path/Bioinformatics/srr_practice

sample.csv 是

sample_name,fq1

A,SRR11412215

B,SRR11412216

C,SRR11412217

D,SRR11412218

E,SRR11412219

蛇文件

import os
import pandas as pd

configfile:"config.yaml"

samples=pd.read_csv(config["samples"], sep=",").set_index("sample_name", drop=False)


def get_fastq(wildcards):

        units=samples.loc[wildcards.sample]
fq=units["fq1"]
        return expand(os.path.join(config["path"], "{fq}.fastq.gz"), fq=fq)

rule all:

    input:

        expand(os.path.join(config["path"], "fastq_output/{sample}.fastqc.html"),sample=samples["fq1"].to_list()),
        expand(os.path.join(config["path"], "fastq_output/{sample}_fastqc.zip"), sample=samples["fq1"].to_list())

rule fastq:

    input:

            get_fastq,
    output:

            zip=os.path.join(config["path"], "fastq_output/{wildcards.sample_name}_fastqc.zip"),
                html=os.path.join(config["path"], "fastq_output/{wildcards.sample_name}_fastqc.html")

    wrapper:
            "0.78.0/bio/fastqc"

snakemake -np 蛇文件

错误构建工作的 DAG... 缺少规则异常：没有生成 Snakefile 的规则（如果您使用输入函数，请确保它们不会引发意外异常）。

【问题讨论】：

标签： python bioinformatics snakemake

【解决方案1】：

通过发出命令snakemake -np Snakefile，您要求snakemake 生成Snakefile，但它不知道该怎么做。

如果您的文件名为Snakefile，则无需指定名称，如果名称不同，则可以使用-s 选项指定。所以现在，运行snakemake -n 应该足以告诉你snakemake 会运行什么。

【讨论】：