pMedian function

pMedian

Pseudomedian

pMedian(x, na.rm = FALSE)

Arguments

• x: a numeric vector
• na.rm: set to TRUE to exclude NAs before computing the pseudomedian

Returns

a scalar numeric value

Details

Uses fast Fortran code to compute the pseudomedian of a numeric vector. The pseudomedian is the median of all possible midpoints of two observations. The pseudomedian is also called the Hodges-Lehmann one-sample estimator. The Fortran code is was originally from JF Monahan, and was converted to C++ in the DescTools package. It has been converted to Fortran 2018 here.

Examples

x <- c(1:4, 10000)
pMedian(x)
# Compare with brute force calculation and with wilcox.test
w <- outer(x, x, '+')
median(w[lower.tri(w, diag=TRUE)]) / 2
wilcox.test(x, conf.int=TRUE)

See Also

https://dl.acm.org/toc/toms/1984/10/3/, https://www4.stat.ncsu.edu/~monahan/jul10/