#jsDisabledContent { display:none; } My Account |  Register |  Help

# Energy distance

Article Id: WHEBN0030316300
Reproduction Date:

 Title: Energy distance Author: World Heritage Encyclopedia Language: English Subject: Collection: Publisher: World Heritage Encyclopedia Publication Date:

### Energy distance

Energy distance is a statistical distance between probability distributions. If X and Y are independent random vectors in Rd with cumulative distribution functions F and G respectively, then the energy distance between the distributions F and G is defined to be the square root of

D^2(F, G) = 2\mathbb E\|X - Y\| - \mathbb E\|X - X'\| - \mathbb E\|Y - Y'\| \geq 0,

where X, X' are independent and identically distributed (iid), Y, Y' are iid, \mathbb E is expected value, and || . || denotes the length of a vector. Energy distance satisfies all axioms of a metric thus energy distance characterizes the equality of distributions: D(F,G) = 0 if and only if F = G. Energy distance for statistical applications was introduced in 1985 by Gábor J. Székely, who proved that for real-valued random variables this distance is exactly twice Harald Cramér's distance:[1]

\int_{-\infty}^\infty (F(x) - G(x))^2 \, dx .

For a simple proof of this equivalence, see Székely and Rizzo (2005).[2] In higher dimensions, however, the two distances are different because the energy distance is rotation invariant while Cramér's distance is not. (Notice that Cramér's distance is not the same as the distribution-free Cramer-von-Mises criterion.)

## Contents

• Generalization to metric spaces 1
• Energy statistics 2
• Testing for equal distributions 2.1
• Goodness-of-fit 2.2
• Applications 3
• References 4

## Generalization to metric spaces

One can generalize the notion of energy distance to probability distributions on metric spaces. Let (M, d) be a metric space with its Borel sigma algebra \mathcal{B} (M). Let \mathcal{P} (M) denote the collection of all probability measures on the measurable space (M, \mathcal{B} (M)). If μ and ν are probability measures in \mathcal{P} (M), then the energy-distance D of μ and ν can be defined as the square root of

D^2(\mu, \nu)= 2 \mathbb E[d(X,Y)] - \mathbb E[d(X,X')] - \mathbb E[d(Y,Y')] .

This is not necessarily non-negative, however. If (M, d) is a strongly negative definite kernel, then D is a metric, and conversely.[3] This condition is expressed by saying that (M, d) has negative type. Negative type is not sufficient for D to be a metric; the latter condition is expressed by saying that (M, d) has strong negative type. In this situation, the energy distance is zero if and only if X and Y are identically distributed. An example of a metric of negative type but not of strong negative type is the plane with the taxicab metric. All Euclidean spaces and even separable Hilbert spaces have strong negative type.[4]

In the literature on kernel methods for machine learning, these generalized notions of energy distance are studied under the name of maximum mean discrepancy.[5]

## Energy statistics

A related statistical concept, the notion of E-statistic or energy-statistic[6] was introduced by Gábor J. Székely in the 1980s when he was giving colloquium lectures in Budapest, Hungary and at MIT, Yale, and Columbia. This concept is based on the notion of Newton’s potential energy.[7] The idea is to consider statistical observations as heavenly bodies governed by a statistical potential energy which is zero only when an underlying statistical null hypothesis is true. Energy statistics are functions of distances between statistical observations.

Energy distance and E-statistic were considered as N-distances and N-statistic in Zinger A.A., Kakosyan A.V., Klebanov L.B. Characterization of distributions by means of mean values of some statistics in connection with some probability metrics, Stability Problems for Stochastic Models. Moscow, VNIISI, 1989,47-55. (in Russian), English Translation: A characterization of distributions by mean values of statistics and certain probabilistic metrics A. A. Zinger, A. V. Kakosyan, L. B. Klebanov in Journal of Soviet Mathematics (1992). In the same paper there was given a definition of strongly negative definite kernel, and provided a generalization on metric spaces, discussed above. The book [3] gives these results and their applications to statistical testing as well. The book contains also some applications to recovering the measure from its potential.

### Testing for equal distributions

Consider the null hypothesis that two random variables, X and Y, have the same probability distributions: μ = v . For statistical samples from X and Y:

x1,…,xn and y1,…,ym,

the following arithmetic averages of distances are computed between the X and the Y samples:

A:= (1/nm)∑|xi – yj|, B:= (1/n2)∑|xi – xj|, C:= (1/m2)∑|yi – yj|.

The E-statistic of the underlying null hypothesis is defined as follows:

Εn,m(X,Y):= 2A – B – C.

One can prove [2][7] that Εn,m(X,Y) ≥ 0 and that the corresponding population value, E(X,Y):= D(μ,ν), is zero if and only if X and Y have the same distribution (μ=ν). Under this null hypothesis the test statistic

T = \frac{nm}{n+m} E_{n,m}(X,Y)

converges in distribution to a quadratic form of independent standard normal random variables. Under the alternative hypothesis T tends to infinity. This makes it possible to construct a consistent statistical test, the energy test for equal distributions.[8]

The E-coefficient of inhomogeneity can also be introduced. This is always between 0 and 1 and is defined as

H = \frac{D^2(F_X,F_Y)}{ 2 \operatorname{\mathbb E}\|X-Y\|} = \frac{2\mathbb E\|X - Y\| - \mathbb E\|X - X'\| - \mathbb E\|Y - Y'\|} {2 \operatorname{\mathbb E}\|X-Y\|},

where \mathbb E denotes the expected value. H = 0 exactly when X and Y have the same distribution.

### Goodness-of-fit

A multivariate goodness-of-fit measure is defined for distributions in arbitrary dimension (not restricted by sample size). The energy goodness-of-fit statistic is

