I have a practical question. Given a phylogenetic tree (i.e. topology plus branch lengths) I’d like to calculate a likelihood of observing that tree shape under a particular generating model, and also find ML estimates of the model parameters.
I don’t know of any package that can do this, but it doesn’t seem impossible. Does anyone have any pointers, or know of existing methods?