'Reduce X to the selected features and then return the score of the
underlying estimator.
Parameters
X : array of shape [n_samples, n_features]
The input samples.
y : array of shape [n_samples]
The target values.'

@if_delegate_has_method(delegate='estimator')
def score(self, X, y):

    check_is_fitted(self, 'estimator_')
    return self.estimator_.score(self.transform(X), y)
