The function normalmixEM in R package mixtools is employed for fitting multi-component normal mixture models.
EMFit(x, k0, mean_constr, sd_constr, npara, d0)
Normal mixture model fit and BIC value of the log-likelihood
A numeric vector
Number of components in the normal mixture model
A constrain on means of components
A constrain on standard deviations of components
Number of parameters
Number of times for fitting mixture model using different starting values