Hello everyone,

Could anyone help me to solve the following problem please?

If the inverse of Hessian matrix (W) , which is estimated directly from the optimization package of Chris Sim, is not positive-definite, how can we calculate its square root to make a jump distribution N(mu, W) for random walk MCMC as suggested in the literature please?

Thank you very much,

Best regards,

A