Post date: Jun 29, 2015 10:43:24 PM
In addition to my original simulations (see draft in previous post) I am now analyzing weak constant and (linear) fluctuating selection (s = 0.01 rather than 0.1).
## weak constant selection (infiles made with mkRSimsWeakSel.pl)
cd /local/scratch/
sleep 10
/home/A01963476/bin/spatpg -g /labs/evolution/projects/spatgen/rsims/mod_rsimWkSelG9.txt -e /labs/evolution/projects/spatgen/rsims/rsimWkSelE9.txt -n 25000 -b 10000 -t 15 -o out_simWkSelG9rep2.hdf5
scp out_simWkSelG9rep2.hdf5 /labs/evolution/projects/spatgen/mcmc/
## weak linear fluctuating selection (infiles made with mkRSimsWeakFluctL.pl)
cd /local/scratch/
sleep 0
/home/A01963476/bin/spatpg -g /labs/evolution/projects/spatgen/rsims/mod_rsimWeakFluctLG8.txt -e /labs/evolution/projects/spatgen/rsims/rsimWeakFluctLE8.txt -n 25000 -b 10000 -t 15 -o out_simWeakFluctLG8rep0.hdf5
scp out_simWeakFluctLG8rep0.hdf5 /labs/evolution/projects/spatgen/mcmc/
I am also adding an analysis with 2N = 50 and 100 (25 or 50 individuals sampled); see mkSub.pl script.
cd /local/scratch/
sleep 5
/home/A01963476/bin/spatpg -g /labs/evolution/projects/spatgen/rsims/sub50mod_rsimFluctLG7.txt -e /labs/evolution/projects/spatgen/rsims/rsimFluctLE7.txt -n 25000 -b 10000 -t 15 -o out_simFluctLG7sub50rep1.hdf5
scp out_simFluctLG7sub50rep1.hdf5 /labs/evolution/projects/spatgen/mcmc/