tweaking to use mvtnorm more