Default method to compute the hessian of the log-likelihood.

In case a hess_loglik method is not provided, this function will be used. It computes the hessian of the log-likelihood function using numerical differentiation.

# S3 method for class 'likelihood_model'
hess_loglik(model, control = list(), ...)

Arguments

model

The likelihood model

control

Custom arguments to pass to numDeriv::hessian.

...

Additional arguments (to pass into loglik)

Value

A function(df, par, ...) that computes the Hessian matrix of the log-likelihood evaluated at par given data df.