Difference between revisions of "Probability density function"
(→Example: Quotient of two standard normals: I removed on instance of repeated identical equation lines) 
en>Omnipaedista (style fix) 

Line 1:  Line 1:  
[[Image:Boxplot vs PDF.svgthumb350px[[Boxplot]] and probability density function of a [[normal distribution]] {{nowrap''N''(0, ''σ''<sup>2</sup>)}}.]]  [[Image:Boxplot vs PDF.svgthumb350px[[Boxplot]] and probability density function of a [[normal distribution]] {{nowrap''N''(0, ''σ''<sup>2</sup>)}}.]]  
−  In [[probability theory]], a '''probability density function''' ('''  +  In [[probability theory]], a '''probability density function''' ('''PDF'''), or '''density''' of a [[continuous random variable]], is a [[Function (mathematics)function]] that describes the relative likelihood for this random variable to take on a given value. The probability of the random variable falling within a particular range of values is given by the [[integral]] of this variable’s density over that range—that is, it is given by the area under the density function but above the horizontal axis and between the lowest and greatest values of the range. The probability density function is nonnegative everywhere, and its integral over the entire space is equal to one. 
The terms "''probability distribution function''"<ref>[http://planetmath.org/?method=png&from=objects&id=2884&op=getobj Probability distribution function] PlanetMath</ref> and "''probability function''"<ref>[http://mathworld.wolfram.com/ProbabilityFunction.html Probability Function] at Mathworld </ref> 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".<ref>Ord, J.K. (1972) ''Families of Frequency Distributions'', Griffin. ISBN 0852641370 (for example, Table 5.1 and Example 5.4)</ref>  The terms "''probability distribution function''"<ref>[http://planetmath.org/?method=png&from=objects&id=2884&op=getobj Probability distribution function] PlanetMath</ref> and "''probability function''"<ref>[http://mathworld.wolfram.com/ProbabilityFunction.html Probability Function] at Mathworld </ref> 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".<ref>Ord, J.K. (1972) ''Families of Frequency Distributions'', Griffin. ISBN 0852641370 (for example, Table 5.1 and Example 5.4)</ref>  
+  
+  ==Example==  
+  Suppose a species of bacteria typically lives 4 to 6 hours. What is the probability that a bacterium lives ''exactly'' 5 hours? The answer is 0%! Lots of bacteria live for ''approximately'' 5 hours, but there is negligible chance that any given bacterium dies at ''exactly'' 5.0000000000... hours.  
+  
+  Instead we might ask: What is the probability that the bacterium dies between 5 hours and 5.01 hours? Let's say the answer is "2%". Next: What is the probability that the bacterium lives between 5 hours and 5.001 hours? The answer is probably around 0.2%, since this is 1/10th of the previous interval. The probability that the bacterium lives between 5 hours and 5.0001 hours is about 0.02%, and so on.  
+  
+  Therefore, in response to the question "What is the probability that the bacterium lives 5 hours?", a literally correct but unhelpful answer is "0", but a better answer is (2 hour<sup>–1</sup>) ''dt''. This is the probability that the bacterium dies within a small ([[infinitesimal]]) window of time around 5 hours, where ''dt'' is the duration of this window. For example, the probability that it lives longer than 5 hours, but shorter than (5 hours + 1 nanosecond), is (2 hour<sup>–1</sup>)×(1 nanosecond) = 6×10<sup>13</sup>.  
+  
+  The quantity 2 hour<sup>–1</sup> is called the ''probability density'' for the bacterium to live 5 hours. There is a ''probability density function'' ''f'' with ''f''(5 hours) = 2 hour<sup>–1</sup>. The [[integral]] of ''f'' over any window of time (not only infinitesimal windows but also large windows) is the probability that the bacterium dies in that window.  
==Absolutely continuous univariate distributions==  ==Absolutely continuous univariate distributions==  
Line 150:  Line 159:  
or  or  
−  : <math>f_Y(y) = \left \frac{dx}{dy} \right f_X(x) = \left \frac{d}{dy} (x) \right f_X(x) = \left \frac{d}{dy} (g^{1}(y)) \rightf_X(g^{1}(y)).</math>  +  : <math>f_Y(y) = \left \frac{dx}{dy} \right f_X(x) = \left \frac{d}{dy} (x) \right f_X(x) = \left \frac{d}{dy} (g^{1}(y)) \rightf_X(g^{1}(y)) = \frac{f_X(g^{1}(y))}{g'(g^{1}(y))}.</math> 
For functions which are not monotonic the probability density function for ''y'' is  For functions which are not monotonic the probability density function for ''y'' is  
Line 178:  Line 187:  
This derives from the following, perhaps more intuitive representation: Suppose '''''x''''' is an ''n''dimensional random variable with joint density ''f''. If {{nowrap'''''y''''' {{=}} ''H''('''''x''''')}}, where ''H'' is a [[bijective]], [[differentiable]] function, then '''''y''''' has density ''g'':  This derives from the following, perhaps more intuitive representation: Suppose '''''x''''' is an ''n''dimensional random variable with joint density ''f''. If {{nowrap'''''y''''' {{=}} ''H''('''''x''''')}}, where ''H'' is a [[bijective]], [[differentiable]] function, then '''''y''''' has density ''g'':  
−  : <math>g(\mathbf{y}) = f(\mathbf{x})\left\vert \det\left(\frac{d\mathbf{x}}{d\mathbf{y}}\right)\right \vert</math>  +  : <math>g(\mathbf{y}) = f(\mathbf{x})\left\vert \det\left(\frac{\mathrm{d}\mathbf{x}}{\mathrm{d}\mathbf{y}}\right)\right \vert</math> 
with the differential regarded as the [[Jacobian matrix and determinantJacobian]] of the inverse of ''H'', evaluated at '''''y'''''.  with the differential regarded as the [[Jacobian matrix and determinantJacobian]] of the inverse of ''H'', evaluated at '''''y'''''.  
Line 256:  Line 265:  
Given two [[standard normal distributionstandard normal]] variables ''U'' and ''V'', the quotient can be computed as follows. First, the variables have the following density functions:  Given two [[standard normal distributionstandard normal]] variables ''U'' and ''V'', the quotient can be computed as follows. First, the variables have the following density functions:  
−  :<math>p(U) = \frac{1}{\sqrt{2\pi}} e^{U^2  +  :<math>p(U) = \frac{1}{\sqrt{2\pi}} e^{\frac{U^2}{2}}</math> 
−  :<math>p(V) = \frac{1}{\sqrt{2\pi}} e^{V^2  +  :<math>p(V) = \frac{1}{\sqrt{2\pi}} e^{\frac{V^2}{2}}</math> 
We transform as described above:  We transform as described above:  
+  
:<math>Y=U/V</math>  :<math>Y=U/V</math>  
:<math>Z=V</math>  :<math>Z=V</math>  
This leads to:  This leads to:  
−  :<math>  +  
−  \begin{align}  +  :<math>\begin{align} 
p(Y) &= \int_{\infty}^{\infty} p_U(YZ)\,p_V(Z)\, Z \, dZ \\  p(Y) &= \int_{\infty}^{\infty} p_U(YZ)\,p_V(Z)\, Z \, dZ \\  
−  &= \int_{\infty}^{\infty} \frac{1}{\sqrt{2\pi}} e^{Y^2Z^  +  &= \int_{\infty}^{\infty} \frac{1}{\sqrt{2\pi}} e^{\frac{1}{2}Y^2Z^2} \frac{1}{\sqrt{2\pi}} e^{\frac{1}{2}Z^2} Z \, dZ \\ 
−  &= \int_{\infty}^{\infty} \frac{1}{2\pi} e^{(Y^2+1)Z^  +  &= \int_{\infty}^{\infty} \frac{1}{2\pi} e^{\frac{1}{2}(Y^2+1)Z^2} Z \, dZ \\ 
−  &= 2\int_{0}^{\infty} \frac{1}{2\pi} e^{(Y^2+1)Z^  +  &= 2\int_{0}^{\infty} \frac{1}{2\pi} e^{\frac{1}{2}(Y^2+1)Z^2} Z \, dZ \\ 
−  &= \int_{0}^{\infty} \frac{1}{\pi} e^{(Y^2+1)u} \, du \  +  &= \int_{0}^{\infty} \frac{1}{\pi} e^{(Y^2+1)u} \, du && u=\tfrac{1}{2}Z^2\\ 
−  &= \frac{1}{\pi(Y^2+1)} e^{(Y^2+1)u}\  +  &= \left. \frac{1}{\pi(Y^2+1)} e^{(Y^2+1)u}\right]_{u=0}^{\infty} \\ 
&= \frac{1}{\pi(Y^2+1)}  &= \frac{1}{\pi(Y^2+1)}  
−  \end{align}  +  \end{align}</math> 
−  </math>  
This is a standard [[Cauchy distribution]].  This is a standard [[Cauchy distribution]].  
==See also==  ==See also==  
+  {  
+    
+  valign="top"  
* [[Density estimation]]  * [[Density estimation]]  
* [[Likelihood function]]  * [[Likelihood function]]  
Line 284:  Line 296:  
* [[Probability mass function]]  * [[Probability mass function]]  
* [[Secondary measure]]  * [[Secondary measure]]  
+    
+   
+  valign="top"  
+  * Uses as ''position probability density'':  
+  ** [[Atomic orbital]]  
+  ** [[Home range]]  
+  }  
==Bibliography==  ==Bibliography== 
Latest revision as of 08:44, 29 December 2014
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 take on a given value. The probability of the random variable falling within a particular range of values is given by the integral of this variable’s density over that range—that is, it is given by the area under the density function but above the horizontal axis and between the lowest and greatest values of the range. 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]}
Contents
 1 Example
 2 Absolutely continuous univariate distributions
 3 Formal definition
 4 Further details
 5 Link between discrete and continuous distributions
 6 Families of densities
 7 Densities associated with multiple variables
 8 Dependent variables and change of variables
 9 Sums of independent random variables
 10 Products and quotients of independent random variables
 11 See also
 12 Bibliography
 13 External links
Example
Suppose a species of bacteria typically lives 4 to 6 hours. What is the probability that a bacterium lives exactly 5 hours? The answer is 0%! Lots of bacteria live for approximately 5 hours, but there is negligible chance that any given bacterium dies at exactly 5.0000000000... hours.
Instead we might ask: What is the probability that the bacterium dies between 5 hours and 5.01 hours? Let's say the answer is "2%". Next: What is the probability that the bacterium lives between 5 hours and 5.001 hours? The answer is probably around 0.2%, since this is 1/10th of the previous interval. The probability that the bacterium lives between 5 hours and 5.0001 hours is about 0.02%, and so on.
Therefore, in response to the question "What is the probability that the bacterium lives 5 hours?", a literally correct but unhelpful answer is "0", but a better answer is (2 hour^{–1}) dt. This is the probability that the bacterium dies within a small (infinitesimal) window of time around 5 hours, where dt is the duration of this window. For example, the probability that it lives longer than 5 hours, but shorter than (5 hours + 1 nanosecond), is (2 hour^{–1})×(1 nanosecond) = 6×10^{13}.
The quantity 2 hour^{–1} is called the probability density for the bacterium to live 5 hours. There is a probability density function f with f(5 hours) = 2 hour^{–1}. The integral of f over any window of time (not only infinitesimal windows but also large windows) is the probability that the bacterium dies in that window.
Absolutely continuous univariate distributions
A probability density function is most commonly associated with absolutely continuous univariate distributions. A random variable X has density f_{X}, where f_{X} is a nonnegative Lebesgueintegrable function, if:
Hence, if F_{X} is the cumulative distribution function of X, then:
and (if f_{X} is continuous at x)
Intuitively, one can think of f_{X}(x) dx as being the probability of X falling within the infinitesimal interval [x, x + dx].
Formal definition
(This definition may be extended to any probability distribution using the measuretheoretic definition of probability.)
A random variable X with values in a measurable space (usually R^{n} with the Borel sets as measurable subsets) has as probability distribution the measure X_{∗}P on : the density of X with respect to a reference measure μ on is the Radon–Nikodym derivative:
That is, f is any measurable function with the property that:
Discussion
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 (e.g. 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
If a random variable X is given and its distribution admits a probability density function f, then the expected value of X (if the expected value exists) can be calculated as
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:
If a probability distribution admits a density, then the probability of every onepoint 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 nonformal 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 (t, t + dt) is equal to f(t) dt, or:
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 having the Rademacher distribution—that is, taking −1 or 1 for values, with probability ½ each. The density of probability associated with this variable is:
More generally, if a discrete variable can take n different values among real numbers, then the associated probability density function is:
where x_{1}, …, x_{n} are the discrete values accessible to the variable and p_{1}, …, p_{n} 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 of the probability.
Families of densities
It is common for probability density functions (and probability mass functions) to be parametrized—that is, to be characterized by unspecified parameters. For example, the normal distribution is parametrized in terms of the mean and the variance, denoted by and respectively, giving the family of densities
It is important to keep in mind the difference between the domain of a family of densities and the parameters of the family. Different values of the parameters describe different distributions of different random variables on the same sample space (the same set of all possible values of the variable); this sample space is the domain of the family of random variables that this family of distributions describes. A given set of parameters describes a single distribution within the family sharing the functional form of the density. From the perspective of a given distribution, the parameters are constants, and terms in a density function that contain only parameters, but not variables, are part of the normalization factor of a distribution (the multiplicative factor that ensures that the area under the density—the probability of something in the domain occurring— equals 1). This normalization factor is outside the kernel of the distribution.
Since the parameters are constants, reparametrizing a density in terms of different parameters, to give a characterization of a different random variable in the family, means simply substituting the new parameter values into the formula in place of the old ones. Changing the domain of a probability density, however, is trickier and requires more work: see the section below on change of variables.
Densities associated with multiple variables
For continuous random variables X_{1}, …, X_{n}, 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 ndimensional space of the values of the variables X_{1}, …, X_{n}, the probability that a realisation of the set variables falls inside the domain D is
If F(x_{1}, …, x_{n}) = Pr(X_{1} ≤ x_{1}, …, X_{n} ≤ x_{n}) is the cumulative distribution function of the vector (X_{1}, …, X_{n}), then the joint probability density function can be computed as a partial derivative
Marginal densities
For i=1, 2, …,n, let f_{Xi}(x_{i}) be the probability density function associated with variable X_{i} alone. This is called the “marginal” density function, and can be deduced from the probability density associated with the random variables X_{1}, …, X_{n} by integrating on all values of the n − 1 other variables:
Independence
Continuous random variables X_{1}, …, X_{n} admitting a joint density are all independent from each other if and only if
Corollary
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
(where each f_{i} 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
Example
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 a 2dimensional random vector of coordinates (X, Y): the probability to obtain in the quarter plane of positive x and y is
Dependent variables and change of variables
If the probability density function of a random variable X is given as f_{X}(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 f_{g(X)} = f_{Y} using a known (for instance uniform) random number generator.
If the function g is monotonic, then the resulting density function is
Here g^{−1} denotes the inverse function.
This follows from the fact that the probability contained in a differential area must be invariant under change of variables. That is,
or
For functions which are not monotonic the probability density function for y is
where n(y) is the number of solutions in x for the equation g(x) = y, and g^{−1}_{k}(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 f_{g(X)} of the new random variable Y = g(X). However, rather than computing
one may find instead
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 onetoone 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(x_{1}, …, x_{n}) shall denote the probability density function of the variables that y depends on, and the dependence shall be y = g(x_{1}, …, x_{n}). Then, the resulting density function is{{ safesubst:#invoke:Unsubstdate=__DATE__ $B= {{#invoke:Category handlermain}}{{#invoke:Category handlermain}}^{[citation needed]} }}
where the integral is over the entire (n1)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 x_{1}, …, x_{n} are then of course functions of this parametrization.
This derives from the following, perhaps more intuitive representation: Suppose x is an ndimensional random variable with joint density f. If y = H(x), where H is a bijective, differentiable function, then y has density g:
with the differential regarded as the Jacobian of the inverse of H, evaluated at y.
Using the deltafunction (and assuming independence) the same result is formulated as follows.
If the probability density function of independent random variables X_{i}, i = 1, 2, …n are given as f_{Xi}(x_{i}), it is possible to calculate the probability density function of some variable Y = G(X_{1}, X_{2}, …X_{n}). The following formula establishes a connection between the probability density function of Y denoted by f_{Y}(y) and f_{Xi}(x_{i}) using the Dirac delta function:
Sums of independent random variables
{{#invoke:see alsoseealso}}
 Not to be confused with Mixture distribution
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:
It is possible to generalize the previous relation to a sum of N independent random variables, with densities U_{1}, …, U_{N}:
This can be derived from a twoway change of variables involving Y=U+V and Z=V, similarly to the example below for the quotient of independent random variables.
Products and quotients of independent random variables
{{#invoke:see alsoseealso}}
Given two independent random variables U and V, each of which has a probability density function, the density of the product Y=UV and quotient Y=U/V can be computed by a change of variables.
Example: Quotient distribution
To compute the quotient Y=U/V of two independent random variables U and V, define the following transformation:
Then, the joint density p(Y,Z) can be computed by a change of variables from U,V to Y,Z, and Y can be derived by marginalizing out Z from the joint density.
The inverse transformation is
The Jacobian matrix of this transformation is
Thus:
And the distribution of Y can be computed by marginalizing out Z:
Note that this method crucially requires that the transformation from U,V to Y,Z be bijective. The above transformation meets this because Z can be mapped directly back to V, and for a given V the quotient U/V is monotonic. This is similarly the case for the sum U+V, difference UV and product UV.
Exactly the same method can be used to compute the distribution of other functions of multiple independent random variables.
Example: Quotient of two standard normals
Given two standard normal variables U and V, the quotient can be computed as follows. First, the variables have the following density functions:
We transform as described above:
This leads to:
This is a standard Cauchy distribution.
See also


Bibliography
 {{#invoke:citation/CS1citation
CitationClass=book }}
 The first major treatise blending calculus with probability theory, originally in French: Théorie Analytique des Probabilités.
 {{#invoke:citation/CS1citation
CitationClass=book }}
 The modern measuretheoretic foundation of probability theory; the original German version (Grundbegriffe der Wahrscheinlichkeitsrechnung) appeared in 1933.
 {{#invoke:citation/CS1citation
CitationClass=book }}
 {{#invoke:citation/CS1citation
CitationClass=book }}
 Chapters 7 to 9 are about continuous variables.
 ↑ Probability distribution function PlanetMath
 ↑ Probability Function at Mathworld
 ↑ Ord, J.K. (1972) Families of Frequency Distributions, Griffin. ISBN 0852641370 (for example, Table 5.1 and Example 5.4)
External links
 {{#invoke:citation/CS1citation
CitationClass=citation }}