THE EMMIX ALGORITHM FOR THE FITTING OF MIXTURES OF NORMAL AND t-COMPONENTS

We consider the tting of normal mixture models to multivariate data, using maximum likelihood via the EM algorithm. This approach requires the initial speciication of an initial estimate of the vector of unknown parameters, or equivalently, of an initial clas-siication of the data with respect to the components of the mixture model under t. We describe an algorithm called EMMIX that automatically undertakes this tting, including the provision of suitable initial values if not supplied by the user. The EMMIX algorithm has several options, including the option to carry out a resampling-based test for the number of components in the mixture model.