Several of you reported that with a small number of populations the coalescence time of the 1-d stepping stone model was less than that of the island model. Since your simulations showed that pattern, I gave you full credit even though I thought the result might be a statistical artifact rather than real. But since it was several of you who saw the phenomenon, I decided to check further by significantly increasing the number of replications, which makes it much less likely that any differences would be a statisctical artifact. I also modified run_simulation() to report both the mean and the median time to coalescence to see if it made a difference which measure of central tendency we used.1


NOTE: I store the result of run_simulation() in result not because I’m going to use it, but to suppress output of all of the individual simulation results. The table that is printed after each simulation provides the mean and median coalescence time for each model.

I set n_reps to 10,000 in the code that’s hidden. That’s 10 times greater than the default.

n_pops = 5

result <- run_simulation(n_pops = 5, m = 0.01, n_reps = n_reps)

result <- run_simulation(n_pops = 5, m = 0.1, n_reps = n_reps)