Resid function

Compute residuals