línea de mejor ajuste en regresión lineal

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