# 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,