Optimal Estimation of Multivariate ARMA Models