Estimating the normal-inverse-Wishart distribution