# fit the model
results = ols_model.fit()