smooth_residuals function

Smooth residuals after model fit