In computation of initial estimate in sparseS() add that the following:
Now initial estimate is computed as a lasso estimate on the observations that correspond to the smallest BDP*n residuals from the lasso estimate on three observations

(This mimics the computation in sparseLTS() and, thus, hopefully gives more stable results)