A proportion or probability constrained to [0, 1], with a shape set by two counts.
Parameters and support
α
first shape parameter - α>0
β
second shape parameter - β>0
Support
0<x<1
The formulas
f(x)
f(x)=Γ(α)Γ(β)Γ(α+β)xα−1(1−x)β−1
F(x)
F(x)=Ix(α,β)
The regularized incomplete beta function. Elementary only when α or β is a small positive integer.
Mean
E[X]=α+βα
Variance
Var(X)=(α+β)2(α+β+1)αβ
MGF
Not examinable for this distribution.
The MGF is a confluent hypergeometric series and is not examinable. Use E[Xᵏ] = ∏(α+j)/(α+β+j) for j = 0…k−1.
Memory hook. Read α and β as pseudo-counts of successes and failures: the mean α/(α+β) is exactly the sample proportion they imply. Beta(1,1) is the standard uniform.
Where the moments come from
The normalising constant is 1/B(α,β) with B(α,β)=Γ(α+β)Γ(α)Γ(β).
E[Xk]=B(α,β)B(α+k,β)=Γ(α)Γ(α+β+k)Γ(α+k)Γ(α+β).
k=1 gives E[X]=α+βα.
k=2 gives E[X2]=(α+β)(α+β+1)α(α+1).
Subtracting the squared mean collapses to (α+β)2(α+β+1)αβ.
Worked example
The proportion of a portfolio that lapses in a year follows a Beta distribution with α = 2 and β = 8. Find the mean and variance of the lapse proportion.
E[X]=α+βα=102=0.2.
Var(X)=(α+β)2(α+β+1)αβ=102×112×8.
=110016=0.014545.
The standard deviation is 0.014545=0.1206.
Answer:0.014545
The mean, variance, CDF and moment generating function above are re-derived numerically from this distribution’s own density on every test run - summed over the support for a discrete distribution, integrated by quadrature for a continuous one - and compared with the closed forms printed here. A typo on this page fails that test run.
Traps
Writing (α+β+1) as (α+β)+1 in the wrong place - the cube-like denominator is (α+β)²(α+β+1).
Forgetting that Beta(1,1) is uniform, which is often the fastest sanity check on an answer.