coaltimeSim: Simulate from inhomogeneous, heterochronous coalescent.
In jrfaulkner/bnps: Shrinkage-Prior Markov Random Field Smoothing

Description Usage Arguments Value Examples

For simulating coalescent times using the method of thinning. This function is based on code written by Michael Karcher for the phylodyn package.

1	coaltimeSim(samp_times, n_sampled, traj, lower_bound, ...)

`samp_times`	numeric vector of sampling times.
`n_sampled`	numeric vector of samples taken per sampling time.
`traj`	function that returns effective population size at time t.
`lower_bound`	numeric lower limit on `traj` function on its support.
`...`	additional arguments to be passed to `traj` function.

A list containing vectors of coalescent times coal_times, intercoalescent times intercoal_times, and number of active lineages lineages, as well as passing along samp_times and n_sampled.