Theil Sen regression involves fitting multiple regression models on subsets of the training data and combining the coefficients together in the end.

The algorithm calculates least square solutions on subsets with size of the samples in . In a final step, the spatial median (or L1 median) is calculated of all least square solutions.

Sources: