which.median.R

which.median <- function (x, ties.method = c('first', 'last', 'random', 'all'), ...) {
  med <- median(x, ...)
  dif <- abs(med - x)
  ind <- which(dif == min(dif))
  if (length(ind) > 1 & ties.method[1] != 'all') {
    ind <- if (ties.method[1] == 'first') ind[1] else {if (ties.method[1] == 'last') tail(ind, 1) else sample(ind, 1)}
  }
  return (ind)
}