Inter-iteration optimization of parallel EM algorithm on message-passing multicomputers

Estimation of the parameters of a probability distribution function is a complicated problem that is frequently encountered in many instances of real world problems. The Expectation Maximization (EM) algorithm often can be employed when there is a many-to-one mapping from all possible distribution patterns to the distribution governing the outcome. With its maximum likelihood (ML) formulation, optimal estimate can be made for the unknown variables after iterations until convergence. A variety of parallel methods have been proposed to boost its performance because of the complexity involved in the algorithm. Despite the efforts, the ML algorithm could not be easily adopted in practice primarily due to both intra- and inter-iteration data dependence problems resulting from the iterative nature of the algorithm. This research builds upon experimentation that demonstrated promising results in speeding up the algorithm in and between iterations using distributed memory message passing architecture.

[1]  Zhiwei Xu,et al.  Modeling communication overhead: MPI and MPL performance on the IBM SP2 , 1996, IEEE Parallel Distributed Technol. Syst. Appl..

[2]  Roy L. Streit,et al.  Maximum likelihood training of probabilistic neural networks , 1994, IEEE Trans. Neural Networks.

[3]  D. Rubin,et al.  On Jointly Estimating Parameters and Missing Data by Maximizing the Complete-Data Likelihood , 1983 .

[4]  L. Shepp,et al.  Maximum Likelihood Reconstruction for Emission Tomography , 1983, IEEE Transactions on Medical Imaging.

[5]  Rolf Hempel,et al.  The MPI Message Passing Interface Standard , 1994 .

[6]  Gordon Bell,et al.  Ultracomputers: a teraflop before its time , 1992, CACM.

[7]  Chung-Ming Chen,et al.  On Parallelizing the EM Algorithm for PET Image Reconstruction , 1994, IEEE Trans. Parallel Distributed Syst..

[8]  Roger W. Hockney,et al.  The Communication Challenge for MPP: Intel Paragon and Meiko CS-2 , 1994, Parallel Computing.

[9]  Scott T. Grafton,et al.  Results of a clinical receiver operating characteristic study comparing filtered backprojection and maximum likelihood estimator images in FDG PET studies. , 1993, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[10]  Simon R. Cherry,et al.  Fast implementations of 3D PET reconstruction using vector and parallel programming techniques , 1992 .

[11]  Álvaro Suárez Sarmiento,et al.  Overlapping computations and communications in torus networks , 1996, Proceedings of 4th Euromicro Workshop on Parallel and Distributed Processing.