Returns the Hessian matrix of the profile log-likelihood.
Notes
The Hessian matrix that is returned is computed with respect
to the parameterization defined by this model’s use_sqrt
attribute. The input value params can be with respect to
any parameterization.