Analytics, Data Science, Statistics, Glm, Poisson Regression.
remaining_candidates and current_score == best_new_score: scores_with_candidates = [] for candidate in remaining_candidates: formula = f"{response} ~ {' + '.join(selected_candidates + [candidate])}" model = smf.glm(formula, data=data, family=sm.families.Poisson()) results = model.fit() score =