crabs.loglikelihood function

Compute likelihood