Spinors in three dimensions

In probability and statistics, the quantile function specifies, for a given probability in the probability distribution of a random variable, the value at which the probability of the random variable will be less than or equal to that probability. It is also called the percent point function or inverse cumulative distribution function. The quantile function is one way of prescribing a probability distribution, and it is an alternative to the probability density function (pdf) or probability mass function, the cumulative distribution function (cdf) and the characteristic function. The quantile function, Q, of a probability distribution is the inverse of its cumulative distribution function F. The derivative of the quantile function, namely the quantile density function, is yet another way of prescribing a probability distribution. It is the reciprocal of the pdf composed with the quantile function.

Definition

Assuming a continuous and strictly monotonic distribution function, $\scriptstyle F\colon R\to (0,1)$, the quantile function returns the value below which random draws from the given distribution would fall, (p×100) percent of the time. That is, it returns the value of x such that

$F(x)=\Pr(X\leq x)=p.\,$

Quantile function for a general distribution function

If the probability distribution is discrete rather than continuous then there may be gaps between values in the domain of its cdf, while if the cdf is only weakly monotonic there may be "flat spots" in its range. In either case, the quantile function is

$Q(p)\,=\,\inf \left\{x\in R:p\leq F(x)\right\}$

for a probability 0 < p < 1, and the quantile function returns the minimum value of x for which the previous probability statement holds.

Simple example

For example, the quantile function for Exponential(λ) (i.e. intensity λ and expected value 1/λ) is

$Q(p;\lambda )={\frac {-\ln(1-p)}{\lambda }},\!$

for 0 ≤ p < 1. The quartiles are therefore:

first quartile

$\ln(4/3)/\lambda \,$

median

$\ln(2)/\lambda \,$

third quartile

$\ln(4)/\lambda .\,$

Applications

Quantile functions are used in both statistical applications and Monte Carlo methods.

For statistical applications, users need to know key percentage points of a given distribution. For example, they require the median and 25% and 75% quartiles as in the example above or 5%, 95%, 2.5%, 97.5% levels for other applications such as assessing the statistical significance of an observation whose distribution is known; see the quantile entry. Before the popularization of computers, it was not uncommon for books to have appendices with statistical tables sampling the quantile function (see, e.g., [1]). Statistical applications of quantile functions are discussed extensively by Gilchrist.

Monte-Carlo simulations employ quantile functions to produce non-uniform random or pseudorandom numbers for use in diverse types of simulation calculations. A sample from a given distribution may be obtained in principle by applying its quantile function to a sample from a uniform distribution. The demands, for example, of simulation methods in modern computational finance are focusing increasing attention on methods based on quantile functions, as they work well with multivariate techniques based on either copula or quasi-Monte-Carlo methods and Monte Carlo methods in finance.

Calculation

The evaluation of quantile functions often involves numerical methods, as the example of the exponential distribution above is one of the few distributions where a closed-form expression can be found (others include the uniform, the Weibull, the Tukey lambda (which includes the logistic) and the log-logistic). When the cdf itself has a closed-form expression, one can always use a numerical root-finding algorithm such as the bisection method to invert the cdf. Other algorithms to evaluate quantile functions are given in the Numerical Recipes series of books. Algorithms for common distributions are built into many statistical software packages.

Quantile functions may also be characterized as solutions of non-linear ordinary and partial differential equations. The ordinary differential equations for the cases of the normal, Student, beta and gamma distributions have been given and solved.

The normal distribution

The normal distribution is perhaps the most important case. Because the normal distribution is a location-scale family, its quantile function for arbitrary parameters can be derived from a simple transformation of the quantile function of the standard normal distribution, known as the probit function. Unfortunately, this function has no closed-form representation using basic algebraic functions; as a result, approximate representations are usually used. Thorough composite rational and polynomial approximations have been given by Wichura and Acklam (see his web site in External Links). Non-composite rational approximations have been developed by Shaw (see Monte Carlo recycling in External Links).

Ordinary differential equation for the normal quantile

A non-linear ordinary differential equation for the normal quantile, w(p), may be given. It is

