Estimation of vector ARMAX models