Calcualate a trimean using Harrell-Davis quantile estimator (default = Tukey's trimean)

trimean.R

trimean <- function(x, p = 0.25, a = 0.25){
  q <- Hmisc::hdquantile(x, probs = c(p, 0.5, 1-p), na.rm = TRUE, names = FALSE)
  tm <- (a * q[1]) + ((1 - (2 * p)) * q[2]) + (a * q[3])
  return(tm)
}