bgc

#!/bin/bash

#SBATCH -p eng
#SBATCH -n 1
#SBATCH -t 24:00:00
#SBATCH --qos=nogpu
#SBATCH -J bgctest

cd ${SLURM_SUBMIT_DIR}

module load bgc/1.03

bgc -a p0in.txt -b p1in.txt -h admixedin.txt -M map.txt -O 0 -x 50000 -n 25000 -p 1 -q 1 -N 1 -m 1 -D 0.5 -t 5 -E 0.0001 -d 1 -s 1 -I 0 -u 0.04

estpost -i mcmcout.hdf5 -p LnL -o ln1 -s 2 -w 0
estpost -i mcmcout.hdf5 -p alpha -o a1 -s 2 -w 0

estpost -i mcmcout.hdf5 -p alpha -o a.out -s 0 -c 0.95 -w 0
estpost -i mcmcout.hdf5 -p beta -o b.out -s 0 -c 0.95 -w 0

estpost -i mcmcout.hdf5 -p gamma-quantile -o qa.out -s 0 -c 0.95 -w 0
estpost -i mcmcout.hdf5 -p zeta-quantile -o qb.out -s 0 -c 0.95 -w 0

exit