mushi.utils.M

M(n, t, y)[source]

The matrix \(\mathbf M\) defined in the paper’s appendix.

Parameters:
  • n (int) – the number of sampled haplotypes \(n\)

  • t (Array) – time grid, starting at zero and ending at np.inf

  • y (Array) – population size in each epoch

Return type:

Array

Returns:

\((n-1)\times m\) matrix, where \(m\) is the number of epochs (the length of the y argument)