Migration Pair Birth MCMC Move

Performs a Migration pair birth move (Ewing et al. 2004). Adds two migration nodes on an edge selected uniformly at random from a structured coalescent process, allocating a deme for the added edge such that a migration event does not target its origin deme.

Usage

migration_pair_birth(ED, n_deme, node_indices)

Arguments

ED

Extended data object; matrix with columns Node ID, parent, child 1, child 2, deme, node age

n_deme

Number of distinct demes in the population

Value

Updated extended data object with the proposal from the migration pair birth move