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 random variable to occur at a given point. The probability for the random variable to fall within a particular region is given by the integral of this variable’s density over the region. The probability density function is nonnegative everywhere, and its integral over the entire space is equal to one.

The terms "probability distribution function"[1] and "probability function"[2] have also sometimes been used to denote the probability density function. However, this use is not standard among probabilists and statisticians. In other sources, "probability distribution function" may be used when the probability distribution is defined as a function over general sets of values, or it may refer to the cumulative distribution function, or it may be a probability mass function rather than the density. Further confusion of terminology exists because density function has also been used for what is here called the "probability mass function".[3]


Absolutely continuous univariate distributions

A probability density function is most commonly associated with absolutely continuous univariate distributions. A random variable X has density f, where f is a non-negative Lebesgue-integrable function, if:

 \operatorname P [a \leq X \leq b] = \int_a^b f(x) \, \mathrm{d}x .

Hence, if F is the cumulative distribution function of X, then:

F(x) = \int_{-\infty}^x f(u) \, \mathrm{d}u ,

and (if f is continuous at x)

 f(x) = \frac{\mathrm{d}}{\mathrm{d}x} F(x) .

Intuitively, one can think of f(x) dx as being the probability of X falling within the infinitesimal interval [xx + dx].

Formal definition

This definition may be extended to any probability distribution using the measure-theoretic definition of probability. A random variable X with values in a measure space (\mathcal{X}, \mathcal{A}) (usually Rn with the Borel sets as measurable subsets) has as probability distribution the measure XP on (\mathcal{X}, \mathcal{A}): the density of X with respect to a reference measure μ on (\mathcal{X}, \mathcal{A}) is the Radon–Nikodym derivative:

f = \frac{\mathrm d X_*P}{\mathrm d \mu} .

That is, f is any measurable function with the property that:

\Pr [X \in A ] = \int_{X^{-1}A} \, \mathrm d P = \int_A f \, \mathrm d \mu

for any measurable set A \in \mathcal{A}.


In the continuous univariate case above, the reference measure is the Lebesgue measure. The probability mass function of a discrete random variable is the density with respect to the counting measure over the sample space (usually the set of integers, or some subset thereof).

