Applications of statistics often require working with

for large *x* and fixed α and β.

The simplest approximation to this ratio of gamma functions is

You can do a little better [1] with

where *c* = (α + β – 1)/2.

To get more accuracy we’ll need more terms. Tricomi and Erdélyi worked out the full asymptotic series 70 years ago:

Here *B*_{n}^{λ} is the “generalized Bernoulli polynomial,” a.k.a. the Nørlund polynomial [3].

There’s no simple way to state the coefficients in the series above. The expression in terms of Nørlund polynomials is compact, but then the Nørlund polynomials are a bit complicated to define. The can be calculated in Mathematica with `NorlundB`

, and so the coefficients would be

Binomial[α-β, n] NorlundB[n, α-β+1, α]

in Mathematica.

Let’s see how the various approximations compare in an example.

{α, β} = {3.4, 5.6} f[x_] := Gamma[x+α]/Gamma[x+β] f1[x_] = x^(α-β) f2[x_] = (x + (α+β-1)/2)^(α-β) f3[x_] = Sum[Binomial[α-β, n] NorlundB[n, α-β+1, α] x^{α-β-n}, {n, 0, 2}]

We can plot the exact function `f`

and its approximations with

LogPlot[{f[x], f1[x], f2[x], f3[x]}, {x, 10, 15}, PlotLabels -> Automatic]

This gives the following.

It’s interesting that the approximation `f2`

works really well and is the most accurate approximation over the range shown. For large enough *x*, say *x* = 1000, the series approximation `f3`

will be more accurate.

***

[1] A. Laforgia, P. Natalini. On the asymptotic expansion of a ratio of gamma functions. Journal of Mathematical Analysis and Applications. Volume 389, Issue 2, 15 May 2012, pp. 833-837.

[2] F. G. Tricomi, A. Erdélyi. The asymptotic expansion of a ratio of gamma functions. Pacific J. Math. 1 (1951), pp. 133-142.

[3] The authors in [1] quote this formula from [2] but leave out the binomial coefficient. Note that the binomial coefficient parameters are not necessarily positive integers and so the binomial coefficients should be interpreted in the sense of the most general definition.