Skip to content

Commit cc9272a

Browse files
ahmed-shuaibiclaude
andcommitted
fix: DIG stage writes bmr_pmfs.dig.csv directly to the cohort dir
The DIG provider names its output bmr_pmfs.dig.csv (not bmr_pmfs.csv) and needs the out dir to exist; point it at output/pancan/<C> directly (no temp dir / cp), which does not clobber CBaSE's bmr_pmfs.csv. Validated on ACC: CBaSE + MutSig2CV (Docker) both run cleanly on the cBioPortal pancan MAF; DIG produces a 21MB per-gene PMF. Co-Authored-By: Claude Opus 4.8 (1M context) <noreply@anthropic.com>
1 parent fda6717 commit cc9272a

1 file changed

Lines changed: 3 additions & 6 deletions

File tree

scripts/run_cohort_pipeline.sh

Lines changed: 3 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -32,14 +32,11 @@ else log "skip cbase"; fi
3232
N=0
3333
[ -f "${ROOT}/${C}/count_matrix.csv" ] && N=$(( $(wc -l < "${ROOT}/${C}/count_matrix.csv") - 1 ))
3434

35-
# 2. DIG --------------------------------------------------------------------
35+
# 2. DIG (writes bmr_pmfs.dig.csv directly; does not clobber CBaSE's bmr_pmfs.csv) ---
3636
if [ ! -f "${ROOT}/${C}/bmr_pmfs.dig.csv" ] && [ "$N" -gt 0 ]; then
3737
log "DIG generate (N=${N})"
38-
rm -rf "/tmp/dig_${C}"
39-
if "$DIALECT" generate -m "$MAF" -o "/tmp/dig_${C}" --bmr dig \
40-
--dig-results "$DIG_RESULTS" --dig-samples "$N" -r hg19; then
41-
cp "/tmp/dig_${C}/bmr_pmfs.csv" "${ROOT}/${C}/bmr_pmfs.dig.csv"
42-
else log "STAGE-FAIL dig"; fi
38+
"$DIALECT" generate -m "$MAF" -o "${ROOT}/${C}" --bmr dig \
39+
--dig-results "$DIG_RESULTS" --dig-samples "$N" -r hg19 || log "STAGE-FAIL dig"
4340
else log "skip dig"; fi
4441

4542
# 3. MutSig2CV (Docker) -----------------------------------------------------

0 commit comments

Comments
 (0)