multivariate GARCH-M models