${\frac {d^{2}w}{dp^{2}}}=w\left({\frac {dw}{dp}}\right)^{2}$

with the centre (boundary) conditions

$w\left(1/2\right)=0,\,$

$w'\left(1/2\right)={\sqrt {2\pi }}.\,$

This equation may be solved by several methods, including the classical power series approach. From this solutions of arbitrarily high accuracy may be developed (see Steinbrecher and Shaw, 2008).

The Student's t-distribution

Mining Engineer (Excluding Oil ) Truman from Alma, loves to spend time knotting, largest property developers in singapore developers in singapore and stamp collecting. Recently had a family visit to Urnes Stave Church. This has historically been one of the more intractable cases, as the presence of a parameter, ν, the degrees of freedom, makes the use of rational and other approximations awkward. Simple formulas exist when the ν = 1, 2, 4 and the problem may be reduced to the solution of a polynomial when ν is even. In other cases the quantile functions may be developed as power series. The simple cases are as follows:

ν = 1 (Cauchy distribution)

Mining Engineer (Excluding Oil ) Truman from Alma, loves to spend time knotting, largest property developers in singapore developers in singapore and stamp collecting. Recently had a family visit to Urnes Stave Church.

$Q(p)=\tan(\pi (p-1/2))\!$

ν = 2

$Q(p)=2(p-1/2){\sqrt {\frac {2}{\alpha }}}\!$

ν = 4

$Q(p)=\operatorname {sign} (p-1/2)\,2\,{\sqrt {q-1}}\!$

where

$q={\frac {\cos \left({\frac {1}{3}}\arccos \left({\sqrt {\alpha }}\,\right)\right)}{\sqrt {\alpha }}}\!$

and

$\alpha =4p(1-p).\!$

In the above the "sign" function is +1 for positive arguments, -1 for negative arguments and zero at zero. It should not be confused with the trigonometric sine function.

Quantile mixtures

Analogously to the mixtures of densities, distributions can be defined as quantile mixtures

$Q(p)=\sum _{i=1}^{m}a_{i}Q_{i}(p)$,

where $Q_{i}(p)$, $i=1,\ldots ,m$ are quantile functions and $a_{i}$, $i=1,\ldots ,m$ are the model parameters. The parameters $a_{i}$ must be selected so that $Q(p)$ is a quantile function. Two four-parametric quantile mixtures, the normal-polynomial quantile mixture and the Cauchy-polynomial quantile mixture, are presented by Karvanen.

Non-linear differential equations for quantile functions

The non-linear ordinary differential equation given for normal distribution is a special case of that available for any quantile function whose second derivative exists. In general the equation for a quantile, Q(p), may be given. It is

${\frac {d^{2}Q}{dp^{2}}}=H(Q)\left({\frac {dQ}{dp}}\right)^{2}$

augmented by suitable boundary conditions, where

$H(x)=-{\frac {d\log[f(x)]}{dx}}$

and ƒ(x) is the probability density function. The forms of this equation, and its classical analysis by series and asymptotic solutions, for the cases of the normal, Student, gamma and beta distributions has been elucidated by Steinbrecher and Shaw (2008). Such solutions provide accurate benchmarks, and in the case of the Student, suitable series for live Monte Carlo use.

See also

• Inverse transform sampling
• Percent point

References

43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.

• Abernathy, Roger W. and Smith, Robert P. (1993) *"Applying series expansion to the inverse beta distribution to find percentiles of the F-distribution", ACM Trans. Math. Softw., 9 (4), 478–480 21 year-old Glazier James Grippo from Edam, enjoys hang gliding, industrial property developers in singapore developers in singapore and camping. Finds the entire world an motivating place we have spent 4 months at Alejandro de Humboldt National Park.

External links

• An algorithm for computing the inverse normal cumulative distribution function
• Refinement of the Normal Quantile
• New Methods for Managing "Student's" T Distribution
• ACM Algorithm 396: Student's t-Quantiles
• Computational Finance: Differential Equations for Monte Carlo Recycling