Machine Learning, Clustering, Expectation Maximization, Gaussian Mixture Model, Image Clustering.
multivariate_normal.pdf(data[i],means[k],covariances[k],allow_singular=True) # Add up responsibilities over each data point and normalize row_sums = resp.sum(axis=1)[:, np.newaxis] resp = resp / row_sums return resp Remember this responsibility is our latent variable. M-Step: Now cluster