# ols = ordinary least squares
ols_model = sm.ols(formula = 'y ~ x', data=syn_data)