plt.plot(fitted_vals, residuals, 'o') # plotting predictions from #fit model vs residuals