import statsmodels.api as sm 

x = sm.add_constant(x)    # 添加截距项
modelDfs = sm.OLS(y, x) 
result = modelDfs.fit()
result.summary()