Snakemake中间规则未执行：依赖未识别问题求助

核心问题原因分析 #

原脚本无法触发melt_bed_null和intersect_scores_null规则，本质是依赖链断裂，具体原因包括：

• checkpoint后续依赖推导错误：fetch_null_scores函数试图从checkpoint输出目录匹配phyloP_scores{i}.txt，但该文件是intersect_scores_null的输出，并非checkpoint的直接产物。Snakemake初始构建DAG时，无法识别这些文件的生成逻辑，直接跳过中间规则。
• 路径与文件匹配歧义：所有中间文件都存放在null/{gene_symbol}目录下，不同规则的输出文件混在一起，加上glob_wildcards匹配目标错误，导致Snakemake无法建立正确的依赖关联。
• 配置与变量冗余：手动加载yaml配置的同时使用configfile指令，造成冗余；gene_symbol提前硬编码，导致规则中的wildcard无法正常传递参数。

修正脚本的关键改进点 #

修正后的脚本通过以下调整修复了依赖链：

• 修正checkpoint wildcard提取逻辑：fetch_null_scores改为从checkpoint直接输出的gene{i}.bed中提取wildcard，再通过expand生成最终需要的null_scores/{gene_symbol}/phyloP_scores{i}.tsv，让Snakemake能反向推导触发中间规则。
• 拆分中间文件目录：将gene{i}.bed、melted_gene{i}.bed、phyloP_scores{i}.tsv分别放在null_beds、null_beds_melted、null_scores独立目录下，避免路径冲突，依赖关系更清晰。
• 简化配置管理：移除手动yaml加载，直接使用Snakemake内置的config变量，避免变量硬编码，保证wildcard正常工作。

技术建议 #

• 使用checkpoint时，必须从checkpoint的直接输出文件中提取wildcard，再通过expand推导后续依赖，不能直接匹配后续规则的输出文件，否则会断裂依赖链。
• 中间文件按规则类型拆分独立目录，避免同目录下多类型文件混淆，帮助Snakemake清晰解析依赖关系。
• 优先使用Snakemake内置的configfile机制加载配置，不要手动重复加载yaml文件，减少变量冲突。
• 调试时可在依赖函数中加入print语句输出中间结果，或执行snakemake --dag | dot -Tpdf > dag.pdf生成依赖图，直观检查DAG结构是否正确。

内容的提问来源于stack exchange，提问作者Rebecca Eliscu