Add the results for the preliminary S-estimate to the return object of sparseMM()
Correct error:
initial lasso estimate needs to be computed on (1-BDP)*n observations (and not on BDP*n observations)
Correct error: Reweighting step needs to be computed on observations with |residuals|<2*scale(residuals)