Swap dependency on MASS with mvtnorm
Test the 2D case