Local DTA

Updates the migration history on a subtree under the DTA model

Usage

local_DTA_subtree_proposal(
  ED,
  st_labels,
  fit_rates,
  n_deme = nrow(fit_rates),
  node_indices = NodeIndicesC(ED),
  eigen_decomp = eigen(fit_rates),
  inverse_vecs = solve(eigen_decomp$vectors)
)

Arguments

ED

ED representation of entire structured phylogenetic tree including virtual migrations needed to isolate subtree

st_labels

ED representation of subtree only

fit_rates

Forwards-in-time migration rates matrix

node_indices

Vector of node indices. Element i gives the row of label i

eigen_decomp

Eigendecomposition of fit_rates (output from eigen(fit_rates))

inverse_vecs

Inverse of matrix of eigenvectors (output from solve(eigen(fit_rates)$vectors))

Value

List of proposed DTA migration histories with the same leaves as input ED