Fitting a robust linear model