Dirichlet average

Dirichlet averages are averages of functions under the Dirichlet density. An important one are dirichlet averages that have a certain argument structure, namely

${\displaystyle F(\mathbf {b}$ ;\mathbf {z} )=\int f(\mathbf {u} \cdot \mathbf {z} )\,d\mu _{b}(\mathbf {u} ),}

where ${\displaystyle \mathbf {u} \cdot \mathbf {z} =\sum _{i}^{N}u_{i}\cdot z_{i}}$ and ${\displaystyle d\mu _{b}(\mathbf {u} )=u_{1}^{b_{1}-1}\cdots u_{N}^{b_{N}-1}d\mathbf {u} }$ is the Dirichlet measure with dimension N. They were introduced by the mathematician Bille C. Carlson in the '70s who noticed that the simple notion of this type of averaging generalizes and unifies many special functions, among them generalized hypergeometric functions or various orthogonal polynomials:[1]. They also play an important role for the solution of elliptic integrals (see Carlson symmetric form) and are connected to statistical applications in various ways, for example in Bayesian analysis [2].

Notable Dirichlet averages

Some Dirichlet averages are so fundamental that they are named. A few are listed below.

R-function

The (Carlson) R-function is the Dirichlet average of ${\displaystyle x^{n}}$,

${\displaystyle R_{n}(\mathbf {b} ,\mathbf {z} )=\int (\mathbf {u} \cdot \mathbf {z} )^{n}\,d\mu _{b}(\mathbf {u} )}$

with ${\displaystyle n}$. Sometimes ${\displaystyle R_{n}(\mathbf {b} ,\mathbf {z} )}$ is also denoted by ${\displaystyle R(-n;\mathbf {b} ,\mathbf {z} )}$.

Exact solutions:

For ${\displaystyle n\geq 0,n\in \mathbb {N} }$ it is possible to write an exact solution in the form of an iterative sum[3]

${\displaystyle R_{n}(\mathbf {b} ,\mathbf {z} )={\frac {\Gamma (n+1)\Gamma (b)}{\Gamma (b+n)}}\cdot D_{n}{\text{ with }}D_{n}={\frac {1}{n}}\sum _{k=1}^{n}\left(\sum _{i=1}^{N}b_{i}\cdot z_{i}^{k}\right)\cdot D_{n-k}}$

where ${\displaystyle D_{0}=1}$, ${\displaystyle N}$ is the dimension of ${\displaystyle \mathbf {b} }$ or ${\displaystyle \mathbf {z} }$ and ${\displaystyle b=\sum b_{i}}$.

S-function

The (Carlson) S-function is the Dirichlet average of ${\displaystyle e^{x}}$,

${\displaystyle S(\mathbf {b} ,\mathbf {z} )=\int \exp(\mathbf {u} \cdot \mathbf {z} )\,d\mu _{b}(\mathbf {u} ).}$

References

1. Carlson, B.C. (1977). Special functions of applied mathematics.
2. Dickey, J.M. (1983). "Multiple hypergeometric functions: Probabilistic interpretations and statistical uses". Journal of the American Statistical Association. 78 (383): 628. doi:10.2307/2288131.
3. Glüsenkamp, T. (2018). "Probabilistic treatment of the uncertainty from the finite size of weighted Monte Carlo data". EPJ Plus. 133 (6): 218. arXiv:1712.01293. doi:10.1140/epjp/i2018-12042-x.