The heavy-tailed Multivariate Normal Inverse Gaussian (MNIG) distribution is a recent variance-mean mixture of a multivariate Gaussian with a univariate inverse Gaussian distribution. Due to the complexity of the likelihood function, parameter estimation by direct maximization is exceedingly difficult. To overcome this problem, we propose a fast and accurate multivariate ExpectationMaximization (EM) ...