linalg
使用Numpy提供的一个线性代数的库
1.1、linalg.det(mat) , 计算矩阵的行列式# 计算回归系数w
def standRegres(xArr,yArr):
'''
计算回归系数
:param xArr: x数据集
:param yArr: y数据集
:return: 回归系数
'''
xMat = mat(xArr)
yMat = mat(yArr).T # 由于yArr是一个列表, 而yMat需要的是一个列向量, 所以需要转置
xTx = xMat.T*xMat
# 前提条件, xTx不可逆 行列式不为0
if linalg.det(xTx) == 0.0:
print("This matrix is singular, cannot do inverse")
return
ws = xTx.I * (xMat.T*yMat)
return ws