Example Scheme code that implements the demographic scenario of Figure 4. The Scheme code specifies both the population structures and the migration rates between populations, and simulates a sample of 100 individuals from population P1 and 50 from each of populations P3 and P4, with the merging of populations P3 and P4 into p2 at time 1.5 and the merging of populations P1 and P2 at time 3. The merge of P1 and P2 is followed by a bottleneck followed by a period of constant population size of f*2N = 10N.
Mailund et al. BMC Bioinformatics 2005 6:252 doi:10.1186/1471-2105-6-252