I used this function when I was working on the paper about online incremental learning two years ago. But cannot recall it when someone asked me on this recently. Put it here to memorize it:
where big theta is the model (represented by a set of parameters). The total number of parameter of the GMM is N(k), and k is the number of Gaussian components. y(obs) are the observations. n is the number of observations.