# 3 values
no_mig_params = [0.357, 0.7876, 0.1456]
#"sym_mig"
#4 Values
sym_mig_params = [0.6832, 1.4416, 0.6072, 2.5373]
#"asym_mig"
#5 Values
asym_mig_params = [0.938, 3.558, 0.5284, 0.1714, 5.8838]
#======================================================================================
# Call the function with the relevant arguments.
# Divergence with no migration
Optimize_Functions.Optimize_Round2(pts, fs, outfile, reps, maxiter, "no_mig",
no_mig_params)
# Split into two populations, with continuous symmetric migration.
Optimize_Functions.Optimize_Round2(pts, fs, outfile, reps, maxiter, "sym_mig",
sym_mig_params)
# Split into two populations, with continuous asymmetric migration.
Optimize_Functions.Optimize_Round2(pts, fs, outfile, reps, maxiter, "asym_mig",
asym_mig_params)
#===========================================================================
#clock the amount of time to complete the script
t_finish = datetime.now()
elapsed = t_finish - t_begin
print '\n', '\n', "-----------------------------------------------------------------------------------------------------"
print "Finished all analyses!"
#======================================================================================
# Now call the function with the relevant arguments.
# There are many models to test here. A brief definition is given for each, but the actual
# models are defined in the Models_2D.py script. The first 15 were implemented in Portik
# et al. 2016 (doi: 10.1111/mec.14266), the following 9 are newer for various projects.
# Here it is set up to call each model one by one sequentially, which could finish relatively quickly.
# If it takes too long, create multiple verisions of this script, block out some models (use hashes or delete),
# and execute one version for every core you have available. It will greatly speed up these steps,
# and sometimes if extrapolations fail the script will crash too and this could prevent it from
# happening too many times.
# Standard neutral model, populations never diverge
Optimize_Functions.Optimize_Round2(pts, fs, outfile, reps, maxiter,
"no_divergence", no_divergence_params)
# Split into two populations, no migration.
Optimize_Functions.Optimize_Round2(pts, fs, outfile, reps, maxiter, "no_mig",
no_mig_params)
# Split into two populations, with continuous symmetric migration.
Optimize_Functions.Optimize_Round2(pts, fs, outfile, reps, maxiter, "sym_mig",
sym_mig_params)
# Split into two populations, with continuous asymmetric migration.
Optimize_Functions.Optimize_Round2(pts, fs, outfile, reps, maxiter, "asym_mig",
asym_mig_params)
# Split with continuous symmetric migration, followed by isolation.
Optimize_Functions.Optimize_Round2(pts, fs, outfile, reps, maxiter,