'Compute the weighted log probabilities for each sample.
Parameters
X : array-like, shape (n_samples, n_features)
List of n_features-dimensional data points. Each row
corresponds to a single data point.
Returns
log_prob : array, shape (n_samples,)
Log probabilities of each data point in X.'

def score_samples(self, X):

    self._check_is_fitted()
    X = _check_X(X, None, self.means_.shape[1])
    return logsumexp(self._estimate_weighted_log_prob(X), axis=1)
