Negative-binomial likelihood function
likelihoodNegbin(tab, mu, sigma)
tab
Table of nodes
mu
Clock rate parameter
sigma
Std of per branch clock rate
log-likelihood