pseudo_res function

Calculate pseudo-residuals