Due to the non-existence of a closed-form expression for sub-Gaussian α-stable densities, the M-step of the Expectation-Maximization (EM) algorithm for the SubGaussian α-Stable Mixture Models (SGαSMMs) is intractable, and the EM algorithm for SGαSMM is still an open problem. SGαSMM can model non-homogeneous Gaussian data, accommodate outliers, and high leverage data points, which are concepts of primary importance in robust mixture models. These models are robust and useful tools in modeling heterogeneous data with outlier observations, such as clustering fnancial or impulsive data. In this paper, a new EM algorithm based on a combination of EM (the frst part) and stochastic EM (the second part) algorithms, is used to obtain the maximum likelihood estimators of the parameters of SGαSMM in the M-step. In the frst part, the model parameters, except αs, are estimated from an analytical form via EM. In the second part, based on a stochastic EM, the maximum likelihood estimator of α, in each component, is calculated from pseudo-simulated data obtained by suitable rejection sampling. The efciency of the proposed algorithm is illustrated by using both real and simulated data.
