Skip to contents

Performs a coalescent node split move (Ewing et al. 2004). Splits a migration node immediately above a coalescent node into two migration nodes placed immediately below the coalescent node

Usage

coalescent_split(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

node_indices

Vector giving row indices for node labels

Value

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