Open hyanwong opened 1 month ago
just a note that
print(
"Mean of cross-coalescence time distr: ",
2 * node_coal_weight @ ts.nodes_time / np.sum(node_coal_weight),
"Branch-mode divergence: ",
ts.divergence(sample_sets, indexes=indexes, mode='branch'),
)
will only match if the tree sequence is trimmed first (which it isn't in this example)
@nspope sent me this, which is a good core for a new tutorial on cross coalescence rates. Probably worth waiting until some of the core functionality (e.g.
rates_from_ecdf
) is implemented in tskit.