inhomogenous_coal.simulate: Simulate inhomogenous coalescent,

In dhelekal/CaveDive: Phylodynamics for Clonal Expansion

Simulate inhomogenous coalescent,

Description

Simulate inhomogenous coalescent,

Usage

inhomogenous_coal.simulate(
  sampling_times,
  Neg.rate,
  Neg.rate.int,
  Neg.rate.int_inv = NULL
)

Arguments

sampling_times	times of leaves.
Neg.rate	function equal to 1/Neg(t).
Neg.rate.int	integral of 1/Neg(t) from t to t+s.
Neg.rate.int_inv	function that inverts the integral of 1/Neg(t) from t to t+s for s.

Value

A list consisting of the simulate coalescent times coalescent_times and the log-likelihood of the simulation log_likelihood.

dhelekal/CaveDive documentation built on June 11, 2024, 4:32 p.m.