Linear regression, with residuals plot