TODO:: Write documentation
.residual( x, y, n, pars, f, scale = 0.07, family = c("normal", "Cauchy"), trunc = FALSE )
A numeric containing the estimated residuals for the model fit
numeric