38 lines
1.1 KiB
Bash
Executable File
38 lines
1.1 KiB
Bash
Executable File
#!/usr/bin/env bash
|
|
# H2 genome-wide: per-chromosome HCT-116 perturbation analysis (chr1..22).
|
|
# Idempotent. Disables per-chrom figures by default.
|
|
#
|
|
# Usage:
|
|
# bash experiments/h2_rewiring/run_genome.sh
|
|
# bash experiments/h2_rewiring/run_genome.sh chr2 chr5
|
|
# SKIP_FIGURES= bash experiments/h2_rewiring/run_genome.sh # render figures
|
|
set -euo pipefail
|
|
|
|
REPO="$(cd "$(dirname "${BASH_SOURCE[0]}")/../.." && pwd)"
|
|
|
|
export SKIP_FIGURES="${SKIP_FIGURES-1}"
|
|
|
|
if [ $# -gt 0 ]; then
|
|
CHROMS=("$@")
|
|
else
|
|
CHROMS=()
|
|
for i in $(seq 1 22); do CHROMS+=("chr$i"); done
|
|
fi
|
|
|
|
START=$(date +%s)
|
|
N=${#CHROMS[@]}
|
|
i=0
|
|
for CHROM in "${CHROMS[@]}"; do
|
|
i=$((i + 1))
|
|
echo ""
|
|
echo "================================================================"
|
|
echo "[H2] $CHROM ($i / $N)"
|
|
echo "================================================================"
|
|
CHROM="$CHROM" bash "$REPO/experiments/h2_rewiring/run.sh"
|
|
done
|
|
|
|
ELAPSED=$(( $(date +%s) - START ))
|
|
echo ""
|
|
echo "=== H2 genome-wide complete in ${ELAPSED}s ==="
|
|
echo "Next: python experiments/h2_rewiring/aggregate.py"
|