REF: NegativeBinomial fit_regularized, try regularized Poisson for start_params