ENH: add untested likelihood calculation functions