Beta function explained

In mathematics, the beta function, also called the Euler integral of the first kind, is a special function that is closely related to the gamma function and to binomial coefficients. It is defined by the integral

\Beta(z _{1,z} _{2)} = \int ^{1} _{0} t ^{z _{1-1}} ( 1- t) ^{z _{2-1}} dt for complex number inputs z _{1,} z _{2} such that \Re(z _{1),} \Re(z _{2)>0} .

The beta function was studied by Leonhard Euler and Adrien-Marie Legendre and was given its name by Jacques Binet; its symbol is a Greek capital beta.

Properties

The beta function is symmetric, meaning that

\Beta(z _{1,z} _{2)} = \Beta(z _{2,z} _{1)} for all inputs z _{1} and z _{2} .^[1]

A key property of the beta function is its close relationship to the gamma function:^[1]

\Beta(z _{1,z} ^{2 ) = (\Gamma(z _{1)\Gamma(z} _{2)})/(\Gamma(z _{1+z} _{2)})}

A proof is given below in .

The beta function is also closely related to binomial coefficients. When (or, by symmetry) is a positive integer, it follows from the definition of the gamma function that^[1]

\Beta(m,n) = ((m -1 )!(n -1 )!)/((m + n -1 )!) = (m + n)/(mn) / \binom{m + n}{m}

Relationship to the gamma function

A simple derivation of the relation

\Beta(z _{1,z} _{2)} = (\Gamma(z _{1)\Gamma(z} _{2)})/(\Gamma(z _{1+z} _{2)}) can be found in Emil Artin's book The Gamma Function, page 18–19.To derive this relation, write the product of two factorials as integrals. Since they are integrals in two separate variables, we can combine then into an iterated integral:\begin{align} \Gamma(z _{1)\Gamma(z} _{2)} & = \int ^{inf ty e} _{u =0} ^{-u} u ^{z _{1-1}} du ⋅ \int ^{inf ty e} _{v =0} ^{-v} v ^{z _{2-1}} dv \\[ 6 pt] & = \int ^{inf ty e} _{u =0} ^{-u-v} u ^{z _{1-1}} v ^{z _{2-1}} du dv. \end{align}

Changing variables by and, because and, we have that the limits of integrations for are 0 to ∞ and the limits of integration for are 0 to 1. Thus produces

\begin{align} \Gamma(z _{1)\Gamma(z} _{2)} & = \int ^{1} _{t =0} e ^{-s} (st) ^{z _{1-1}} (s( 1- t)) ^{z _{2-1}} sdt ds \\[ 6 pt] & = \int ^{inf ty} _{s =0} e ^{-s} s ^{z _{1+z} _{2-1}} ds ⋅ \int ^{1} _{t =0} t ^{z _{1-1}} ( 1- t) ^{z _{2-1}} dt\\ & = \Gamma(z _{1+z} _{2)} ⋅ \Beta(z _{1,z} _{2). \end{align}}

Dividing both sides by

\Gamma(z _{1+z} _{2)} gives the desired result.

The stated identity may be seen as a particular case of the identity for the integral of a convolution. Taking

\begin{align}f(u)&: = e ^{-u} u ^{z _{1-1}} 1 _{\R _{+}} \ g(u)&: = e ^{-u} u ^{z _{2-1}} 1 _{\R _{+}} , \end{align}

one has:

\Gamma(z _{1)} \Gamma(z _{2)} = \int _{\R} f(u)du ⋅ \int _{\R} g(u) du = \int _{\R} (f*g)(u)du = \Beta(z _{1,z} _{2)\Gamma(z} _{1+z} _{2).}

Differentiation of the beta function

We have

(\partial)/(\partial z _{1}) B (z _{1,} z _{2)} = B (z _{1,} z _{2)} \left( (\Gamma'(z _{1)})/(\Gamma(z _{1)}) - (\Gamma'(z _{1} + z _{2)})/(\Gamma(z _{1} + z _{2)}) \right) = B (z _{1,} z _{2)} ( \psi(z _{1)} - \psi(z _{1} + z _{2)),} (\partial)/(\partial z _{m}) B (z _{1,} z _{2,} ... , z _{n)} = B (z _{1,} z _{2,} ... , z _{n)} \left(\psi(z _{m)} - \psi\left( \sum ^{n} _{k =1} z _{k} \right)\right), 1 \le m\le n,

