English

A stochastic diffusion process for the Dirichlet distribution

Mathematical Physics 2013-03-05 v2 math.MP Probability Data Analysis, Statistics and Probability

Abstract

The method of potential solutions of Fokker-Planck equations is used to develop a transport equation for the joint probability of N coupled stochastic variables with the Dirichlet distribution as its asymptotic solution. To ensure a bounded sample space, a coupled nonlinear diffusion process is required: the Wiener-processes in the equivalent system of stochastic differential equations are multiplicative with coefficients dependent on all the stochastic variables. Individual samples of a discrete ensemble, obtained from the stochastic process, satisfy a unit-sum constraint at all times. The process may be used to represent realizations of a fluctuating ensemble of N variables subject to a conservation principle. Similar to the multivariate Wright-Fisher process, whose invariant is also Dirichlet, the univariate case yields a process whose invariant is the beta distribution. As a test of the results, Monte-Carlo simulations are used to evolve numerical ensembles toward the invariant Dirichlet distribution.

Keywords

Cite

@article{arxiv.1303.0217,
  title  = {A stochastic diffusion process for the Dirichlet distribution},
  author = {J. Bakosi and J. R. Ristorcelli},
  journal= {arXiv preprint arXiv:1303.0217},
  year   = {2013}
}

Comments

Accepted in International Journal of Stochastic Analysis, March 1, 2013

R2 v1 2026-06-21T23:35:06.714Z