Local MTT Likelihood

Computes the proposal kernel for a MultiTypeTree NodeRetype operator

Usage

local_MTT_transition_kernel(
  ED,
  st_labels,
  bit_rates,
  node_indices = NodeIndicesC(ED),
  eigen_vals = eigen(bit_rates)$values,
  eigen_vecs = eigen(bit_rates)$vectors,
  inverse_vecs = solve(eigen_vecs)
)

Arguments

ED

Structured phylogeny using internal data structure

st_labels

subset of rows from an ED object giving all nodes in a subtree

bit_rates

Backwards-in-time migration rates matrix

Value

Likelihood & log_likelihood between events