Quantile function

Quantile function

: "See also quantile."

In probability theory, a quantile function of a
probability distribution is the inverse "F" −1 of its cumulative distribution function (cdf) "F". Assuming a continuous and strictly monotonic distribution function, scriptstyle Fcolon R o (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

:Pr(X le x) = p.,

If the probability distribution is discrete rather than continuous thenthere 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

:F^{-1}(p) = infleft{ xin R : p le F(x) ight}

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(&lambda;) is

:F^{-1}(p;lambda) = frac{-ln(1-p)}{lambda}, !

for 0 &le; "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. Statistical applications of quantile functions are discussed extensively by Gilchrist (2000).

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 (see Jackel, 2002) 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, Weibull, logistic and 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 in to 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 (see Steinbrecher and Shaw, 2008).

The normal distribution

The normal distribution is perhaps the most important case, and, in the absence of a simple formula, approximate representations are usually used. Thorough composite rational and polynomial approximations have been given by Wichura (1988) and Acklam (see his web site in External Links). Also see the entry on the probit function.

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}{d p^2} = w left(frac{d w}{d p} ight)^2

with the centre (boundary) conditions

:wleft(1/2 ight) = 0,,

:w'left(1/2 ight) = sqrt{2pi}.,

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

This has historically been one of the more intractable cases, as the presence of a parameter, &nu;, the degrees of freedom, makes the use of rational and other approximations awkward. Simple formulas exist when the &nu; = 1, 2, 4 and the problem may be reduced to the solution of a polynomial when &nu; is even. In other cases the quantile functions may be developed as power series (see Shaw (2006) for details). The simple cases are as follows:


= &nu; = 1 (Cauchy distribution)=

:F^{-1}(p) = an (pi(p-1/2)) !


= &nu; = 2=

:F^{-1}(p) = frac{2p-1}{sqrt{2p(1-p) !


= &nu; = 4 =

:F^{-1}(p) = operatorname{sign}(p-1/2)sqrt{q-4}!

where

:q = frac{4}{sqrt{alpha cos left( frac{1}{3} arccos left( sqrt{alpha} , ight) ight)!

and

:alpha = 4p(1-p).!

See also

* Inverse transform sampling

References

*cite book
author=Gilchrist, W.
year=2000
title=Statistical Modelling with Quantile Functions

*cite book
author=Jaeckel, P.
year=2002
title=Monte Carlo methods in finance

*cite journal
author=Wichura, M.J.
year=1988
title=Algorithm AS241: The Percentage Points of the Normal Distribution
journal=Applied Statistics
volume=37
pages=477–484
doi=10.2307/2347330

*cite journal
author=Shaw, W.T.
year=2006
title=Sampling Student’s T distribution – use of the inverse cumulative distribution function.
journal=Journal of Computational Finance
volume=9
issue=4
pages=37–73

*cite journal
author=Steinbrecher, G., Shaw, W.T.
year=2008
title=Quantile mechanics
journal=European Journal of Applied Mathematics
volume=19
issue=2
pages=87–112
doi=10.1017/S0956792508007341

External links

* [http://home.online.no/~pjacklam/notes/invnorm/] An algorithm for computing the inverse normal cumulative distribution function.

* [http://www.mth.kcl.ac.uk/~shaww/web_page/papers/NormalQuantile1.pdf] Refinement of the Normal Quantile

* [http://www.mth.kcl.ac.uk/~shaww/web_page/papers/Tdistribution06.pdf] New Method's for Managing "Student's" T Distribution

* [http://portal.acm.org/citation.cfm?id=355600] ACM Algorithm 396: Student's t-Quantiles

* [http://portal.acm.org/citation.cfm?id=168387] Applying series expansion to the inverse beta distribution to find percentiles of the F-distribution


Wikimedia Foundation. 2010.

Игры ⚽ Поможем сделать НИР

Look at other dictionaries:

  • Quantile — Quantiles are points taken at regular intervals from the cumulative distribution function of a random variable. Dividing ordered data into q essentially equal sized data subsets is the motivation for q quantiles; the quantiles are the data values …   Wikipedia

  • quantile — The argument of the cumulative distribution function corresponding to a specified probability; (of a sample) the value below which occur a specified proportion of the observations in the ordered set of observations …   Dictionary of automotive terms

  • Cumulative distribution function — for the normal distributions in the image below …   Wikipedia

  • Error function — Plot of the error function In mathematics, the error function (also called the Gauss error function) is a special function (non elementary) of sigmoid shape which occurs in probability, statistics and partial differential equations. It is defined …   Wikipedia

  • Characteristic function (probability theory) — The characteristic function of a uniform U(–1,1) random variable. This function is real valued because it corresponds to a random variable that is symmetric around the origin; however in general case characteristic functions may be complex valued …   Wikipedia

  • Probability density function — Boxplot and probability density function of a normal distribution N(0, σ2). In probability theory, a probability density function (pdf), or density of a continuous random variable is a function that describes the relative likelihood for this… …   Wikipedia

  • Moment-generating function — In probability theory and statistics, the moment generating function of any random variable is an alternative definition of its probability distribution. Thus, it provides the basis of an alternative route to analytical results compared with… …   Wikipedia

  • Q-Q plot — Not to be confused with P P plot. A normal Q Q plot of randomly generated, independent standard exponential data, (X   Exp(1)). This Q Q plot compares a sample of data on the vertical axis to a statistical population on the horizontal… …   Wikipedia

  • Normal distribution — This article is about the univariate normal distribution. For normally distributed vectors, see Multivariate normal distribution. Probability density function The red line is the standard normal distribution Cumulative distribution function …   Wikipedia

  • Cauchy distribution — Not to be confused with Lorenz curve. Cauchy–Lorentz Probability density function The purple curve is the standard Cauchy distribution Cumulative distribution function …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”