Please note: I posted this first on Mathoverflow. Someone there advised me that on stats.stackexchange the question might fit better here. This is the link to the original post. I currently have to ...
I know that Mclust does the fit on its own but I am trying to implement an optimization with the aim to generate a mixture of 2 gaussians with the combine moments as closed as possible to the moment ...