import statsmodels.api as sm
X_stats = sm.add_constant(X) 
stats = sm.OLS(y, X_stats).fit()
print(separator)
print( stats.summary() )