计算相关系数
import numpy
temp=numpy.corrcoef(X,Y)
r=temp[0][1]
reference
https://jingyan.baidu.com/article/60ccbceba822f464cab1972a.html
tips
1.相关系数对异常值敏感
ECratio[120,3]=-99,第一个相关系数为0.051
ECratio[120,3]=0,第一个相关系数为0.074
回归系数 t检验
python 线性回归 统计检验 p值_线性回归模型理论&案例
使用 sklearn
https://scikit-learn.org/stable/modules/linear_model.html#ordinary-least-squares
中文解释
使用statsmodels
基于表格数据
1.
https://blog.csdn.net/weixin_44835596/article/details/89716684
2.
https://www.cnblogs.com/dhcn/p/13474762.html
基于随机生成的数据
https://www.cnblogs.com/youcans/archive/2021/05/06/14734197.html
Python数模笔记-StatsModels 统计回归(2)线性回归
https://www.cnblogs.com/youcans/archive/2021/05/06/14734197.html
import statsmodels.api as sm
xin=SRratio
yin=ECratio[:,1]
# 创建线性回归最小二乘法模型
x = np.array(xin).reshape(-1,1)
y = np.array(yin).reshape(-1,1)
x = np.asarray(x,'int64')
y = np.asarray(y,'int64')
xx= sm.add_constant(x)
#xx=np.hstack((np.ones((65,1)),x))
model = sm.OLS(y,xx)
results = model.fit()
results.summary()