Extracts standard errors from the Hessian of a log-likelihood.
SE(theta_i) = sqrt(diag(solve(I(theta)))) where I = -H.
Arguments
- hess
The Hessian matrix (from
hessian())- is_loglik
If TRUE, treats hess as Hessian of log-likelihood and uses -H⁻¹. If FALSE, uses H⁻¹ directly.