python scikit線性回歸怪異結果

Question

我是python的新手。

使用matplotlib繪制線性回歸的結果。

我已經嘗試了一些基本數據並且可以正常工作，但是當我嘗試實際數據時，回歸線完全錯誤。 我認為我在fit（）或predict（）函數中做錯了什么。

這是代碼：

import matplotlib.pyplot as plt
from sklearn import linear_model
import scipy
import numpy as np
regr=linear_model.LinearRegression()
A=[[69977, 4412], [118672, 4093], [127393, 12324], [226158, 15453], [247883, 8924], [228057, 6568], [350119, 4040], [197808, 6793], [205989, 8471], [10666, 632], [38746, 1853], [12779, 611], [38570, 1091], [38570, 1091], [95686, 8752], [118025, 17620], [79164, 13335], [83051, 1846], [4177, 93], [29515, 1973], [75671, 5070], [10077, 184], [78975, 4374], [187730, 17133], [61558, 2521], [34705, 1725], [206514, 10548], [13563, 1734], [134931, 7117], [72527, 6551], [16014, 310], [20619, 403], [21977, 437], [20204, 258], [20406, 224], [20551, 375], [38251, 723], [20416, 374], [21125, 429], [20405, 235], [20042, 431], [20016, 366], [19702, 200], [20335, 420], [21200, 494], [22667, 487], [20393, 405], [20732, 414], [20602, 393], [111705, 7623], [112159, 5982], [6750, 497], [59624, 418], [111468, 10209], [40057, 1484], [435, 0], [498848, 17053], [26585, 1390], [75170, 3883], [139146, 3540], [84931, 7214], [19144, 3125], [31144, 2861], [66573, 818], [114253, 4155], [15421, 2094], [307497, 5110], [484904, 10273], [373476, 36365], [128152, 10920], [517285, 106315], [453483, 10054], [270763, 17542], [9068, 362], [61992, 1608], [35791, 1747], [131215, 6227], [4314, 191], [16316, 2650], [72791, 2077], [47008, 4656], [10853, 1346], [66708, 4855], [214736, 11334], [46493, 4236], [23042, 737], [335941, 11177], [65167, 2433], [94913, 7523], [454738, 12335]]
#my data are selected from a Mysql DB  and stored in np array like this one above.



regr.fit(A,A[:,1]) 
plt.scatter(A[:,0],A[:,1], color='black')
plt.plot(A[:,1],regr.predict(A), color='blue',linewidth=3)
plt.show()

所需的是使用A的第一列和第二列的數據的回歸線。 結果如下：

我知道異常值的存在會嚴重影響輸出，但是我嘗試使用其他收費進行回歸，回歸線距離點所在的區域更近，所以我肯定會丟失一些東西。

謝謝。

編輯1：按照建議，我再次嘗試僅更改plot（）參數。 我使用A [：，0]而不是A [：，1]，這是結果：

scikit-learn.org/stable/modules/linear_model.html上的一個簡單示例看起來像我的。 我不需要預測，所以我沒有在訓練和測試集中切出我的數據……也許是我誤解了“ X，y”的含義，但是再次看一下鏈接中的示例，它看起來像我的。

編輯2：終於成功了。

X=A[:,0]
X=X[:,np.newaxis]
regr=linear_model.LinearRegression()
regr.fit(X,A[:,1])
plt.plot(X,regr.predict(X))

X參數僅需要是2 Dim數組。 編輯1中的示例確實誤導了我:(。

Answer 1

您似乎在訓練數據中包括目標值A[:, 1] 。 擬合命令的格式為regr.fit(X, y) 。

您似乎也對這一行有疑問：

plt.plot(A[:,1],regr.predict(A), color='blue',linewidth=3)

我認為，如果要針對預測變量值繪制預測，應該將A[:, 1]替換為A[:, 1] A[:, 0] 。

您可能會發現更容易在開始時將數據分為X和y ，這可能會使事情變得更清楚。