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 [：，1]我使用了A [：，0]，这就是结果：

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 param只需要是一个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[:, 0]。

您可能会发现在开始时将数据拆分为X和y更容易 - 这可能会使事情更加清晰。