If only mode_compute=6 can be used to find mode, since 9 or 4 has ‘‘hessian matrix at the “mode” is not positive definite’’ issue. Then after I start at the mode found by 6 and run MCMC, say , 1 million draws, and finally all
estimated variables converge. Then would this result reliable? Should I do other things to enhance the probability of finding a global mode?