一、标准线性回归
在线性回归中我们要求的参数为:
详细的推导可以参见:http://blog.csdn.net/weiyongle1996/article/details/73727505
所以代码实现主要就是实现上式,python代码如下:
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
# implement stand regress
def standRegress(xArr,yArr):
# 将数组转换为矩阵
xMat = np.mat(xArr)
yMat = np.mat(yArr)
xTx = xMat.T * xMat # 计算xTx的
if np.linalg.det(xTx) == 0.0:
print('xTx不能求逆矩阵')
return
theta = xTx.I * (xMat.T * yMat)
yHat = xMat*theta
return yHat
# import data
ex0 = np.loadtxt('ex0.txt',delimiter='\t')
# deal with data
xArr = []
yArr = []
for data in ex