package phylogenetics

  1. Overview
  2. Docs
type vector
type matrix
type symbol
type t = matrix
val make : (symbol -> symbol -> float) -> t
val jc69 : unit -> t
val gtr : equilibrium_frequencies:vector -> transition_rates:Linear_algebra.vec -> t
val stationary_distribution : t -> vector

stationary_distribution r numerically computes the asymptotic probability distribution pi of the CTMC defined by r.

val scaled_rate_matrix : vector -> t -> t

scaled_rate_matrix pi r is a new matrix rate such that the corresponding CTMC has one expected transition per unit of time. In addition, if r is symetrical, the result has pi as stationary distribution.

val scale : t -> t

rescale matrix such that the sum of off-diagonal elements is 1.