For the multivariate equation-error moving average system, a multivariate maximum likelihood multi-innovation extended stochastic gradient (M-ML-MIESG) algorithm is delivered. The key is to decompose the system into several regressive identification subsystems according to the number of the system outputs. Then, a multivariate maximum likelihood extended stochastic gradient algorithm is presented to estimate the parameters of these subsystems. The M-ML-MIESG algorithm has higher parameter estimation accuracy than the multivariate extended stochastic gradient algorithm. The simulation examples indicate that the proposed methods work well.