Линия наилучшего соответствия в линейной регрессии

fx = np.array([xs.min(), xs.max()])
fy = res.intercept + res.slope * fx
plt.plot(fx, fy, '-')
josh.ipynb