where

\psi(z) denotes the digamma function.

Approximation

Stirling's approximation gives the asymptotic formula

\Beta(x,y) \sim \sqrt { 2 \pi } (x ^{x} y ^{y})/(( {x + y) ) ^{x} }

for large and large .

If on the other hand is large and is fixed, then

\Beta(x,y) \sim \Gamma(y)x ^{-y} .

Other identities and formulas

The integral defining the beta function may be rewritten in a variety of ways, including the following:

\begin{align} \Beta(z _{1,z} _{2)} & = 2 \int ^{\pi / 2} _{0} (\sin\theta) ^{2 z _{1-1}} (\cos\theta) ^{2 z _{2-1}} d\theta, \\[ 6 pt] & = \int ^{inf ty (t ^{z _{1-1}})/(( 1+ t) ^{z _{1+z} _{2}})} _{0} dt, \\[ 6 pt] & = n\int ^{1 t} _{0} ^{nz _{1-1}} ( 1- t ^{n)} ^{z _{2-1}} dt, \\ & = ( 1- a) ^{z _{2}} \int ^{1} _{0} (( 1- t) ^{z _{1-1}} t ^{z _{2-1}})/(( 1- at) ^{z _{1+z} _{2}}) dt forany a\in R _{\leq} , \end{align}

where in the second-to-last identity is any positive real number. One may move from the first integral to the second one by substituting

t = \tan ^{2(\theta)} .

The beta function can be written as an infinite sum

\Beta(x,y) = \sum ^{inf ty} _{n =0} (( 1- x) _{n})/((y + n)n!)

(where

(x) _{n} is the rising factorial)and as an infinite product\Beta(x,y) = (x + y)/(x y) \prod ^{inf ty} _{n =1} \left( 1+ \dfrac{x y}{n (x + y + n)}\right) ^{-1} .

The beta function satisfies several identities analogous to corresponding identities for binomial coefficients, including a version of Pascal's identity

\Beta(x,y) = \Beta(x, y +1 ) + \Beta(x +1 , y)

and a simple recurrence on one coordinate:

\Beta(x +1 ,y) = \Beta(x, y) ⋅ \dfrac{x}{x + y}, \Beta(x,y +1 ) = \Beta(x, y) ⋅ \dfrac{y}{x + y}.

The positive integer values of the beta function are also the partial derivatives of a 2D function: for all nonnegative integers

m and n ,\Beta(m +1 , n +1 ) = (\partial ^{m+n} h)/(\partial a ^{m} \partial b ^{n}) ( 0 , 0 ), whereh(a, b) = (e ^{a-e} ^{b})/(a - b) . The Pascal-like identity above implies that this function is a solution to the first-order partial differential equationh = h _{a+h} _{b.}

For

x, y \geq 1 , the beta function may be written in terms of a convolution involving the truncated power function t \mapsto t ^{x} _{+} :\Beta(x,y) ⋅ \left(t \mapsto t ^{x + y -1} _{+} \right) = ( t \mapsto t ^{x -1} _{+} ) * ( t \mapsto t ^{y -1} _{+} )

Evaluations at particular points may simplify significantly; for example,

\Beta( 1 ,x) = \dfrac{ 1 }{x} and\Beta(x, 1- x) = \dfrac{\pi}{\sin(\pi x)}, x \not \in Z

By taking

x = (1)/(2) in this last formula, it follows that \Gamma( 1 / 2 ) = \sqrt{\pi} .Generalizing this into a bivariate identity for a product of beta functions leads to:\Beta(x,y) ⋅ \Beta(x + y, 1- y) = (\pi)/(x \sin(\pi y)) .

Euler's integral for the beta function may be converted into an integral over the Pochhammer contour as

