'Initialization of the mixture parameters.
Parameters
X : array-like, shape (n_samples, n_features)
resp : array-like, shape (n_samples, n_components)'

def _initialize(self, X, resp):

    (nk, xk, sk) = _estimate_gaussian_parameters(X, resp, self.reg_covar, self.covariance_type)
    self._estimate_weights(nk)
    self._estimate_means(nk, xk)
    self._estimate_precisions(nk, xk, sk)
