[MRG+1] TheilSen robust linear regression

A multiple linear Theil-Sen regression for the Scikit-Learn toolbox. The implementation is based on the algorithm from the paper "Theil-Sen Estimators in a Multiple Linear Regression Model" of Xin Dang, Hanxiang Peng, Xueqin Wang and Heping Zhang. It is parallelized with the help of joblib.
On a personal note, I think that the popular Theil-Sen regression would be a nice addition to Scikit-Learn.
I am looking forward to your feedback. 

Florian
dumb question, why not materialising the indices in _lstsq? Is the indices list huge? 

Sorry if this is obvious. Could you briefly explained the difference between both index generation strategy? Is it wanted to sample with replacement?