This would be a great feature and would make the method more consistent with the other linear regression estimators.

I think the only fix that is needed would be to add these lines to the BayesianRidge.fit() method, and update a couple kwargs and test examples?