黑松山资源网 Design By www.paidiu.com
本文实例讲述了Python实现的线性回归算法。分享给大家供大家参考,具体如下:
用python实现线性回归
Using Python to Implement Line Regression Algorithm
小菜鸟记录学习过程
代码:
#encoding:utf-8 """ Author: njulpy Version: 1.0 Data: 2018/04/09 Project: Using Python to Implement LineRegression Algorithm """ import numpy as np import pandas as pd from numpy.linalg import inv from numpy import dot from sklearn.model_selection import train_test_split import matplotlib.pyplot as plt from sklearn import linear_model # 最小二乘法 def lms(x_train,y_train,x_test): theta_n = dot(dot(inv(dot(x_train.T, x_train)), x_train.T), y_train) # theta = (X'X)^(-1)X'Y #print(theta_n) y_pre = dot(x_test,theta_n) mse = np.average((y_test-y_pre)**2) #print(len(y_pre)) #print(mse) return theta_n,y_pre,mse #梯度下降算法 def train(x_train, y_train, num, alpha,m, n): beta = np.ones(n) for i in range(num): h = np.dot(x_train, beta) # 计算预测值 error = h - y_train.T # 计算预测值与训练集的差值 delt = 2*alpha * np.dot(error, x_train)/m # 计算参数的梯度变化值 beta = beta - delt #print('error', error) return beta if __name__ == "__main__": iris = pd.read_csv('iris.csv') iris['Bias'] = float(1) x = iris[['Sepal.Width', 'Petal.Length', 'Petal.Width', 'Bias']] y = iris['Sepal.Length'] x_train, x_test, y_train, y_test = train_test_split(x, y, test_size=0.2, random_state=5) t = np.arange(len(x_test)) m, n = np.shape(x_train) # Leastsquare theta_n, y_pre, mse = lms(x_train, y_train, x_test) # plt.plot(t, y_test, label='Test') # plt.plot(t, y_pre, label='Predict') # plt.show() # GradientDescent beta = train(x_train, y_train, 1000, 0.001, m, n) y_predict = np.dot(x_test, beta.T) # plt.plot(t, y_predict) # plt.plot(t, y_test) # plt.show() # sklearn regr = linear_model.LinearRegression() regr.fit(x_train, y_train) y_p = regr.predict(x_test) print(regr.coef_,theta_n,beta) l1,=plt.plot(t, y_predict) l2,=plt.plot(t, y_p) l3,=plt.plot(t, y_pre) l4,=plt.plot(t, y_test) plt.legend(handles=[l1, l2,l3,l4 ], labels=['GradientDescent', 'sklearn','Leastsquare','True'], loc='best') plt.show()
输出结果
sklearn: [ 0.65368836 0.70955523 -0.54193454 0. ]
LeastSquare: [ 0.65368836 0.70955523 -0.54193454 1.84603897]
GradientDescent: [ 0.98359285 0.29325906 0.60084232 1.006859 ]
附:上述示例中的iris.csv文件点击此处本站下载。
更多关于Python相关内容感兴趣的读者可查看本站专题:《Python数学运算技巧总结》、《Python数据结构与算法教程》、《Python函数使用技巧总结》、《Python字符串操作技巧汇总》及《Python入门与进阶经典教程》
希望本文所述对大家Python程序设计有所帮助。
黑松山资源网 Design By www.paidiu.com
广告合作:本站广告合作请联系QQ:858582 申请时备注:广告合作(否则不回)
免责声明:本站资源来自互联网收集,仅供用于学习和交流,请遵循相关法律法规,本站一切资源不代表本站立场,如有侵权、后门、不妥请联系本站删除!
免责声明:本站资源来自互联网收集,仅供用于学习和交流,请遵循相关法律法规,本站一切资源不代表本站立场,如有侵权、后门、不妥请联系本站删除!
黑松山资源网 Design By www.paidiu.com
暂无评论...
更新日志
2024年10月07日
2024年10月07日
- 《丁当 20首重量级歌曲 Fu Good 下一站天后 2CD》[WAV/分轨][650MB]
- 歌莉雅.2012-My.Voice.Story【环星】【FLAC分轨】
- 群星.1993-一曲成名·青春无悔【飞碟】【WAV+CUE】
- 刘德华.2000-爱无知(金装版)【NEWMELODY】【WAV+CUE】
- 许巍《时光漫步》1CD[FLAC/分轨][804.3MB]
- 群星《2024第一季度百度歌曲排行抖音榜top100》1CD[FLAC/分轨][2.3GB]
- 张柏芝《全新经验》[FLAC/分轨][450MB]
- 唐朝.2008-浪漫骑士【音乐家】【WAV+CUE】
- 张芸京.2016-失败的高歌【泡耳音乐】【WAV+CUE】
- 群星.1991-音乐工厂1·皇后大道东【音乐工厂】【WAV+CUE】
- 群星.1992-音乐工厂2·首都【音乐工厂】【WAV+CUE】
- 群星.1994-音乐工厂3·儿童乐园【音乐工厂】【WAV+CUE】
- 胡鸿钧.2021-ex:CHANGE【星梦娱乐】【WAV+CUE】
- 班得瑞原装进口《第一张新世纪专辑:仙境》1CD[APE/CUE分轨][292.3MB]
- 伍佰《摇滚教父 伍佰黄金精选 2CD》[WAV+CUE][990MB]