Open In App

Mathematics | Beta Distribution Model

Last Updated : 18 Feb, 2021
Improve
Improve
Like Article
Like
Save
Share
Report

Introduction :

Suppose an event can occur several times within a given unit of time. When the total number of occurrences of the event is unknown, we can think of it as a random variable. When a random variable X takes on values on the interval from 0 to 1, one choice of a probability density is the beta distribution whose probability density function is given by as follows.

Representation of probability density function –

f(x) = \frac{1}{\Beta(m.n)}x^{m-1}(1-x)^{n-1}

It will be applicable only when the given below given condition will pass.

x > 0, m >0, n >0.  
f(x) = 0 , Otherwise

Here, you will see the meaning of function as you have shown in representation of probability density function where, B(m.n) is the value of the beta function.

Representation of B(m.n) –

\Beta(m,n) = \int^{1}_{0} x^{m- 1}(1-x)^{n-1}dx = \Beta(n,m)

Integrating it by parts, we will get the following expression as given below. 

\Beta(m,n) = \frac{\Gamma(m)\Gamma(n)}{\Gamma(m+n)}

Where, Γ(x) is the gamma function of x, calculated as –

Γ(x) = (x-1)Γ(x-1)    for x>1
=> \Gamma(x) = (x-1)!    when x is an integer

The random variable X is represented as follows.

Representation of  random variable X –

X ~ BETA(m,n)

Expected Value:

The Expected Value of the Beta distribution can be found by summing up products of Values with their respective probabilities.

\mu = E(X) = \int^{\infty}_{-\infty} x.f(x) dx
\mu = \frac{1}{\Beta(m,n)} \int^{1}_{0}x. x^{m - 1}(1-x)^{n-1} dx
\mu = \frac{1}{\Beta(m,n)} \int^{1}_{0} x^{m}(1-x)^{n-1} dx

Upon using the value of beta function, we will get the following expression a follows.

\mu = \frac{1}{\Beta(m,n)} \Beta(m+1,n)
but, \Beta(m,n) = \frac{\Gamma(m)\Gamma(n)}{\Gamma(m+n)}
also, \Beta(m+1,n) = \frac{\Gamma(m+1)\Gamma(n)}{\Gamma(m+n+1)}
=> \mu = \frac{\Gamma(m+1)\Gamma(n)}{\Gamma(m+n+1)} * \frac{\Gamma(m+n)}{\Gamma(m)\Gamma(n)}

Upon using the property of gamma function, that is Γ(x) = (x-1)!, we will get the following expression as follows.

 So, E(x) = \frac{m! (n-1)! (m+n-1)!}{(m-1)! (n-1)! (m+n)!}
here, m, n -> integers

   where provided m and n are integers.

Variance and Standard Deviation :

The Variance of the Beta distribution can be found using the Variance Formula.

σ^2 = E( X − μ )^2 = E( X^2 ) − μ^2
E(X^2) = \int^{\infty}_{-\infty} x^2.f(x) dx
E(X^2) = \int^{1}_{0}\frac{1}{\Beta(m,n)}x^2. x^{m - 1}(1-x)^{n-1} dx
E(X^2) = \frac{1}{\Beta(m,n)}\int^{1}_{0}x^{m + 1}(1-x)^{n-1} dx

Upon using the value of beta function, we will get the following expression as follows.

E(X^2) = \frac{1}{\Beta(m,n)} \Beta(m+2,n+2)
but, \Beta(m,n) = \frac{\Gamma(m)\Gamma(n)}{\Gamma(m+n)}
also, \Beta(m+1,n) = \frac{\Gamma(m+2)\Gamma(n)}{\Gamma(m+n+2)}
=> \mu = \frac{\Gamma(m+2)\Gamma(n)}{\Gamma(m+n+2)} * \frac{\Gamma(m+n)}{\Gamma(m)\Gamma(n)}

Upon using the property of gamma function, that is Γ(x) = (x-1)!, we will get the following expression as follows.

E(X^2) = \frac{(m+1)!(n-1)!(m+n-1)!}{(m+n+1)!(m-1)!(n-1)!}
E(X^2) = \frac{m(m+1)}{(m+n+1)(m+n)}
So, Var(X) = E(X^2) - \mu^2
Var(X) = \frac{m(m+1)}{(m+n+1)(m+n)} - \frac{m^2}{(m+n)^2}
Var(X) = \sigma^2= \frac{m.n}{(m+n+1)(m+n)^2}

Standard Deviation is given by as follows.

\sigma = \sqrt{\frac{m.n}{(m+n+1)(m+n)^2}} = \frac{1}{(m+ n)}\sqrt{\frac{m.n}{(m+n+1)}}

Example – 

In a certain county, the proportion of highway sections requiring repairs in any given year is a random variable having the beta distribution with m = 3 and n = 2.

(a) On the average, what percentage of the highway sections require repairs in any given year?

(b) Find the probability that at most half of the highway sections will require repairs in any given year.

Solution –

(a) \mu = \frac{3}{3+2} = 0.60

which means that on the average 60% of the highway sections require repairs in any given year.

(b) Substituting m=3 and n=2 in the probability density function, and 
using Γ(3) = 2! = 2, Γ(2) = 1! = 1, Γ(5) = 4! = 24, 

we will get the following expression as follows.

f(x) = 12 x^2 (1 − x) for 0<x<1. Otherwise f(x) = 0

Thus, the desired probability is given by as follows.

\int^{\frac{1}{2}}_0 12x^2 (1-x) dx = \frac{5}{16}


Like Article
Suggest improvement
Previous
Next
Share your thoughts in the comments

Similar Reads