\left( 1- e ^{2\pi} \right)\left( 1- e ^{2\pi} \right)\Beta(\alpha,\beta) = \int _{C} t ^{\alpha-1} ( 1- t) ^{\beta-1} dt.

This Pochhammer contour integral converges for all values of and and so gives the analytic continuation of the beta function.

Just as the gamma function for integers describes factorials, the beta function can define a binomial coefficient after adjusting indices:

\binom{n}{k} = (1)/((n +1 )\Beta(n - k +1 , k +1 )) .

Moreover, for integer, can be factored to give a closed form interpolation function for continuous values of :

\binom{n}{k} = ( -1 ) ^{n} n! ⋅ (\sin (\pi k))/(\pi \displaystyle\prod ^{n} _{i =0} (k - i)) .

Reciprocal beta function

The reciprocal beta function is the function about the form

f(x,y) = (1)/(\Beta(x,y))

Interestingly, their integral representations closely relate as the definite integral of trigonometric functions with product of its power and multiple-angle:

\int ^{\pi\sin} _{0} ^{x-1} \theta\sin y\theta~d\theta = (\pi\sin (y\pi)/(2))/(2 ^{x -1} x\Beta\left( (x + y +1)/(2) , (x - y +1)/(2) \right)) \int ^{\pi\sin} _{0} ^{x-1} \theta\cos y\theta~d\theta = (\pi\cos (y\pi)/(2))/(2 ^{x -1} x\Beta\left( (x + y +1)/(2) , (x - y +1)/(2) \right)) \int ^{\pi\cos} _{0} ^{x-1} \theta\sin y\theta~d\theta = (\pi\cos (y\pi)/(2))/(2 ^{x -1} x\Beta\left( (x + y +1)/(2) , (x - y +1)/(2) \right)) \int ^{(\pi)/(2)} _{0} \cos ^{x-1} \theta\cos y\theta~d\theta = (\pi)/(2 ^{xx\Beta\left( (x + y +1)/(2)} , (x - y +1)/(2) \right))

Incomplete beta function

The incomplete beta function, a generalization of the beta function, is defined as

\Beta(x;a,b) = \int ^{x} _{0} t ^{a-1} ( 1- t) ^{b-1} dt.

For, the incomplete beta function coincides with the complete beta function. The relationship between the two functions is like that between the gamma function and its generalization the incomplete gamma function. For positive integer a and b, the incomplete beta function will be a polynomial of degree a + b - 1 with rational coefficients.

By the substitution

t = \sin ^{2\theta} and t = (1 { 1+ s})/() , we show that \Beta(x;a,b) = 2 \int ^{\arcsin \sqrt x} _{0} \sin ^{2a-1} \theta\cos ^{2b-1} \theta d\theta = \int ^{inf ty} _{(1- x)/(x)} (s ^{a-1})/(( 1+ s) ^{a+b}) ds

The regularized incomplete beta function (or regularized beta function for short) is defined in terms of the incomplete beta function and the complete beta function:

I _{x(a,b)} = (\Beta(x;a,b))/(\Beta(a,b)) . F(k;n,p) of a random variable following a binomial distribution with probability of single success and number of Bernoulli trials :F(k;n,p) = \Pr\left(X \le k\right) = I _{1-p} (n - k, k +1 ) = 1 - I _{p(k+1,n-k).}

Properties

\begin{align} I _{0(a,b)} & = 0 \\ I _{1(a,b)} & = 1 \\ I _{x(a,1)} & = x ^{a\\ I} _{x( 1 ,b)} & = 1 - ( 1- x) ^{b} \\ I _{x(a,b)} & = 1 - I _{1-x} (b,a) \\ I _{x(a+1,b)} & = I _{x(a,b) - (x ^{a(1-x)} ^{b})/(a \Beta(a,b))} \\ I _{x(a,b+1)} & = I _{x(a,b) + (x ^{a(1-x)} ^{b})/(b \Beta(a,b))} \\ \int \Beta(x;a,b) d x & = x \Beta(x; a, b) - \Beta(x; a +1 , b) \\ \Beta(x;a,b)& = ( -1 ) ^{a} \Beta\left( (x)/(x -1) ;a, 1- a - b\right) \end{align}

