loglik_HR function

Compute Huesler-Reiss log-likelihood, AIC, and BIC