The computations of likelihood or posterior distribution of parameters of complex population genetics models are common tasks in computational biology. The numerical results of these approaches are often found by Monte Carlo simulations. Much of the recent work of Monte Carlo approaches to...

This paper is concerned with hierarchical clustering of long binary sequence data. We propose two alternative improvements of the EM algorithm used in Chen and Lindsay (2006). One is the FixEM. It is just the regular EM but we no longer update the weights πs used in the ancestral mixture...