'Estimate the parameters of the Gaussian distribution.
Parameters
nk : array-like, shape (n_components,)
xk : array-like, shape (n_components, n_features)'

def _estimate_means(self, nk, xk):

    self.mean_precision_ = (self.mean_precision_prior_ + nk)
    self.means_ = (((self.mean_precision_prior_ * self.mean_prior_) + (nk[:, np.newaxis] * xk)) / self.mean_precision_[:, np.newaxis])
