How to find the best fit of a GMM trying multiple components in Matlab

In R, there is Mclust, which would by default run the EM algorithm considering G=9 components. In Matlab, it seems that it only performs the EM algorithm for a single given number of components. Is there either: A function that is like

Continue reading →