Beta distribution
The beta distribution is a bounded continuous distribution. If is often used to express an uncertainty in a proportion, frequency, or percentage, which are all quantities between 0 and 1. With positive «x» and «y» parameters and arbitrary «lower» and «upper» bounds, it is also called a Pert distribution, and in this form is sometimes used as a smooth bell-shaped variation appropriate where a Triangular might otherwise be used. The beta distribution can take on several shapes including bell-shaped unimodal (when a,b>1), bimodal (when 0<a,b<1), uniform (when a=b=1), exponentially decaying (when 0<a<1<b) or exponentially increasing (when 0<b<1<a).
A Beta( n+1, m+1 )
is sometimes used as an estimate for the proportion of individuals with a given trait after observing n individuals with the trait and m individuals without the trait.
Functions
Beta(a, b, lower, upper)
The Beta distribution.
Creates a continuous distribution of numbers between 0 and 1 with a/(a + b) representing the mean, if the optional parameters «lower» and «upper» are omitted. For bounds other than 0 and 1, specify the optional «lower» and «upper» bounds to offset and expand the distribution.
«a» and «b» must be positive.
The probability density of the beta distribution is given by
- [math]\displaystyle{ p(z) = {1\over{BetaFn(a, b)}} z^{a-1} (1-z)^{b-1} }[/math]
where [math]\displaystyle{ 0 \le z \le 1 }[/math].
Dens_Beta(x, a, b, lower, upper)
The probability density at «x».
CumBeta(x, a, b)
The cumulative density, i.e., the probability that the outcome is less than or equal to «x».
CumBetaInv(p, a, b)
The inverse cumulative density, also called the quantile function, which returns the value x with a «p» probability of being larger than the true value.
When to use
Use a beta distribution if the uncertain quantity is bounded by 0 and 1 (or 100%), is continuous, and has a single mode. This distribution is particularly useful for modeling an opinion about the fraction of a population that has some characteristic. For example, if you have observed n members of the population, of which r display the characteristic c, you can represent the uncertainty about the true fraction with c using a beta distribution with parameters «x» = r+1
and «y» = n - r + 1
.
If the uncertain quantity has lower and upper bounds other than 0 and 1, include the «lower» and «upper» bounds parameters to obtain a transformed beta distribution. The transformed beta is a very flexible distribution for representing a wide variety of bounded quantities.
Parameter Estimation
Suppose D
contains sampled historical data indexed by I
, and you want to estimate the «X» and «Y» parameters of the beta distribution from this historical data. With your data in D
normalized to be between the known bounds of 0 and 1, the parameters can be obtained from the following estimation formulas:
X := Var m := Mean(D, I);
Var v := Variance(D, I);
(m^2 - m^3 - v*m) / v
Y := Var m := Mean(D, I);
Var v := Variance(D, I);
(m*(1 - m)^2 - v * (1 - m)) / v
When the range is given and over something other than [0, 1]
, the above estimation formula apply with D
replaced with (D - «lower»)/(«upper» - «lower»)
. Maximum likelihood estimation of all four parameters when «lower» and «upper» are not known is difficult, but worked out in Johnson, Kotz and Balakrishan (1994), Continuous Univariate Distributions, 2ne ed., Volume II, p. 221-235, John Wiley & sons.
See Also
- BetaFn -- the complete beta function
- BetaI -- the incomplete beta function, gives the cumulative density analytically.
- BetaIInv -- the inverse of BetaI.
- Dens_Beta -- the probability density of Beta(a, b) at x
- Beta_m_sd
- Pert -- A parametric variation on the beta distribution
- Variance
- Parametric continuous distributions
- Distribution Densities Library
Enable comment auto-refresher