add zero.add and zero.keep.margins argument to all fitting functions to address zero cells in bivariate frequency tables