Pivotal quantity

From HandWiki
Short description: Function of observations and unobservable parameters

In statistics, a pivotal quantity or pivot is a function of observations and unobservable parameters such that the function's probability distribution does not depend on the unknown parameters (including nuisance parameters).[1] A pivot quantity need not be a statistic—the function and its value can depend on the parameters of the model, but its distribution must not. If it is a statistic, then it is known as an ancillary statistic.

More formally,[2] let [math]\displaystyle{ X = (X_1,X_2,\ldots,X_n) }[/math] be a random sample from a distribution that depends on a parameter (or vector of parameters) [math]\displaystyle{ \theta }[/math]. Let [math]\displaystyle{ g(X,\theta) }[/math] be a random variable whose distribution is the same for all [math]\displaystyle{ \theta }[/math]. Then [math]\displaystyle{ g }[/math] is called a pivotal quantity (or simply a pivot).

Pivotal quantities are commonly used for normalization to allow data from different data sets to be compared. It is relatively easy to construct pivots for location and scale parameters: for the former we form differences so that location cancels, for the latter ratios so that scale cancels.

Pivotal quantities are fundamental to the construction of test statistics, as they allow the statistic to not depend on parameters – for example, Student's t-statistic is for a normal distribution with unknown variance (and mean). They also provide one method of constructing confidence intervals, and the use of pivotal quantities improves performance of the bootstrap. In the form of ancillary statistics, they can be used to construct frequentist prediction intervals (predictive confidence intervals).


Normal distribution

One of the simplest pivotal quantities is the z-score; given a normal distribution with mean [math]\displaystyle{ \mu }[/math] and variance [math]\displaystyle{ \sigma^2 }[/math], and an observation x, the z-score:

[math]\displaystyle{ z = \frac{x - \mu}{\sigma}, }[/math]

has distribution [math]\displaystyle{ N(0,1) }[/math] – a normal distribution with mean 0 and variance 1. Similarly, since the n-sample sample mean has sampling distribution [math]\displaystyle{ N(\mu,\sigma^2/n), }[/math] the z-score of the mean

[math]\displaystyle{ z = \frac{\overline{X} - \mu}{\sigma/\sqrt{n}} }[/math]

also has distribution [math]\displaystyle{ N(0,1). }[/math] Note that while these functions depend on the parameters – and thus one can only compute them if the parameters are known (they are not statistics) – the distribution is independent of the parameters.

Given [math]\displaystyle{ n }[/math] independent, identically distributed (i.i.d.) observations [math]\displaystyle{ X = (X_1, X_2, \ldots, X_n) }[/math] from the normal distribution with unknown mean [math]\displaystyle{ \mu }[/math] and variance [math]\displaystyle{ \sigma^2 }[/math], a pivotal quantity can be obtained from the function:

[math]\displaystyle{ g(x,X) = \frac{x - \overline{X}}{s/\sqrt{n}} }[/math]


[math]\displaystyle{ \overline{X} = \frac{1}{n}\sum_{i=1}^n{X_i} }[/math]


[math]\displaystyle{ s^2 = \frac{1}{n-1}\sum_{i=1}^n{(X_i - \overline{X})^2} }[/math]

are unbiased estimates of [math]\displaystyle{ \mu }[/math] and [math]\displaystyle{ \sigma^2 }[/math], respectively. The function [math]\displaystyle{ g(x,X) }[/math] is the Student's t-statistic for a new value [math]\displaystyle{ x }[/math], to be drawn from the same population as the already observed set of values [math]\displaystyle{ X }[/math].

Using [math]\displaystyle{ x=\mu }[/math] the function [math]\displaystyle{ g(\mu,X) }[/math] becomes a pivotal quantity, which is also distributed by the Student's t-distribution with [math]\displaystyle{ \nu = n-1 }[/math] degrees of freedom. As required, even though [math]\displaystyle{ \mu }[/math] appears as an argument to the function [math]\displaystyle{ g }[/math], the distribution of [math]\displaystyle{ g(\mu,X) }[/math] does not depend on the parameters [math]\displaystyle{ \mu }[/math] or [math]\displaystyle{ \sigma }[/math] of the normal probability distribution that governs the observations [math]\displaystyle{ X_1,\ldots,X_n }[/math].

This can be used to compute a prediction interval for the next observation [math]\displaystyle{ X_{n+1}; }[/math] see Prediction interval: Normal distribution.

Bivariate normal distribution

In more complicated cases, it is impossible to construct exact pivots. However, having approximate pivots improves convergence to asymptotic normality.

Suppose a sample of size [math]\displaystyle{ n }[/math] of vectors [math]\displaystyle{ (X_i,Y_i)' }[/math] is taken from a bivariate normal distribution with unknown correlation [math]\displaystyle{ \rho }[/math].

An estimator of [math]\displaystyle{ \rho }[/math] is the sample (Pearson, moment) correlation

[math]\displaystyle{ r = \frac{\frac1{n-1} \sum_{i=1}^n (X_i - \overline{X})(Y_i - \overline{Y})}{s_X s_Y} }[/math]

where [math]\displaystyle{ s_X^2, s_Y^2 }[/math] are sample variances of [math]\displaystyle{ X }[/math] and [math]\displaystyle{ Y }[/math]. The sample statistic [math]\displaystyle{ r }[/math] has an asymptotically normal distribution:

[math]\displaystyle{ \sqrt{n}\frac{r-\rho}{1-\rho^2} \Rightarrow N(0,1) }[/math].

However, a variance-stabilizing transformation

[math]\displaystyle{ z = \rm{tanh}^{-1} r = \frac12 \ln \frac{1+r}{1-r} }[/math]

known as Fisher's z transformation of the correlation coefficient allows creating the distribution of [math]\displaystyle{ z }[/math] asymptotically independent of unknown parameters:

[math]\displaystyle{ \sqrt{n}(z-\zeta) \Rightarrow N(0,1) }[/math]

where [math]\displaystyle{ \zeta = {\rm tanh}^{-1} \rho }[/math] is the corresponding distribution parameter. For finite samples sizes [math]\displaystyle{ n }[/math], the random variable [math]\displaystyle{ z }[/math] will have distribution closer to normal than that of [math]\displaystyle{ r }[/math]. An even closer approximation to the standard normal distribution is obtained by using a better approximation for the exact variance: the usual form is

[math]\displaystyle{ \operatorname{Var}(z) \approx \frac1{n-3} . }[/math]


Main page: Robust statistics

From the point of view of robust statistics, pivotal quantities are robust to changes in the parameters – indeed, independent of the parameters – but not in general robust to changes in the model, such as violations of the assumption of normality. This is fundamental to the robust critique of non-robust statistics, often derived from pivotal quantities: such statistics may be robust within the family, but are not robust outside it.

See also


  1. Shao, J. (2008). "Pivotal quantities". Mathematical Statistics (2nd ed.). New York: Springer. pp. 471–477. ISBN 978-0-387-21718-5. https://books.google.com/books?id=_bEPBwAAQBAJ&pg=PA471. 
  2. DeGroot, Morris H.; Schervish, Mark J. (2011). Probability and Statistics (4th ed.). Pearson. p. 489. ISBN 978-0-321-70970-7.