In this paper we develop an approximate methodology for estimating the multidimensional quantum density associated with a statistical bundle of de Broglie-Bohm trajectories. The quantum density is constructed as a discrete sum of nonequivalent Gaussian components. We incorporate the ideas of Bayesian statistical analysis and an expectation-maximization procedure to compute an approximate quantum force that drives the statistical ensemble quantum trajectories. (C) 2003 American Institute of Physics.