blr_step_aic_both function

Stepwise AIC selection