compute_betas_from_shifts
computes the optimal values at the nodes and tips of the
tree, given the value at the root and the list of shifts occuring in the tree.
It assumes an OU model.
compute_betas_from_shifts(phylo, optimal.value, shifts)
phylo
.