Is there any way to increase the accuracy of "lncdfmvn" when you are dealing with a 4-dimensional normal distribution?
lncdfmvn uses a global variable
_lncdfn_eps to control convergence. The default value is 0.0001. Setting it to something smaller might help get you a little more accuracy.