R/scaleResidualsByMAD.R

scaleResidualsByMAD <- function (residuals) {
  if (length(residuals) / 2 > sum(residuals == 0)) {
    residuals <- abs(residuals) / (1.4826 * median(abs(residuals)))
  } else {
    residuals <- abs(residuals) / (1.4826 * median(abs(residuals[residuals != 0])))
  } # all scaled residuals are non-negative
  return(residuals)
}

Try the crmReg package in your browser

Any scripts or data that you put into this service are public.

crmReg documentation built on June 8, 2025, 10:03 a.m.