U-statistic explained

In statistical theory, a U-statistic is a class of statistics defined as the average over the application of a given function applied to all tuples of a fixed size. 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.

History

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.

Pranab K. 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

K

be either the real or complex numbers, and let

f\colon(Kd)r\toK

be a

K

-valued function of

r

d

-dimensional variables.For each

n\ger

the associated U-statistic

fn\colon(Kd)n\toK

is defined to be the average of the values
f(x
i1

,...c,

x
ir

)

over the set

Ir,

of

r

-tuples of indices from

\{1,2,...c,n\}

with distinct entries.Formally,

fn(x1,...c,xn)=

1
r-1
\prod(n-i)
i=0
\sum
(i1,...c,ir)\inIr,
f(x
i1

,...c,

x
ir

)

.In particular, if

f

is symmetric the above is simplified to

fn(x1,...c,xn)=

1
\binom{n

{r}}

\sum
(i1,...c,ir)\inJr,
f(x
i1

,...c,

x
ir

)

,where now

Jr,

denotes the subset of

Ir,

of increasing tuples.

Each U-statistic

fn

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() is exactly equal to the population value ƒN(x).

Examples

Some examples:If

f(x)=x

the U-statistic

fn(x)=\barxn=(x1++xn)/n

is the sample mean.

If

f(x1,x2)=|x1-x2|

, the U-statistic is the mean pairwise deviation

fn(x1,\ldots,xn)=2/(n(n-1))\sumi>|xi-xj|

, defined for

n\ge2

.

If

f(x1,x2)=(x1-

2/2
x
2)
, the U-statistic is the sample variance

fn(x)=\sum(xi-\bar

2/(n-1)
x
n)
with divisor

n-1

, defined for

n\ge2

.

The third

k

-statistic

k3,n(x)=\sum(xi-\bar

3
x
n)

n/((n-1)(n-2))

,the sample skewness defined for

n\ge3

,is a U-statistic.

The following case highlights an important point. If

f(x1,x2,x3)

is the median of three values,

fn(x1,\ldots,xn)

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, 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

References

Notes and References

  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 Book: Koroljuk . V. S. . Borovskich . Yu. V. . Theory of U-statistics . Translated by P. V. Malyshev and D. V. Malyshev from the 1989 Russian original . Mathematics and its Applications . 273 . Kluwer Academic Publishers Group . Dordrecht . 1994 . x+552 . 0-7923-2608-3 . 1472486 .
  5. Pages 381–382 in Book: Borovskikh, Yu. V. . U-statistics in Banach spaces . VSP . Utrecht . 1996 . xii+420 . 90-6764-200-2 . 1419498 .
  6. Page xii in Book: Kwapień . Stanisƚaw . Woyczyński . Wojbor A. . Random series and stochastic integrals: Single and multiple . Probability and its Applications . Birkhäuser Boston, Inc.. Boston, MA . 1992 . xvi+360 . 0-8176-3572-6 . 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).