Post date: Dec 16, 2014 12:4:52 AM
DESCRIBE PROGRAM HERE
Genotype infiles in projects/popanc_sims/sims/.
10, 100 or 1000 generations, currently all with F = 0.3, no selection, ancestral allele frequencies beta(0.8,0.8), no drift in terms of allele frequencies (might want to change this).
Get these by first running simGenotypes.pl on the gen_demong_gens*.txt file. Then reHeader.pl to add the number of individual and loci to each file.
I used the wrap_qsub* script to try short cross-validation runs on a few reps. just to make sure everything works (see below).
cd /local/scratch/
sleep 5
popanc -o out_r2f0.3gen_demog_gens10.hdf5 -m 10000 -b 10000 -t 5 -f 0 -w 0 -v 1 -c 0.1,1,10 /labs/evolution/projects/popanc_sims/sims/genoP0_r2f0.3gen_demog_gens10.txt /labs/evolution/projects/popanc_sims/sims/genoP1_r2f0.3gen_demog_gens10.txt /labs/evolution/projects/popanc_sims/sims/genoAdmx_r2f0.3gen_demog_gens10.txt > /labs/evolution/projects/popanc_sims/mcmc/log_r2f0.3gen_demog_gens10
scp out_r2f0.3gen_demog_gens10.hdf5 /labs/evolution/projects/popanc_sims/mcmc/