This paper addresses the problem of estimating the parameters of multivariate ARMA processes by using higher-order statistics called cumulants. The main objective in this paper is to extend the idea of the q-slice algorithm in univariate ARMA processes to multivariate ARMA processes. It is shown for a multivariate ARMA process that the MA coefficient matrices can be estimated up to postmultiplication of a permutation matrix by using the third-order cumulants and of an extended permutation matrix by using the fourth-order cumulants. Simulation examples are included to demonstrate the effectiveness of the proposed method.