'Compute the per-sample average log-likelihood of the given data X.
Parameters
X : array-like, shape (n_samples, n_dimensions)
List of n_features-dimensional data points. Each row
corresponds to a single data point.
Returns
log_likelihood : float
Log likelihood of the Gaussian mixture given X.'

def score(self, X, y=None):

    return self.score_samples(X).mean()