Q_n = n \left( \frac{2}{n} \sum_{i=1}^n \mathbb E \|x_i - X\|^\alpha - \mathbb E\|X - X'\|^\alpha - \frac{1}{n^2} \sum_{i=1}^n \sum_{j=1}^n \|x_i - x_j\|^\alpha \right),

where X and X' are independent and identically distributed according to the hypothesized distribution, and \alpha \in (0,2). The only required condition is that X has finite \alpha moment under the null hypothesis. Under the null hypothesis \mathbb EQ_n=\mathbb E\|X-X'\|^\alpha, and the asymptotic distribution of Qn is a quadratic form of centered Gaussian random variables. Under an alternative hypothesis, Qn tends to infinity stochastically, and thus determines a statistically consistent test. For most applications the exponent 1 (Euclidean distance) can be applied. The important special case of testing multivariate normality[2] is implemented in the energy package for R. Tests are also developed for heavy tailed distributions such as Pareto (power law), or stable distributions by application of exponents in (0,1).

## Applications

Applications include

Gneiting and Raftery[17] apply energy distance to develop a new and very general type of proper scoring rule for probabilistic predictions, the energy score.
• Robust statistics [18]
• Gene selection [19]
• Microarray data analysis[20]
• Material structure analysis[21]

Applications of energy statistics are implemented in the open source energy package[22] for R.

## References

1. ^ Cramér, H. (1928) On the composition of elementary errors, Skandinavisk Aktuarietidskrift, 11, 141–180.
2. ^ a b c d Reprint
3. ^ a b Klebanov, L. B. (2005) N-distances and their Applications, Karolinum Press, Charles University, Prague.
4. ^ [1]
5. ^ [2]
6. ^ G. J. Szekely and M. L. Rizzo (2013). Energy statistics: statistics based on distances. Journal of Statistical Planning and Inference Volume 143, Issue 8, August 2013, pp. 1249-1272. [3]
7. ^ a b Székely, G.J. (2002) E-statistics: The Energy of Statistical Samples, Technical Report BGSU No 02-16.
8. ^ G. J. Szekely and M. L. Rizzo (2004). Testing for Equal Distributions in High Dimension, InterStat, Nov. (5). Reprint.
9. ^ Székely, G. J. and Rizzo, M. L. (2005) Hierarchical Clustering via Joint Between-Within Distances: Extending Ward's Minimum Variance Method, Journal of Classification, 22(2) 151–183
10. ^ "eprint".
11. ^ M. L. Rizzo and G. J. Székely (2010). DISCO Analysis: A Nonparametric Extension of Analysis of Variance, Annals of Applied Statistics Vol. 4, No. 2, 1034–1055. PDF
12. ^ Szekely, G. J. and Rizzo, M. L. (2004) Testing for Equal Distributions in High Dimension, InterStat, Nov. (5). Reprint.
13. ^ PDF
14. ^ [4] Preprint:TR534.
15. ^ Székely, G. J., Rizzo M. L. and Bakirov, N. K. (2007). "Measuring and testing independence by correlation of distances", The Annals of Statistics, 35, 2769–2794. PDF
16. ^ Székely, G. J. and Rizzo, M. L. (2009). "Brownian distance covariance", The Annals of Applied Statistics, 3/4, 1233–1308. PDF
17. ^ Reprint
18. ^ Klebanov L.B. A class of Probability Metrics and its Statistical Applications, Statistics in Industry and Technology: Statistical Data Analysis, Yadolah Dodge, Ed. Birkhauser, Basel, Boston, Berlin, 2002,241-252.
19. ^ Statistics and Data Analysis, 2006, 50, 12, 3619-3628Rui Hu, Xing Qiu, Galina Glazko, Lev Klebanov, Andrei Yakovlev Detecting intergene correlation changes in microarray analysis: a new approach to gene selection, BMCBioinformatics, Vol.10, 20 (2009), 1-15.
20. ^ Yuanhui Xiao, Robert Frisina, Alexander Gordon, Lev Klebanov, Andrei Yakovlev Multivariate Search for Diferentially Expressed Gene Combinations BMC Bioinformatics, 2004, 5:164; Antoni Almudevar, Lev Klebanov, Xing Qiu, Andrei Yakovlev Utility of correlation measures in analysis of gene expression, In: NeuroRX, 2006, 3, 3, 384-395; Klebanov Lev, Gordon Alexander, Land Hartmut, Yakovlev Andrei A permutation test motivated by microarray data analysis
21. ^ Viktor Benes, Radka Lechnerova, Lev Klebanov, Margarita Slamova, Peter Slama Statistical comparison of the geometry of second-phase particles, Materials Characterization , Vol. 60 (2009 ), 1076 - 1081.
22. ^
This article was sourced from Creative Commons Attribution-ShareAlike License; additional terms may apply. World Heritage Encyclopedia content is assembled from numerous content providers, Open Access Publishing, and in compliance with The Fair Access to Science and Technology Research Act (FASTR), Wikimedia Foundation, Inc., Public Library of Science, The Encyclopedia of Life, Open Book Publishers (OBP), PubMed, U.S. National Library of Medicine, National Center for Biotechnology Information, U.S. National Library of Medicine, National Institutes of Health (NIH), U.S. Department of Health & Human Services, and USA.gov, which sources content from all federal, state, local, tribal, and territorial government publication portals (.gov, .mil, .edu). Funding for USA.gov and content contributors is made possible from the U.S. Congress, E-Government Act of 2002.

Crowd sourced content that is contributed to World Heritage Encyclopedia is peer reviewed and edited by our editorial staff to ensure quality scholarly research articles.