explore function

Run first step of model fitting to find good penalization interval