Returns a function that computes the Hessian matrix of the log-likelihood. Uses user-provided Hessian function if available, otherwise falls back to numerical differentiation via numDeriv::hessian.
Usage
# S3 method for class 'dfr_dist'
hess_loglik(model, ...)