step_with_na function

Apply step(), taking into account missing values