Sufficient Statistic

Jump to: navigation, search

A Sufficient Statistic is a statistical inference based on the sufficiency principle.



  • (Wikipedia, 2015) ⇒ Retrieved 2016-07-30
    • In statistics, a statistic is sufficient with respect to a statistical model and its associated unknown parameter if "no other statistic that can be calculated from the same sample provides any additional information as to the value of the parameter". In particular, a statistic is sufficient for a family of probability distributions if the sample from which it is calculated gives no additional information than does the statistic, as to which of those probability distributions is that of the population from which the sample was taken.

      Roughly, given a set [math] \mathbf{X}[/math] of independent identically distributed data conditioned on an unknown parameter [math]\theta[/math], a sufficient statistic is a function [math]T(\mathbf{X})[/math] whose value contains all the information needed to compute any estimate of the parameter (e.g. a maximum likelihood estimate). Due to the factorization theorem (see below), for a sufficient statistic [math]T(\mathbf{X})[/math], the joint distribution can be written as [math]p(\mathbf{X}) = h(\mathbf{X}) \, g(\theta, T(\mathbf{X}))\,[/math]. From this factorization, it can easily be seen that the maximum likelihood estimate of [math]\theta[/math] will interact with [math]\mathbf{X}[/math] only through [math]T(\mathbf{X})[/math]. Typically, the sufficient statistic is a simple function of the data, e.g. the sum of all the data points.

      More generally, the "unknown parameter" may represent a vector of unknown quantities or may represent everything about the model that is unknown or not fully specified. In such a case, the sufficient statistic may be a set of functions, called a jointly sufficient statistic. Typically, there are as many functions as there are parameters. For example, for a Gaussian distribution with unknown mean and variance, the jointly sufficient statistic, from which maximum likelihood estimates of both parameters can be estimated, consists of two functions, the sum of all data points and the sum of all squared data points (or equivalently, the sample mean and sample variance).

      The concept, due to Ronald Fisher, is equivalent to the statement that, conditional on the value of a sufficient statistic for a parameter, the joint probability distribution of the data does not depend on that parameter. Both the statistic and the underlying parameter can be vectors.