Expansion/migration script and random seeding questions

[Cameron Grey]

Hi Ryan,

I've been using dadi to model the initial expansion and gene flow of a recent invasion, and I've started on three populations. I have SNP data for the 3 populations and generation times for how long they have (probably) been there. I wanted a model that could examine likely initial expansion events as well as migration post founding of each population. I have scripts for each scenario, and have been comparing AIC values calculated from output LL. 

Specifically -  I tested different expansion scenarios first with the nuPOP parameters, and found the best fit expansion scenario. From there, I have added migration overtop of the expansion with the m_POP parameters in my scripts.

1. Question about migration -  I have included a script where I have both the expansion scenario and one of the migration scenarios I am testing overtop the expansion results. I was wondering if this was a reasonable/informative approach to expansion and migration. Because the data is from one time point, how does dadi differentiate initial expansion events from migration that occurs later on given that these shared allele frequency changes can be quite subtle for recent events? I'm curious because I just got a couple results suggesting migration in the same direction as my likely expansion scenario (I tested just expansion first).

2. Question about seeding - I realized that I was missing a random seeding for the different replicates of the model, and I am thinking of adding in np.random.seed() to the run_opt() portion. Any advice on this approach? I appreciate it.

Thank you for all of your help and for making such an interesting model!

-Cameron