def fit(self, X, y):
X, y, X_mean, y_mean, X_std = self._center_data(X, y, self.fit_intercept, self.normalize, self.copy_X)<EOL>self.coef_, self.shrink_coef_, self.rss_ = non_negative_garotte(X, y, self.alpha)<EOL>self._set_intercept(X_mean, y_mean, X_std)<EOL>
:param X:
:param y:
:return: