Algebraic statistics

From testwiki
Jump to navigation Jump to search

Template:Short description Algebraic statistics is a branch of mathematical statistics that focuses on the use of algebraic, geometric, and combinatorial methods in statistics. While the use of these methods has a long history in statistics, algebraic statistics is continuously forging new interdisciplinary connections.

This growing field has established itself squarely at the intersection of several areas of mathematics, including, for instance, multilinear algebra, commutative algebra, algebraic geometry, convex geometry, combinatorics, theoretical problems in statistics, and their practical applications. For example, algebraic statistics has been useful for experimental design, parameter estimation, and hypothesis testing.

History

Algebraic statistics can be traced back to Karl Pearson, who used polynomial algebra to study Gaussian mixture models. Subsequently, Ronald A. Fisher, Henry B. Mann, and Rosemary A. Bailey applied Abelian groups to the design of experiments. Experimental designs were also studied with affine geometry over finite fields and then with the introduction of association schemes by R. C. Bose. Orthogonal arrays were introduced by C. R. Rao also for experimental designs.

The field experienced a major revitalization in the 1990s. In 1998, Diaconis and Sturmfels introduced Gröbner bases for constructing Markov chain Monte Carlo algorithms for conditional sampling from discrete exponential families. Pistone and Wynn, in 1996, applied computational commutative algebra to the design and analysis of experiments, providing new tools for understanding confounding and identifiability in complex experimental settings. These works, along with the monograph by Giovanni Pistone, Eva Riccomagno, and Henry P. Wynn, in which the term “algebraic statistics” was first used, played a pivotal role in establishing this field as a unified area of research.

Modern researchers in algebraic statistics explore a wide range of topics, including computational biology, graphical models, and statistical learning.

Active Research Areas

  • Phylogenetics
  • Maximum likelihood estimation
  • Method of moments
  • Graphical models
  • Tropical statistics

Introductory Example

Consider a random variable X which can take on the values 0,1,2. Such a variable is completely characterized by the three probabilities

pi=Pr(X=i),i=0,1,2

and these numbers satisfy

∑i=02pi=1and0≤pi≤1.

Conversely, any three such numbers unambiguously specify a random variable, so we can identify the random variable X with the tuple (p0,p1,p2)∈ℝ3.

Now suppose X is a binomial random variable with parameter q and n=2, i.e. X represents the number of successes when repeating a certain experiment two times, where each experiment has an individual success probability of q. Then

pi=Pr(X=i)=(2i)qi(1−q)2−i

and it is not hard to show that the tuples (p0,p1,p2) which arise in this way are precisely the ones satisfying

4p0p2−p12=0. 

The latter is a polynomial equation defining an algebraic variety (or surface) in ℝ3, and this variety, when intersected with the simplex given by

∑i=02pi=1and0≤pi≤1,

yields a piece of an algebraic curve which may be identified with the set of all 3-state Bernoulli variables. Determining the parameter q amounts to locating one point on this curve; testing the hypothesis that a given variable X is Bernoulli amounts to testing whether a certain point lies on that curve or not.

References

Template:Reflist

  1. ↑ Template:Cite web
  2. ↑ A gap in Garrett Birkhoff's original proof was filled by Alexander Ostrowski. Garrett Birkhoff (1967), Lattice Theory, 3rd ed. Vol. 25 of AMS Colloquium Publications, American Mathematical Society.