weighted_median: Weighted Median
In gamboostLSS: Boosting Methods for GAMLSS

Description Usage Arguments Details References See Also Examples

Function to compute the weighted median.

1	weighted.median(x, w = 1, na.rm = FALSE)

`x`	a numeric vector containing the values whose median is to be computed.
`w`	weights that are used to compute the median. This can be either a single value (which will be used as weight for all observations) or a numeric vector of the same length as `x`.
`na.rm`	logical. Should `NA`s (from weights `w` and/or data `x`) be removed?

The weighted median is computed as the value where the cumulative relative weights (relative to the sum of all weights) crosses 0.5.

This function is used in the stabilization of the negative gradient via the meadian absolute deviation (MAD). For details see Hofner et al (2014).

B. Hofner, A. Mayr, M. Schmid (2014). gamboostLSS: An R Package for Model Building and Variable Selection in the GAMLSS Framework. Technical Report, arXiv:1407.1774.

glmboostLSS, gamboostLSS and blackboostLSS for fitting of GAMLSS where the standardization is explained in more detail.

## compute the weighted median with case weights
x <- c(1, 2, 3, 4)
w <- c(0, 1, 2, 3)
weighted.median(x, w)

## compute the weighted median with arbitrary weights
x <- rnorm(100)
w <- runif(100)
weighted.median(x, w)