# final formula
 formula = "{} ~ {} + 1".format(response, ' + '.join(possible_factors))
 # model object
 model = sm.ols(formula, data).fit()
 return model