Depiction of the multiple population-pair divergence model used for the ABC estimates of Ψ, E(τ), and Ω. (A): The white lines depict a gene tree with TMRCA being the time to the gene sample's most recent common ancestor, and the black tree containing the gene tree is the population/species tree. (B): Parameters in the multiple population-pair divergence model. The population mutation parameter, θ, is 2Nμ where 2N is the summed haploid effective female population size of each pair of daughter populations (μ is the per gene per generation mutation rate). The time since isolation of each population pair is denoted by τ (in units of 2NAVE generations, where NAVE is the parametric expectation of N across Y population pairs given the prior distribution). Population mutation parameters for daughter populations a and b are θa and θb, whereas θ 'a and θ'b are the population mutation parameters for the sizes of daughter populations a and b at the time of divergence until τ' (length of bottleneck). The daughter populations θ 'a and θ'b then grow exponentially to sizes θa and θb. The population mutation parameter for each ancestral population is depicted as θA. The migration rate between each pair of daughter populations is depicted as M (number of effective migrants per generation). (C): Example of four population-pairs where parameters in (B) are drawn from uniform priors.
Hickerson et al. BMC Bioinformatics 2007 8:268 doi:10.1186/1471-2105-8-268