Spectral Gap for Measure-Valued Diffusion Processes
Panpan Ren, Feng-Yu Wang

TL;DR
This paper estimates the spectral gap for measure-valued diffusion processes on Riemannian manifolds, providing explicit convergence rates to distributions relevant in population genetics.
Contribution
It introduces explicit spectral gap estimates for measure-valued diffusions, enhancing understanding of their convergence behavior on Riemannian manifolds.
Findings
Explicit spectral gap estimates derived.
Convergence rates to Dirichlet and Gamma distributions established.
Applications to population genetics models.
Abstract
The spectral gap is estimated for measure-valued diffusion processes induced by the intrinsic/extrinsic derivatives on the space of finite measures over a Riemannian manifold. This provides explicit exponential convergence rate for these processes to approximate the Dirichlet and Gamma distributions arising from population genetics.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsBayesian Methods and Mixture Models · Genetic and phenotypic traits in livestock · Genetic Mapping and Diversity in Plants and Animals
**Spectral Gap for Measure-Valued Diffusion Processes 111Supported in
part by NNSFC (11771326, 11831014, 11726627).**
**Panpan Renb,c), Feng-Yu Wanga,b)
a)** Center for Applied Mathematics, Tianjin University, Tianjin 300072, China
b) Department of Mathematics, Swansea University, Singleton Park, SA2 8PP, United Kingdom
c) Mathematical Institute,Woodstock Road, OX2 6GG, University of Oxford
[email protected], [email protected] [email protected], [email protected]
Abstract
The spectral gap is estimated for some measure-valued processes, which are induced by the intrinsic/extrinsic derivatives on the space of finite measures over a Riemannian manifold. These processes are symmetric with respect to the Dirichlet and Gamma distributions arising from population genetics. In addition to the evolution of allelic frequencies investigated in the literature, they also describe stochastic movements of individuals.
AMS subject Classification: 60J60, 58J65.
Keywords: Extrinsic derivative, weighted Gamma distribution, Poincaré inequality, weak Poincaré inequality, super Poincaré inequality.
1 Introduction
The Dirichlet distribution arises naturally in Bayesian inference as conjugate priors for categorical distribution and infinite non-parametric discrete distributions respectively. In population genetics, it describes the distribution of allelic frequencies (see for instance [3, 10, 14]). To simulate the Dirichlet distribution using stochastic dynamic systems, some diffusion processes generalized from the Wright-Fisher diffusion have been considered, see [4, 5, 6, 7, 20] and references within. In this paper, we investigate diffusion processes induced by the Dirichlet distribution and the intrinsic/extrinsic derivatives, where the extrinsic derivative term determines the evolution of allelic frequencies, and the intrinsic derivative drives the movement of individuals.
In the following three subsections, we introduce the reference measures, intrinsic and extrinsic derivatives, and the main results of the paper respectively. We will take the notation for a measurable space and .
1.1 Reference measures
Let be a complete Riemannian manifold. Consider the space of all nonnegative finite measures on , and let be the set of all probability measures on . According to [13, Theorem 3.2], both spaces are Polish under the weak topology. In general, for an ergodic Markov process with stationary distribution , one may simulate by using the empirical measures where is the Dirac measure at point . In practice, one may also approximate using the discrete time empirical measures
[TABLE]
See for instance [8] and references within for the study of the convergence rate.
For , the Dirichlet distribution with shape is the unique probability measure on such that for any measurable partition of ,
[TABLE]
obeys the Dirichlet distribution with parameter . Recall that for any , the Dirichlet distribution with parameter is the following probability measure on the simplex :
[TABLE]
where in case we set , and denotes the Dirac measure at point in a measurable space. If refers to the distribution of population on , then under the state , stand for the proportions of population located in the areas respectively.
We will also consider the Gamma distribution on with shape , whose marginal distribution on coincides with the Dirichlet distribution . Recall that is the unique probability measure on such that for any finitely many disjoint measurable subsets of ,
[TABLE]
are independent Gamma random variables with shape parameters and scale parameter ; that is,
[TABLE]
holds for any where for a constant ,
[TABLE]
and we set , the Dirac measure at point [math].
The Gamma distribution is supported on the class of finite discrete measures
[TABLE]
Moreover, under the random variables and are independent with
[TABLE]
Consequently,
[TABLE]
Both and are images of the Poisson measure with intensity
[TABLE]
on the product manifold . Recall that is the unique probability measure on the configuration space
[TABLE]
equipped with the vague topology, such that for any disjoint compact subsets of , are independent Poisson random variables of parameters By [9, Theorem 6.2], we have
[TABLE]
and
[TABLE]
where
[TABLE]
Combining (1.4) with (1.7), we obtain
[TABLE]
1.2 Intrinsic and extrinsic derivatives
These derivatives were introduced in [1] and [15] on the configuration space and the space of probability measures respectively, which can be extended to under the map , see for instance [11].
To introduce the intrinsic derivative for a function on (or ), we let be the class of smooth vector fields with compact supports on . For any , let
[TABLE]
where is the exponential map on . Then . For a function on , we define its directional derivative along by
[TABLE]
if it exists. Let be the space of all measurable vector fields on with . When exists for all such that
[TABLE]
holds for some constant , then by Riesz representation theorem there exists a unique such that
[TABLE]
In this case, we call intrinsically differentiable at with derivative . If is intrinsically differentiable at all (or ), we call it intrinsically differentiable on (or ).
Next, a measurable real function on is called extrinsically differentiable at , if
[TABLE]
such that
[TABLE]
When a function on is considered, it is called intrinsically differentiable if
[TABLE]
and . If is extrinsically differentiable at all (or ), we call it extrinsically differentiable on (or ). Let (respectively ) denote the set of functions which are intrinsically and extrinsically differentiable on (respectively ).
A typical subclass of and is the set of cylindrical functions
[TABLE]
where This class is dense in and , and the cylindrical function is differentiable with
[TABLE]
where is the gradient operator on . Restricting on we will consider
[TABLE]
which is the centered extrinsic derivative of at since
[TABLE]
See [16] for general results on the relations of and
1.3 The main result
Now, for any , we consider the square fields for :
[TABLE]
which lead to the following bilinear forms on and respectively:
[TABLE]
To ensure the closability of these bilinear forms, we take
[TABLE]
where is the Riemannian volume measure. Then the integration by parts formula gives
[TABLE]
So, the bilinear form is closable in , and the closure is a Dirichlet form.
We will prove the closability of and , and calculate the spectral gaps for the corresponding Dirichlet forms.
Recall that for a probability space and a symmetric Dirichlet form on with and , the spectral gap of the Dirichlet form is given by
[TABLE]
By the spectral theorem, is the exponential convergence rate of the associated Markov semigroup , i.e.
[TABLE]
Let
[TABLE]
be the spectral gap of the Dirichlet form in . The main result of this paper is the following.
Theorem 1.1**.**
Let be in and let .
* is closable in and the closure is a quasi-regular symmetric Dirichlet form with spectral gap * 2.
* is closable in and the closure is a quasi-regular symmetric Dirichlet form with spectral gap satisfying*
[TABLE]
Consequently, if then Moreover, for and for any ,
[TABLE]
The formula (1.17) will play a crucial role in the proof of the closability of . When , the exact value of is unknown. Since in this case the intrinsic derivative part will play a non-trivial role, we believe that is strictly larger than hopefully our upper bound could be sharp.
When and without the intrinsic derivative part, the Dirichlet form reduces to
[TABLE]
which is associated with the Fleming-Viot process. It has been derived in [18] that
[TABLE]
see [20] for the study of log-Sobolev inequality in finite-dimensions as well as [5, 21, 22] for functional inequalities of a modified Fleming-Viot process.
When and , [11, Theorem 14] presents an integration by parts formula for in the class
[TABLE]
where
[TABLE]
Therefore, letting be the generator of the Dirichlet form , we have However, in general is not included in . Indeed, according to [11] we have the integration by parts formula
[TABLE]
for B_{\varepsilon,v}(\eta):=\sum_{\eta(\{x\})\geq\varepsilon}\big{\{}{\rm div}(v)+\langle v,\nabla V\rangle\big{\}}(x), where is the divergence operator in . This formula makes sense because
[TABLE]
So, for any and . However, since , (1.20) does not make sense for .
To prove Theorem 1.1, we will formulate the bilinear form as the image of the Dirichlet form on the configuration space constructed in [1], for which the (weak) Poincaré inequality has been established in [17, Section 7]. To this end, we first recall in Section 2 some known results on the configuration space, then prove Theorem 1.1(1) in Section 3 by transforming these results to the Gamma process on , and finally prove Theorem 1.1(2) in Section 4 by mapping the Gamma process to the subclass .
2 Analysis on the configuration space
In this section, we first recall the diffusion process on the configuration space constructed in [1, 2], then calculate the spectral gap.
For where and , let
[TABLE]
For , we take the following Riemannian metric on the manifold
[TABLE]
Let and be the Laplacian, gradient and volume measure on respectively.
Consider the bilinear form
[TABLE]
To formulate the integration by parts formula of this form, we intend to find out a function on such that
[TABLE]
where is given by (1.5). So, for
[TABLE]
we have the integration by parts formula
[TABLE]
Therefore, letting
[TABLE]
we have (see [2, Theorem 4.3])
[TABLE]
which implies the closability of , so that the closure is a symmetric Dirichlet form in . Moreover, as explained in [13, Section 4.5.1], the result [13, Corollary 4.9] applies to this situation, so that the Dirichlet form is quasi-regular and local, and hence is associated with a diffusion process on . Recall that is equipped with the vague topology.
To calculate the generator defined in (2.7), it suffices to figure out the operator To this end, for a fixed point , we take the normal coordinates in a neighbourhood of such that
[TABLE]
satisfy
[TABLE]
Let
[TABLE]
Then the Riemannian volume measure on is
[TABLE]
so that (2.4) holds for
[TABLE]
Letting , we derive from (2.9) that . So,
[TABLE]
By the same reason, for any we have
[TABLE]
This together with (2.10) implies that at point ,
[TABLE]
Since is arbitrary, this formula holds for all points
Theorem 2.1**.**
Let and be as in . Then
Proof.
According to [19], see also [17, Theorem 7.1], we have
[TABLE]
where is the closure of under the Sobolev norm Let
[TABLE]
By (2.11) we have
[TABLE]
Combining this with (2.6), for any we have
[TABLE]
Therefore,
On the other hand, since is complete, there exists a sequence such that
[TABLE]
For any let
[TABLE]
Then with So, for fixed , by the dominated convergence theorem
[TABLE]
[TABLE]
Thus, with
[TABLE]
Combining this with
[TABLE]
we obtain
[TABLE]
This together with derived above gives So, the proof is finished by (2.12). ∎
3 Proof of Theorem 1.1(1)
Theorem 3.1**.**
Let be as in . Then is closable in and the closure is a quasi-regular symmetric Dirichlet form. Moreover,
Proof.
Due to (1.7), we may prove this result by using (2.8) and Theorem 2.1. For any , let Then
[TABLE]
Let . We have
[TABLE]
such that
[TABLE]
Because of (2.2) and (3.1), we have
[TABLE]
Thus,
[TABLE]
Combining this with (1.7), (1.14) and (2.3), we obtain
[TABLE]
Below we prove the closability, quasi-regularity, and the spectral gap bounds respectively.
(a) The closability. Let such that
[TABLE]
It remains to show that By (1.7) and (3.3), (3.4) implies
[TABLE]
so that the closability of and (3.4) imply
[TABLE]
(b) The quasi-regularity. According to [12, Chap. IV, Def. 3.1], the Dirichlet form is quasi-regular if and only if there exist a sequences of compact subsets of such that the class
[TABLE]
is dense in under the Sobolev norm
[TABLE]
In this case, the sequence is called a -nest.
To construct such a -nest, we choose a function with as and such that , where is the Riemannian distance from a fixed point . Thus,
[TABLE]
Since is complete and as , the level sets
[TABLE]
are compact in . It remains to show that is a -nest. Since is dense in , it suffices to show that for any , there exist a sequence , where is defined in (4.3) for the present , such that
[TABLE]
We will prove this formula for
[TABLE]
To this end, we first confirm that , so that by the definitions of and . By (3.6), there exist functions such that
[TABLE]
Noting that obeys the Gamma-distribution with parameter , we have
[TABLE]
By the dominated convergence theorem and (3.9), the functions satisfy
[TABLE]
Moreover, (1.11), (1.13), (1.14) and (3.9) imply
[TABLE]
Therefore, Now, let be defined in (3.8). Then in , and (1.11), (1.13) and (1.14) imply
[TABLE]
where Moreover, the dominated convergence theorem implies that as . Therefore, (3.7) holds as desired.
(c) Spectral gap estimates. By Theorem 2.1 and (3.3), for any with , we have
[TABLE]
This implies On the other hand, we take and intend to show that with
[TABLE]
so that by definition and hence the proof is finished.
Let such that for For in (2.14), define
[TABLE]
Then and as . By (1.3) we have
[TABLE]
So, the dominated convergence theorem implies as , and
[TABLE]
Therefor, with
[TABLE]
Combining this with (3.12) we derive (3.10) and hence finish the proof.
∎
4 Proof of Theorem 1.1(2)
The quasi-regularity can be proved by the same means as in the step (b) of the proof of Theorem 1.1(1). So, we need only to prove the closability, the formula (1.17), and the claimed spectral gap bounds.
(1) Closability and (1.17). Let . It suffices to prove that for any , one has such that (1.17) holds. Indeed, since is closed, (1.17) implies the closability of .
Similarly to (b) in the proof of Theorem 3.1, and noting that , we see that for any , in with
[TABLE]
However, since for general we do not have , this does not imply as desired.
To approximate using functions in , we write for some and . Since the Riemannian manifold is complete, we may construct such that
[TABLE]
For any let
[TABLE]
So, (4.1) implies . Obviously, . By (4.2) and noting that
[TABLE]
we may find out a constant such that
[TABLE]
Since , by (4.1) and Fatou’s lemma, we arrive at
[TABLE]
Then . Similarly, for , we define in the same way. By the above formulas of intrinsic and extrinsic derivatives for and , it is easy to see that
[TABLE]
and the same holds for replacing . Therefore, by the dominated convergence theorem and using (1.3) and (1.4), we obtain
[TABLE]
Therefore, (1.17) holds.
(2) Spectral gap estimates. Since
[TABLE]
where is given in (1.18) as the Dirichlet form of the Fleming-Viot process, it follows from (1.19) that
[TABLE]
To verify the spectral gap upper bound, for any with and , let Then It suffices to show that
[TABLE]
To this end, we recall that (see for instance the proof of [17, Lemma 7.2])
[TABLE]
Letting for and applying (1.3), (1.4), and (1.7), we deduce from (4.4) that
[TABLE]
and similarly,
[TABLE]
Thus, for with and , we have
[TABLE]
and
[TABLE]
Therefore, (4.3) holds.
Acknowledgement.
We would like to thank the referees for helpful comments on an earlier version of the paper.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] S. Albeverio,Y. G. Kondratiev, M. Röckner, Differential geometry of Poisson spaces, C. R. Acad. Sci. Paris Sér. I Math. 323(1996), 1129–1134.
- 2[2] S. Albeverio,Y. G. Kondratiev, M. Röckner, Analysis and geometry on configuration spaces, J. Funct. Anal. 154(1998), 444–500.
- 3[3] R. J. Connor, J. E. Mosimann, Concepts of independence for proportions with a generalization of the Dirichlet distribution, J. Amer. Statist. Assoc. 64(1969), 194-206.
- 4[4] C. L. Epstein, R. Mazzeo, Wright-Fisher diffusion in one dimension, SIAM J. Math. Anal. 42(2010), 568–608.
- 5[5] S. Feng, L. Miclo, F.-Y. Wang, Poincaré inequality for Dirichlet distributions and infinite-dimensional generalizations, Lat. Am. J. Probab. Math. Stat. 14(2017), 361–380.
- 6[6] S. Feng, F.-Y. Wang, A class of infinite-dimensional diffusion processes with connection to population genetics, J. Appl. Probab. 44(2007), 938–949.
- 7[7] S. Feng, F.-Y. Wang, Harnack inequality and applications for infinite-dimensional GEM processes, Potential Anal. 44(2016), 137–153
- 8[8] F. Q. Gao, Long time asymptotics of unbounded additive functionals of Markov processes, Electron J. Probab. 22(2017), 1–21.