Note that it is not possible to define a density with reference to an arbitrary measure (i.e. one can't choose the counting measure as a reference for a continuous random variable). Furthermore, when it does exist, the density is almost everywhere unique.

Further details

Unlike a probability, a probability density function can take on values greater than one; for example, the uniform distribution on the interval [0, ½] has probability density f(x) = 2 for 0 ≤ x ≤ ½ and f(x) = 0 elsewhere.

The standard normal distribution has probability density

    f(x) = \frac{1}{\sqrt{2\pi}}\; e^{-x^2/2}.

If a random variable X is given and its distribution admits a probability density function f, then the expected value of X (if it exists) can be calculated as

    \operatorname{E}[X] = \int_{-\infty}^\infty x\,f(x)\,dx.

Not every probability distribution has a density function: the distributions of discrete random variables do not; nor does the Cantor distribution, even though it has no discrete component, i.e., does not assign positive probability to any individual point.

A distribution has a density function if and only if its cumulative distribution function F(x) is absolutely continuous. In this case: F is almost everywhere differentiable, and its derivative can be used as probability density:

    \frac{d}{dx}F(x) = f(x).

If a probability distribution admits a density, then the probability of every one-point set {a} is zero; the same holds for finite and countable sets.

Two probability densities f and g represent the same probability distribution precisely if they differ only on a set of Lebesgue measure zero.

In the field of statistical physics, a non-formal reformulation of the relation above between the derivative of the cumulative distribution function and the probability density function is generally used as the definition of the probability density function. This alternate definition is the following:

If dt is an infinitely small number, the probability that X is included within the interval (tt + dt) is equal to f(tdt, or:

    \Pr(t<X<t+dt) = f(t)\,dt.

Link between discrete and continuous distributions

It is possible to represent certain discrete random variables as well as random variables involving both a continuous and a discrete part with a generalized probability density function, by using the Dirac delta function. For example, let us consider a binary discrete random variable taking −1 or 1 for values, with probability ½ each.

The density of probability associated with this variable is:

f(t) = \frac{1}{2}(\delta(t+1)+\delta(t-1)).

More generally, if a discrete variable can take n different values among real numbers, then the associated probability density function is:

f(t) = \sum_{i=1}^np_i\, \delta(t-x_i),

where x1, …, xn are the discrete values accessible to the variable and p1, …, pn are the probabilities associated with these values.

This substantially unifies the treatment of discrete and continuous probability distributions. For instance, the above expression allows for determining statistical characteristics of such a discrete variable (such as its mean, its variance and its kurtosis), starting from the formulas given for a continuous distribution.

Densities associated with multiple variables

For continuous random variables X1, …, Xn, it is also possible to define a probability density function associated to the set as a whole, often called joint probability density function. This density function is defined as a function of the n variables, such that, for any domain D in the n-dimensional space of the values of the variables X1, …, Xn, the probability that a realisation of the set variables falls inside the domain D is

\Pr \left( X_1,\ldots,X_N \isin D \right) 
 = \int_D f_{X_1,\dots,X_N}(x_1,\ldots,x_N)\,dx_1 \cdots dx_N.

If F(x1, …, xn) = Pr(X1 ≤ x1, …, Xn ≤ xn) is the cumulative distribution function of the vector (X1, …, Xn), then the joint probability density function can be computed as a partial derivative

    f(x) = \frac{\partial^n F}{\partial x_1 \cdots \partial x_n} \bigg|_x

Marginal densities

For i=1, 2, …,n, let fXi(xi) be the probability density function associated with variable Xi alone. This is called the “marginal” density function, and can be deduced from the probability density associated with the random variables X1, …, Xn by integrating on all values of the n − 1 other variables:

f_{X_i}(x_i) = \int f(x_1,\ldots,x_n)\, dx_1 \cdots dx_{i-1}\,dx_{i+1}\cdots dx_n .


Continuous random variables X1, …, Xn admitting a joint density are all independent from each other if and only if

f_{X_1,\dots,X_n}(x_1,\ldots,x_n) = f_{X_1}(x_1)\cdots f_{X_n}(x_n).


If the joint probability density function of a vector of n random variables can be factored into a product of n functions of one variable

f_{X_1,\dots,X_n}(x_1,\ldots,x_n) = f_1(x_1)\cdots f_n(x_n),

(where each fi is not necessarily a density) then the n variables in the set are all independent from each other, and the marginal probability density function of each of them is given by

f_{X_i}(x_i) = \frac{f_i(x_i)}{\int f_i(x)\,dx}.


This elementary example illustrates the above definition of multidimensional probability density functions in the simple case of a function of a set of two variables. Let us call \vec R a 2-dimensional random vector of coordinates (X, Y): the probability to obtain \vec R in the quarter plane of positive x and y is

\Pr \left( X > 0, Y > 0 \right)
 = \int_0^\infty \int_0^\infty f_{X,Y}(x,y)\,dx\,dy.

Sums of independent random variables

The probability density function of the sum of two independent random variables U and V, each of which has a probability density function, is the convolution of their separate density functions:

f_{U+V}(x) = \int_{-\infty}^\infty f_U(y) f_V(x - y)\,dy
= \left( f_{U} * f_{V} \right) (x)

It is possible to generalize the previous relation to a sum of N independent random variables, with densities U1, …, UN:

f_{U_{1} + \ldots + U_{N}}(x)
= \left( f_{U_{1}} * \ldots * f_{U_{N}} \right) (x)

Dependent variables and change of variables

If the probability density function of a random variable X is given as fX(x), it is possible (but often not necessary; see below) to calculate the probability density function of some variable Y = g(X). This is also called a “change of variable” and is in practice used to generate a random variable of arbitrary shape fg(X) = fY using a known (for instance uniform) random number generator.

If the function g is monotonic, then the resulting density function is

f_Y(y) = \left| \frac{1}{g'(g^{-1}(y))} \right| \cdot f_X(g^{-1}(y)).

Here g−1 denotes the inverse function and g' denotes the derivative.

This follows from the fact that the probability contained in a differential area must be invariant under change of variables. That is,

\left| f_Y(y)\, dy\right| = \left| f_X(x)\, dx\right|,


f_Y(y) = \left| \frac{dx}{dy} \right| f_X(x) = \left| \frac{1}{g'(x)} \right| f_X(x) = \left| \frac{1}{g'(g^{-1}(y))} \right|f_X(g^{-1}(y)).

For functions which are not monotonic the probability density function for y is

\sum_{k=1}^{n(y)} \left| \frac{1}{g'(g^{-1}_{k}(y))} \right| \cdot f_X(g^{-1}_{k}(y))

where n(y) is the number of solutions in x for the equation g(x) = y, and g−1k(y) are these solutions.

It is tempting to think that in order to find the expected value E(g(X)) one must first find the probability density fg(X) of the new random variable Y = g(X). However, rather than computing

 E(g(X)) = \int_{-\infty}^\infty y f_{g(X)}(y)\,dy,

one may find instead

E(g(X)) = \int_{-\infty}^\infty g(x) f_X(x)\,dx.

The values of the two integrals are the same in all cases in which both X and g(X) actually have probability density functions. It is not necessary that g be a one-to-one function. In some cases the latter integral is computed much more easily than the former.

Multiple variables

The above formulas can be generalized to variables (which we will again call y) depending on more than one other variable. f(x0, x1, …, xm−1) shall denote the probability density function of the variables y depends on, and the dependence shall be y = g(x0, x1, …, xm−1). Then, the resulting density function is

 \int\limits_{y = g(x_0, x_1, \dots, x_{m-1})} \frac{f(x_0, x_1,\dots, x_{m-1})}{\sqrt{\sum_{j=0}^{j<m}} (\frac{\partial g}{\partial x_j}(x_0, x_1, \dots , x_{m-1}))^2} \; dV

where the integral is over the entire (m-1)-dimensional solution of the subscripted equation and the symbolic dV must be replaced by a parametrization of this solution for a particular calculation; the variables x0, x1, …, xm−1 are then of course functions of this parametrization.

This derives from the following, perhaps more intuitive representation: Suppose x is an n-dimensional random variable with joint density f. If y = H(x), where H is a bijective, differentiable function, then y has density g:

g(\mathbf{y}) = f(\mathbf{x})\left\vert \det\left(\frac{d\mathbf{x}}{d\mathbf{y}}\right)\right \vert

with the differential regarded as the Jacobian of the inverse of H, evaluated at y.

Using the delta-function (and assuming independence) the same result is formulated as follows.

If the probability density function of independent random variables Xi, i = 1, 2, …n are given as fXi(xi), it is possible to calculate the probability density function of some variable Y = G(X1, X2, …Xn). The following formula establishes a connection between the probability density function of Y denoted by fY(y) and fXi(xi) using the Dirac delta function:

f_Y(y) = \int_{-\infty}^\infty \int_{-\infty}^\infty \ldots \int_{-\infty}^\infty f_{X_1}(x_1)f_{X_2}(x_2) \ldots f_{X_n}(x_n)\delta(y-G(x_1,x_2,\ldots x_n))\,dx_1\,dx_2\,\ldots dx_n

See also


  1. ^ Probability distribution function PlanetMath
  2. ^ Probability Function at Mathworld
  3. ^ Ord, J.K. (1972) Families of Frequency Distributions, Griffin. ISBN 0-85264-137-0 (for example, Table 5.1 and Example 5.4)


  • Pierre Simon de Laplace (1812). Analytical Theory of Probability. 
The first major treatise blending calculus with probability theory, originally in French: Théorie Analytique des Probabilités.
The modern measure-theoretic foundation of probability theory; the original German version (Grundbegriffe der Wahrscheinlichkeitrechnung) appeared in 1933.
  • Patrick Billingsley (1979). Probability and Measure. New York, Toronto, London: John Wiley and Sons. ISBN 0471007102. 
  • David Stirzaker (2003). Elementary Probability. ISBN 0521420288. 
Chapters 7 to 9 are about continuous variables. This book is filled with theory and mathematical proofs.

External links

Wikimedia Foundation. 2010.

Look at other dictionaries:

  • probability density function — tikimybės tankio funkcija statusas T sritis automatika atitikmenys: angl. probability density function vok. Wahrscheinlichkeitsdichtefunktion, f rus. функция плотности вероятности, f pranc. fonction de densité de probabilité, f …   Automatikos terminų žodynas

  • probability density function — Statistics. 1. a function of a continuous variable whose integral over a region gives the probability that a random variable falls within the region. 2. Also called frequency function. a function of a discrete variable whose sum over a discrete… …   Universalium

  • probability density function — noun Date: 1957 1. probability function 2. a function of a continuous random variable whose integral over an interval gives the probability that its value will fall within the interval …   New Collegiate Dictionary

  • probability density function — probabil′ity den sity func tion n. 1) sta a function of a continuous variable whose integral over a region gives the probability that a random variable falls within the region 2) sta a function of a discrete variable whose sum over a discrete set …   From formal English to slang

  • probability density function — noun 1. : probability function 2. : a function of a continuous random variable whose integral over an interval gives the probability that its value will fall within the interval …   Useful english dictionary

  • probability density function — in statistics, a mathematical function that describes the distribution of measurements for a population; a curve that describes a population. Its integral is the cumulative distribution function, so the probability that an individual measurement… …   Medical dictionary

  • probability density function — noun Any function whose integral over a set gives the probability that a random variable has a value in that set …   Wiktionary

  • Probability density function — The probability function for a continuous random variable. The New York Times Financial Glossary …   Financial and business terms

  • probability density function — The function that describes the change of certain realizations for a continuous random variable. Bloomberg Financial Dictionary …   Financial and business terms

  • Probability Density Function - PDF — A statistical measure that defines a probability distribution for a random variable and is often denoted as f(x). When the PDF function is graphically portrayed, the area under the graph will indicate the interval under which the variable with… …   Investment dictionary

Share the article and excerpts

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

We are using cookies for the best presentation of our site. Continuing to use this site, you agree with this.