bgc implements Bayesian estimation of genomic clines to quantify introgression at many loci.
| Version | module |
|---|
| 1.03 | bgc/1.03 |
#!/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
For more information visit https://sites.google.com/site/bgcsoftware/
