Snakemake如何在dry run时仅显示Job counts信息？

我完全理解你的痛点——大流水线用-n做dry run时，Job counts藏在几百行任务列表的末尾，中途中断真实运行又太麻烦。给你几个实用的方案，用Snakemake自带参数就能直接拿到干净的Job counts输出：

方案1：直接输出简洁的Job counts #

用--quiet参数抑制dry run的任务列表输出，只保留统计信息：

snakemake -n --quiet

运行后会直接得到类似这样的输出：

Job counts:
	count jobs	504
	BMO	1
	all	504
	fit_nbinoms	517
	motifs_in_peaks	503
	motifs_outside_peaks	2029
Nothing to be done.

如果连最后的状态提示（比如Nothing to be done）都不想看到，再加个过滤：

snakemake -n --quiet | grep -v "Nothing to be done"

就能得到你想要的纯统计内容。

方案2：结构化提取（适合脚本自动化） #

如果需要把Job counts数据用于后续脚本处理，可以用--stats生成JSON格式的统计文件，再用工具解析：

# 用jq工具快速提取（需提前安装jq）
snakemake -n --stats job_stats.json && jq '.job_counts' job_stats.json

# 没有jq的话，用Python解析输出成易读格式
snakemake -n --stats job_stats.json && python -c "
import json
with open('job_stats.json') as f:
    stats = json.load(f)
print('Job counts:')
for rule, count in stats['job_counts'].items():
    print(f'{rule}\t{count}')
"

这个方法能拿到结构化的统计数据，方便后续自动化分析。

补充说明 #

这些方法本质上还是会完成dry run的依赖解析（毕竟要计算准确的Job counts必须先遍历所有规则和依赖），但--quiet避免了输出大量冗余的任务列表，让你不用在几百行输出里翻找统计结果，效率提升很多。

内容的提问来源于stack exchange，提问作者Ricardo A.