U-statistic

In statistical theory, a U-statistic is a class of statistics that is especially important in estimation theory; the letter "U" stands for unbiased. In elementary statistics, U-statistics arise naturally in producing minimum-variance unbiased estimators.

The theory of U-statistics allows a minimum-variance unbiased estimator to be derived from each unbiased estimator of an estimable parameter (alternatively, statistical functional) for large classes of probability distributions.[1][2] An estimable parameter is a measurable function of the population's cumulative probability distribution: For example, for every probability distribution, the population median is an estimable parameter. The theory of U-statistics applies to general classes of probability distributions.

Many statistics originally derived for particular parametric families have been recognized as U-statistics for general distributions. In non-parametric statistics, the theory of U-statistics is used to establish for statistical procedures (such as estimators and tests) and estimators relating to the asymptotic normality and to the variance (in finite samples) of such quantities.[3] The theory has been used to study more general statistics as well as stochastic processes, such as random graphs.[4][5][6]

Suppose that a problem involves independent and identically-distributed random variables and that estimation of a certain parameter is required. Suppose that a simple unbiased estimate can be constructed based on only a few observations: this defines the basic estimator based on a given number of observations. For example, a single observation is itself an unbiased estimate of the mean and a pair of observations can be used to derive an unbiased estimate of the variance. The U-statistic based on this estimator is defined as the average (across all combinatorial selections of the given size from the full set of observations) of the basic estimator applied to the sub-samples.

Sen (1992) provides a review of the paper by Wassily Hoeffding (1948), which introduced U-statistics and set out the theory relating to them, and in doing so Sen outlines the importance U-statistics have in statistical theory. Sen says[7] "The impact of Hoeffding (1948) is overwhelming at the present time and is very likely to continue in the years to come". Note that the theory of U-statistics is not limited to[8] the case of independent and identically-distributed random variables or to scalar random-variables.[9]

Definition

The term U-statistic, due to Hoeffding (1948), is defined as follows.

Let f\colon R^r\to R be a real-valued or complex-valued function of r variables. For each n\ge r the associated U-statistic f_n\colon R^n \to R is equal to the average over ordered samples \varphi(1),\ldots, \varphi(r) of size r of the sample values f(x_\varphi). In other words, f_n(x_1,\ldots, x_n) = \operatorname{ave} f(x_{\varphi(1)},\ldots, x_{\varphi(r)}), the average being taken over distinct ordered samples of size r taken from \{1,\ldots, n\}. Each U-statistic f_n(x_1,\ldots, x_n) is necessarily a symmetric function.

U-statistics are very natural in statistical work, particularly in Hoeffding's context of independent and identically-distributed random variables, or more generally for exchangeable sequences, such as in simple random sampling from a finite population, where the defining property is termed `inheritance on the average'.

Fisher's k-statistics and Tukey's polykays are examples of homogeneous polynomial U-statistics (Fisher, 1929; Tukey, 1950). For a simple random sample φ of size n taken from a population of size N, the U-statistic has the property that the average over sample values ƒn(xφ) is exactly equal to the population value ƒN(x).

Examples

Some examples: If f(x) = x the U-statistic f_n(x) = \bar x_n = (x_1 + \cdots + x_n)/n is the sample mean.

If f(x_1, x_2) = |x_1 - x_2|, the U-statistic is the mean pairwise deviation f_n(x_1,\ldots, x_n) = \sum_{i\neq j} |x_i - x_j| / (n(n-1)), defined for n\ge 2.

If f(x_1, x_2) = (x_1 - x_2)^2/2, the U-statistic is the sample variance f_n(x) = \sum(x_i - \bar x_n)^2/(n-1) with divisor n-1, defined for n\ge 2.

The third k-statistic k_{3,n}(x) = \sum(x_i - \bar x_n)^3 n/((n-1)(n-2)), the sample skewness defined for n\ge 3, is a U-statistic.

The following case highlights an important point. If f(x_1, x_2, x_3) is the median of three values, f_n(x_1,\ldots, x_n) is not the median of n values. However, it is a minimum variance unbiased estimate of the expected value of the median of three values and in this application of the theory it is the population parameter defined as "the expected value of the median of three values" which is being estimated, not the median of the population. Similar estimates play a central role where the parameters of a family of probability distributions are being estimated by probability weighted moments or L-moments.

See also

Notes

  1. Cox & Hinkley (1974),p. 200, p. 258
  2. Hoeffding (1948), between Eq's(4.3),(4.4)
  3. Sen (1992)
  4. Page 508 in Koroljuk, V. S.; Borovskich, Yu. V. (1994). Theory of U-statistics. Mathematics and its Applications 273 (Translated by P. V. Malyshev and D. V. Malyshev from the 1989 Russian original ed.). Dordrecht: Kluwer Academic Publishers Group. pp. x+552. ISBN 0-7923-2608-3. MR 1472486.
  5. Pages 381–382 in Borovskikh, Yu. V. (1996). U-statistics in Banach spaces. Utrecht: VSP. pp. xii+420. ISBN 90-6764-200-2. MR 1419498.
  6. Page xii in Kwapień, Stanisƚaw; Woyczyński, Wojbor A. (1992). Random series and stochastic integrals: Single and multiple. Probability and its Applications. Boston, MA: Birkhäuser Boston, Inc. pp. xvi+360. ISBN 0-8176-3572-6. MR 1167198.
  7. Sen (1992) p. 307
  8. Sen (1992), p306
  9. Borovskikh's last chapter discusses U-statistics for exchangeable random elements taking values in a vector space (separable Banach space).

References

This article is issued from Wikipedia - version of the Wednesday, November 11, 2015. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.