computeDR function

Compute deviance residuals