test.sh

#!/bin/bash

set -euxo pipefail

SEED1=123
SEED2=456
BASE1="test/out_${SEED1}"
BASE2="simulation/${BASE1}_stats/rep${SEED2}"
BASEX="${BASE2}_xval_"

export PYTHONPATH="$PWD:$PWD/simulation${PYTHONPATH:+:${PYTHONPATH}}"
echo "PYTHONPATH=$PYTHONPATH"
pushd simulation
rm -rf test/out_*
slim -d 'OUTDIR="test"' -s $SEED1 fecundity_regulation.slim
python compute_stats.py ${BASE1}.trees 20 $SEED2
python plot_maps.py ${BASE1}.spstats.csv
popd
python setup_crossval_grid.py simulation/test/crossvalidation.json 3 0.1 1 3 0.1 1.0 ${BASE2}
for j in $BASEX*/xval_params.json;
do
    python crossvalidation.py --json $j --max_iter 5
done

python crossvalidation.py -H -e 1 -i 1 --json `ls $BASEX*/xval_params.json | head -n1` --max_iter 5

for p in $(ls $BASEX*/*.pkl | head -n 2)
do
    python plot_crossvalidation.py $p
    q=${p%pkl}truth.qmd
    cp simulation/compute_truth.qmd $q
    quarto render $q -P results_file:$(basename $p)
done
QMD=${BASE2}.summary.qmd
cp simulation/summarise_crossvalidation.qmd $QMD
quarto render $QMD -P basename:$(basename $BASE2)

echo "All done! Produced:"
find simulation/${BASE1}_stats -not -path "${QMD%.qmd}_files/*"