ELynx.Simulate.MarkovProcess

Creation date: Thu Jan 24 09:02:25 2019.

Documentation

type ProbMatrix = Matrix R

A probability matrix, P_ij(t) = Pr (X_t = j | X_0 = i).

type State = Int

Make type signatures a little clearer.

probMatrix :: RateMatrix -> Double -> ProbMatrix

The important matrix that gives the probabilities to move from one state to another in a specific time (branch length).

jump :: PrimMonad m => State -> ProbMatrix -> Gen (PrimState m) -> m State

Move from a given state to a new one according to a transition probability matrix .

This function is the bottleneck of the simulator and takes up most of the computation time.