sklearn和自己的线性回归的实现

小锟

1. import pandas as pd
   
2. import matplotlib.pyplot as plt
   
3. import numpy as np
   
4. import os
   
5. os.chdir(r'C:\Users\Administrator\Desktop')
   
6. 
   
7. data = pd.read_table('1.txt',sep = '\s*')
   
8. 
   
9. plt.scatter(data['Height'],data['Weight'])
   
10. #进行标准化
    
11. a = data['Height'].values
    
12. a_ = (a - a.mean()) / a.std()
    
13. b = data['Weight'].values
    
14. b_ = (b - b.mean()) / b.std()
    
15. 
    
16. ###Python线性回归的两种实现
    
17. import lines
    
18. theat = np.array([0,0])
    
19. y = lines.function1(data['Height'].values[:,np.newaxis],data['Weight'].values[:,np.newaxis].T)
    
20. h = lines.function2(a_[:,np.newaxis],
    
21.                       b_[:,np.newaxis].T,
    
22.                       rate=0.01,theat=theat,num=19,minnum=0.000001)
    
23. ##进行逆标准化
    
24. h[0],h[1] = h[0] * b.std() - h[1] * b.std() / a.std() * a.mean() + b.mean(),\
    
25.             h[1] * b.std() / a.std()
    
26. print y , h
    
27. 
    
28. ##和sklearn库进行比较
    
29. from sklearn import linear_model
    
30. linreg = linear_model.LinearRegression()
    
31. linreg.fit(data['Height'].values[:,np.newaxis],data['Weight'].values[:,np.newaxis])
    
32. print linreg.coef_
    
33. print linreg.intercept_
    
34. linre = linear_model.LinearRegression()
    
35. linre.fit(a_[:,np.newaxis],
    
36.                       b_[:,np.newaxis])
    
37. print linreg.coef_
    
38. print linreg.intercept_
    
39. y_ = linreg.predict(a[:,np.newaxis])
    
40. plt.plot(a,y_)
    
41. 
    
42. 
    
43. plt.subplots()
    
44. plt.scatter(a_,b_)
    
45. y__ = linre.predict(a_[:,np.newaxis])
    
46. plt.plot(a_,y__)
    
47. plt.show()

复制代码

最小二乘[[-143.02691844]
[   3.89903027]]
梯度下降 [[-142.99965677]
[   3.89859294]]
sklearn库
[[ 3.89903027]]
[-143.02691844]
[[ 3.89903027]]
[-143.02691844]

数据在附件中

shooter110

这是线性回归吧，这不是向量自回归模型

小锟

shooter110 发表于 2017-6-8 22:26
这是线性回归吧，这不是向量自回归模型

这不是var模型啊，var用Eviews好一点

Z101

你好，请问为什么报错了？
FileNotFoundError: [WinError 3] 系统找不到指定的路径。: 'C:\Users\Administrator\Desktop'