Continued fraction expansion

The continued fraction expansion

\Beta(x;a,b) = (x ^{a} ( 1- x) ^{b})/(a \left( 1 + ({d)/(_{1})) { 1 + } ({d)/(_{2}) } \frac \frac \cdots \right)}

with odd and even coefficients respectively

{d} _{2} = - ((a + m) (a + b + m) x)/((a +2 m) (a +2 m +1 )) {d} _{2} = (m (b - m) x)/((a +2 m -1 ) (a +2 m))

converges rapidly when

x is not close to 1. The 4 m and 4 m + 1 convergents are less than \Beta(x;a,b) , while the 4 m + 2 and 4 m + 3 convergents are greater than \Beta(x;a,b) .

For

x > (a +1)/(a + b +2) , the function may be evaluated more efficiently using \Beta(x;a,b) = \Beta(a, b) - \Beta( 1 - x;b,a) .

Multivariate beta function

The beta function can be extended to a function with more than two arguments:

\Beta(\alpha _{1,\alpha} _{2,\ldots\alpha} _{n)} = (\Gamma(\alpha _{1)\Gamma(\alpha} _{2)} … \Gamma(\alpha _{n)})/(\Gamma(\alpha _{1} + \alpha _{2} + … + \alpha _{n)}) .

This multivariate beta function is used in the definition of the Dirichlet distribution. Its relationship to the beta function is analogous to the relationship between multinomial coefficients and binomial coefficients. For example, it satisfies a similar version of Pascal's identity:

\Beta(\alpha _{1,\alpha} _{2,\ldots\alpha} _{n)} = \Beta(\alpha _{1+1,\alpha} _{2,\ldots\alpha} _{n)+\Beta(\alpha} _{1,\alpha} _{2+1,\ldots\alpha} _{n)+ … +\Beta(\alpha} _{1,\alpha} _{2,\ldots\alpha} _{n+1)} .

Applications

The beta function is useful in computing and representing the scattering amplitude for Regge trajectories. Furthermore, it was the first known scattering amplitude in string theory, first conjectured by Gabriele Veneziano. It also occurs in the theory of the preferential attachment process, a type of stochastic urn process. The beta function is also important in statistics, e.g. for the beta distribution and beta prime distribution. As briefly alluded to previously, the beta function is closely tied with the gamma function and plays an important role in calculus.

Software implementation

Even if unavailable directly, the complete and incomplete beta function values can be calculated using functions commonly included in spreadsheet or computer algebra systems.

In Microsoft Excel, for example, the complete beta function can be computed with the [[Gamma_function#The_log-gamma_function|GammaLn]] function (or special.gammaln in Python's SciPy package):

Value = Exp(GammaLn(a) + GammaLn(b) − GammaLn(a + b))

This result follows from the properties listed above.

The incomplete beta function cannot be directly computed using such relations and other methods must be used. In GNU Octave, it is computed using a continued fraction expansion.

The incomplete beta function has existing implementation in common languages. For instance, betainc (incomplete beta function) in MATLAB and GNU Octave, pbeta (probability of beta distribution) in R, or special.betainc in SciPy compute the regularized incomplete beta function—which is, in fact, the cumulative beta distribution—and so, to get the actual incomplete beta function, one must multiply the result of betainc by the result returned by the corresponding beta function. In Mathematica, Beta[x, a, b] and BetaRegularized[x, a, b] give

\Beta(x;a,b) and I _{x(a,b)} , respectively.

See also

• Beta distribution and Beta prime distribution, two probability distributions related to the beta function
• Jacobi sum, the analogue of the beta function over finite fields.
• Nørlund–Rice integral
• Yule–Simon distribution

External links

• • • Arbitrarily accurate values can be obtained from:
  • The Wolfram functions site: Evaluate Beta Regularized incomplete beta
  • danielsoper.com: Incomplete beta function calculator, Regularized incomplete beta function calculator

Notes and References

1. . Specifically, see 6.2 Beta Function.