'Compute the log probability under the model.
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
logprob : array_like, shape (n_samples,)
Log probabilities of each data point in X'

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

    (logprob, _) = self.score_samples(X)
    return logprob
