Estimation in Mixtures of Two Normal Distributions

Abstract

This paper is concerned primarily with the method of moments in dissecting a mixture of two normal distributions. In the general case, with two means, two standard deviations, and a proportionality factor to be estimated, the first five sample moments are required, and it becomes necessary to find a particular solution of a ninth degree polynomial equation that was originally derived by Karl Pearson [10]. A procedure which circumvents solution of the nonic equation and thereby considerably reduces the total computational effort otherwise required, is presented. Estimates obtained in the simpler special case in which the two standard deviations are assumed to be equal, are employed as first approximations in an iterative method for simultaneously solving the basic system of moment equations applicable in the more general case in which the two standard deviations are unequal. Conditional maximum likelihood and conditional minimum chi-square estimation subject to having the first four sample moments equated to corresponding population moments, are also considered. An illustrative example is included.