how to compute likelihood for phylogenetic trees