sm.OLS(y, sm.add_constant(x)).fit()
