logLik_hapassoc function

Return log-likelihood