'Compute the average log-likelihood of the samples
Parameters
X : array, shape (n_samples, n_features)
The data
Returns
ll : float
Average log-likelihood of the samples under the current model'

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

    return np.mean(self.score_samples(X))
