ENH/REF add fit_regularized to Poisson, returns results, with tests for params, bse