In financial time series analysis we encounter several instances of non-negative valued processes (volumes, trades, durations, realized volatility, daily range, and so on) which exhibit clustering and can be modeled as the product of a vector of conditionally autoregressive scale factors and a multivariate iid innovation process (vector Multiplicative Error Model).Two novel points are introduced in this paper relative to previous suggestions: a more general specification which sets this vector MEM apart from an equation by equation specification; and the adoption of a GMM-based approach which bypasses the complicated issue of specifying a general multivariate non-negative valued innovation process.A vMEM for volumes, number of trades and realized volatility reveals empirical support for a dynamically interdependent pattern of relationships among the variables on a number of NYSE stocks. * This paper develops some ideas introduced in Cipollini, Engle and Gallo (2006) where estimation was based in the framework of estimating functions. Without implicating, we acknowledge comments by Nour Meddahi and Kevin Sheppard which led us to present the Estimating Functions approach in a more familiar GMM notation. We acknowledge financial support from the Italian MIUR under grant PRIN 2006131140 004.