step_optim function

Forward and backward model selection