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 (ndarray) – time grid, starting at zero and ending at np.inf

  • y (ndarray) – population size in each epoch

Return type

ndarray

Returns

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