xfit = np.linspace(0, 10, 1000) yfit = model.predict(xfit[:,np.newaxis])
plt.scatter(x, y) plt.plot(xfit, yfit, 'r');
model.coef_[0]
3.0154779977316206
model.intercept_
-3.9468162623147478