Maximum likelihood estimation for vector autoregressions with multivariate stochastic volatility