Instatistics, theHodges–Lehmann estimator is arobust andnonparametricestimator of a population'slocation parameter. For populations that are symmetric about onemedian, such as the Gaussian ornormal distribution or theStudentt-distribution, the Hodges–Lehmann estimator is aconsistent and median-unbiased estimate of the population median. For non-symmetric populations, the Hodges–Lehmann estimator estimates the "pseudo–median", which is closely related to the population median.
The Hodges–Lehmann estimator was proposed originally for estimating the location parameter of one-dimensional populations, but it has been used for many more purposes. It has been used to estimate thedifferences between the members of two populations. It has been generalized from univariate populations tomultivariate populations, which produce samples ofvectors.
It is based on theWilcoxon signed-rank statistic. In statistical theory, it was an early example of arank-based estimator, an important class of estimators both in nonparametric statistics and in robust statistics. The Hodges–Lehmann estimator was proposed in 1963 independently byPranab Kumar Sen and byJoseph Hodges andErich Lehmann, and so it is also called the "Hodges–Lehmann–Sen estimator".[1]
In the simplest case, the "Hodges–Lehmann" statistic estimates the location parameter for a univariate population.[2][3] Its computation can be described quickly. For a dataset withn measurements, the set of all possible two-element subsets of it such that ≤ (i.e. specifically including self-pairs; many secondary sources incorrectly omit this detail), which set hasn(n + 1)/2 elements. For each such subset, the mean is computed; finally, the median of thesen(n + 1)/2 averages is defined to be the Hodges–Lehmann estimator of location.
The two-sample Hodges–Lehmann statistic is an estimate of a location-shift typedifference between two populations. For two sets of data withm andn observations, the set of two-element sets made of them is their Cartesian product, which containsm × n pairs of points (one from each set); each such pair defines one difference of values. The Hodges–Lehmann statistic is themedian of them × n differences.[4]
In the general case the Hodges-Lehmann statistic estimates the population'spseudomedian,[5] alocation parameter that is closely related to themedian. The difference between the median and pseudo-median is relatively small, and so this distinction is neglected in elementary discussions. Like thespatial median,[6] the pseudo–median is well defined for all distributions of random variables having dimension two or greater; for one-dimensional distributions, there exists some pseudo–median, which need not be unique, however. Like the median, the pseudo–median is defined for even heavy–tailed distributions that lack any (finite)mean.[7]
For a population that is symmetric, the Hodges–Lehmann statistic also estimates the population's median. It is a robust statistic that has abreakdown point of 0.29, which means that the statistic remains bounded even if nearly 30 percent of the data have been contaminated. This robustness is an important advantage over the sample mean, which has a zero breakdown point, being proportional to any single observation and so liable to being misled by even oneoutlier. Thesample median is even more robust, having a breakdown point of 0.50.[8] The Hodges–Lehmann estimator is much better than the sample mean when estimating mixtures of normal distributions, also.[9]
For symmetric distributions, the Hodges–Lehmann statistic sometimes has greaterefficiency at estimating the center of symmetry (population median) than does the sample median. For the normal distribution, the Hodges-Lehmann statistic is nearly as efficient as the sample mean. For theCauchy distribution (Student t-distribution with one degree of freedom), the Hodges-Lehmann is infinitely more efficient than the sample mean, which is not a consistent estimator of the median,[8] but it is not more efficient than the median in that instance.
The one-sample Hodges–Lehmann statistic need not estimate any population mean, which for many distributions does not exist. The two-sample Hodges–Lehmann estimator need not estimate the difference of two means or the difference of two (pseudo-)medians; rather, it estimates the median of the distribution of the difference between pairs of random–variables drawn respectively from the two populations.[4]
The Hodges–Lehmannunivariate statistics have several generalizations inmultivariate statistics:[10]