import statsmodels.api as sm

X_stats = sm.add_constant( X ) 
stats = sm.OLS( y, X_stats ).fit()
print( stats.summary() )