find_modes function

Find modes of a distribution

find_modes(anc, node, mode_proba_method = "integral")

Arguments

• anc: an object of class ancestralCauchy
• node: a vector of node number in the tree

Returns

a list, with opt_vals the abscissa values where the distribution is maximum, and prop_vals the approximate density under each optimal value, normalized to sum to one.