R/weighted.mad.R

In descriptio: Descriptive Statistical Analysis

weighted.mad <- function(x, weights = NULL, na.rm = FALSE) {
  if(is.null(weights)) weights <- rep(1, length(x))
  if(any(is.na(weights))) stop("There are empty values in weights.")
  if(na.rm) {
    complete <- !is.na(x)
    x <- x[complete]
    weights <- weights[complete]
  } else {
    if(any(is.na(x))) stop("There are empty values in x. \nPlease consider transforming your data (filtering, recoding, imputation, etc.) or set na.rm to TRUE.")
  }
  if(length(x)==1) {
    mad <- 0 
  } else {
    med <- weighted.quantile(x=x, weights=weights, probs = .5)
    ad <- abs(x-med)
    mad <- weighted.quantile(x=ad, weights=weights, probs = .5)
  }
  return(mad)
}