be_zeroinfl function

conduct backward stepwise variable elimination for zero inflated count regression