arxiv: v1 [math.pr] 20 Dec 2017

Transcription

1 On the time to absorption in Λ-coalescents Götz Kersting and Anton Wakolbinger arxiv: v [math.pr] 2 Dec 27 Abstract We present a law of large numbers and a central limit theorem for the time to absorption of Λ-coalescents, started from n blocks, as n. The proofs rely on an approximation of the logarithm of the block-counting process of Λ-coalescents with a dust component by means of a drifted subordinator. AMS 2 subject classification: 6J75 primary, 6J27, 6F5 secondary Keywords: coalescents, time to absorption, law of large numbers, central limit theorem, subordinator with drift Introduction and main results How long does it take for the ancestral lineages of a large sample of individuals back to its common ancestor? For population of constant size this turns into a question on the absorption time of a coalescents, which describes the genealogical tree of n individuals by means of merging partitions. Here we consider coalescent with multiple mergers, also known as Λ-coalescents, which were introduced in 999 by Pitman [6] and Sagitov [7]. If Λ is a finite, non-zero measure on [,], then the Λ-coalescent started with n blocks is a continuous-time Markov chain Π n t,t taking its values in the set of partitions of {,...,n}. It has the property that whenever there are b blocks, each possible transition that involves merging k 2 of the blocks into a single block happens at rate λ b,k = p k b k Λdp p p 2, [,] and these are the only possible transitions. Let N n t be the number of blocks in the partition Π n t, t. Then τ n := inf{t : N n t = } is the time of the last merger, also called the absorption time of the coalescent started in n blocks. We will investigate the asymptotic distribution of τ n as n. Our first result is a law of large numbers for the times τ n. Let µ := [,] log Λdp p p 2, in particular µ = in case of Λ{} > or Λ{} >. Institut für Mathematik, Goethe Universität, Frankfurt am Main, Germany kersting@math.uni-frankfurt.de, wakolbinger@math.uni-frankfurt.de Work partially supported by the DFG Priority Programme SPP 59 Probabilistic Structures in Evolution

2 Theorem. For any Λ-coalescent, τ n logn µ. in probability as n. This theorem says that in a Λ-coalescent the number of blocks decays at least at an exponential rate. If µ =, then the right-hand limit is, and the coalescent decreases even super-exponentially fast. The case µ < is equivalently captured by the simultaneous validity of the conditions [,] Λdp p < and [,] log p Λdp <. The first one is a requirement on Λ in the vicinity of, it prohibits a swarm of small mergers as they occur in coalescents coming down from infinity, meaning that the τ n are bounded in probability uniformly in n. The second is a condition on Λ in the vicinity of. It rules out the possibility of mergers which, although appearing only every now and then, are so vast that they make the coalescent collapse. A counterpart to Theorem, with τ n in replaced by its expectation, was obtained by Herriger and Möhle []. Our second result is a central limit theorem. Here we confine ourselves to coalescents with µ <. Then the function fy := [,] p ey e y Λdp p 2, y R 2 is everywhere finite. Also f is a positive, monotone decreasing, continuous function with the property fy for y. Let where we choose κ such that b n := logn κ dy µ fy, fy µ 2 for all y κ. Theorem 2. Assume that µ < and moreover σ 2 2 Λdp := log p p 2 <. Then as n. [,] τ n b n logn d N,σ 2 /µ 3 3 2

3 Under the additional condition [,] log p Λdp p <. 4 the CLT 3 has been obtained by Gnedin, Iksanov and Marynych [2], with b n replaced by logn/µ. Their condition 9 is equivalent to the above condition 4, see Remark 3 in [4]. Thus the question arises, whether the simplified centering by log n/µ is always feasible. The next proposition shows that this can be done under a condition that is weaker than 4, but not in any case. Proposition 3. Let c <. Then b n = logn µ + 2c µ 2 logn+o logn 5 as n, if and only if log Λdp r c 6 [,r] p as r. Example. We consider for γ R the finite measures Λdp = +log p γdp, p. For γ = this gives the Bolthausen-Sznitman coalescent. For γ > it leads to coalescents with µ,σ 2 <. Note that 4 is satisfied iff γ > 2, and 6 is fulfilled iff γ > 3/2. Thus within the range < γ 3/2 one has to come back to the constants b n in the central limit theorem. The law of large numbers from Theorem holds for all γ >. For the regime γ, Theorem just tells us that τ n = o P logn. For γ =, the Bolthausen-Sznitman coalescent, it is known that τ n is already down to the order loglogn [3]. For γ <, applying Schweinsberg s criterion [8], it can be shown that the coalescents come down from infinity. There remains the gap < γ. It is tempting to conjecture that τ n is of order logn γ for < γ <. If equation 6 is violated then the subsequent approximation to b n may be practical. Starting from the identity we obtain the expansion b n = logn µ + µ 2 µ fy = µ + fy µ 2 + f2 y µ fk y µ k+ + f k+ y µ k+ µ fy logn fydy + + µ k+ logn f k ydy +O logn f k+ ydy. 3

4 Let us now explain the method of proving Theorems and 2. We are mainly dealing with Λ-coalescents having a dust component. Shortly speaking these are the coalescents for which the rate, at which a single lineage merges with some others from the sample, stays bounded as the sample size tends to infinity. As is well-known this property is characterized by the condition Λdp <. 7 p [,] An established tool for the analysis of a Λ-coalescent with dust is the subordinator S = S t t, which is used to approximate the logarithm of its block-counting process N n = N n t t see e.g. Pitman [6], Möhle [5], and theabove mentioned paperby Gnedinet al [2]. We will recall this subordinator in Sec. 3. Indeed, analogues of Theorems and 2 are well-known for first-passage times of subordinators with finite first resp. second moment. However, this approximation neglects the subtlety that a coalescent of b lineages results in a downward jump of size b and not b for the process N n. This effect becomes significant when many small jumps accumulate over time, as it happens close to the dustless case and as it becomes visible in Proposition 3 and in the above example. Then the appropriate approximation is provided by a drifted subordinator Y n = Y n t t, given by the SDE Y n t = logn S t + t fy n sds, t, with initial value Y n = logn. The drift compensates the just mentioned difference between b and b. In Kersting et al [4] it is shown that sup Yn t logn n t = OP t<τ n as n, that is, these random variables are bounded in probability. In Sec. 3 we suitably strengthen this result. In Sec. 2 we provide the required limit theorems for passage times for a more general class of drifted subordinators. The above results are then proved in Sec. 4. It turns out that the regime considered by Gnedin et al [2] is the one in which the random variables τ n fy n sds are bounded in probability uniformly in n. This can be seen to be equivalent to the requirement fydy <, which likewise is equivalent to 4 see the proof of Corollary 2 in [4]. Under this assumption Gnedin et al [2] proved their central limit theorem also with non-normal stable or Mittag-Leffler limiting distributions of τ n. A similar generalization of Theorem 2 is feasible in the general dust case, without the requirement 4. 2 Limit theorems for a drifted subordinator Let S = S t t be a pure jump subordinator with Lévy measure λ on,. Recall that this requires y λdy <. With regard to the mentioned properties of the function in 2, let f : R R be an arbitrary positive, non-increasing, continuous function with lim fy =. y 4

5 Let the process Y z = Y z t t denote the unique solution of the SDE Y z t = z S t + t fys z ds 8 with initial value z >. We will investigate the asymptotic behaviour of its passage times across x R, T z x := inf{t : Y z t < x}, in the limit z. The first result provides a law of large numbers. Denote µ := yλdy. 9, Proposition 4. Assume that µ <. Then for any x R in probability as z. Proof. Let z > x. Then {T z x z Tz x µ t} = {Y z s x for all s t} = {S s z x+ s } fyu z du for all s t. By positivity of the function f it follows PTx z t PS t z x, thus for any ε > P Tx z ε z P S µ εz/µ z x. Now µ = E[S ], thus by the law of large numbers S t t µ a.s., hence the right-hand term in converges to for z and also P Tx z ε z. µ On the other hand, { {Tx z t} = {Ys z x for all s t} = Ys z x for all s t, S t z x+ t } fys z ds. Monotonicity of f implies PTx z t P S t z x+tfx}. Therefore, since fx as x, P Tx z +ε z P S µ +εz/µ z x++ε z µ fx P S +εz/µ z+ε/2 x 5

6 if only x is sufficiently large. Now the right-hand term converges to, thus it follows that P Tx z +ε z. µ Note that we proved this result only for x sufficiently large, depending on ε. However, this restriction may be skipped, since for fixed x < x 2 the random variables Tx z Tx z 2 are bounded in probability uniformly in z. Thus altogether we have for any x P ε z µ Tz x < +ε z µ as z, which since ε > was arbitrary is our assertion. Now we turn to a central limit theorem for passage times of the processes Y z. Let the function β z, z κ, be given by z dy β z := µ fy, where we choose κ so large that Proposition 5. Suppose that σ 2 := κ supfy µ y κ 2., y 2 λdy <. 2 Then T z x β z z d N,σ 2 /µ 3 as z. Proof. i Note again that for x < x 2 the random variables T z x T z x 2 are bounded in probability uniformly in z. Thus it suffices to prove our theorem for all x x for some x R. Therefore we may change fx for all x < x. We do it in such a way that fx µ/2 for all x R, without touching the other properties of f. Thus we assume from now that fy µ 2 for all y R 3 and set κ = in. Consequently, z µ β z 2z µ, z >. 4 For any z > we define the function ρ z t = ρ z t, t β z, such that β ρ z t = β z t, t β z, in particular ρ z = z and ρ z β z =. This means that ρ z arises by first inverting the function β restricted to the interval [,z], and then reversing the time parameter on its domain [,β z ]. By differentiation we obtain ρ z t = fρz t µ, 6

7 consequently ρ t µ/2 and ρ z t = z µt+ t fρ z sds. ii A glimpse on 8 suggests that ρ z will make a good approximation for the process Y z. In order to estimate their difference observe that For given t > define u t = Y z t ρ z t = S t µt+ t fy z s fρz s ds. { sup{s < t : Y z s ρ z s} on the event Y z t > ρ z t sup{s < t : Y z s ρz s } on the event Y z t < ρ z t and u t := t on the event Yt z = ρ z t. We have u t t, since Y z = z = ρz. Because f is a decreasing function, the event Yt z > ρ z t implies that On the event Y z t < ρ z t Y z t ρ z t Y z t ρ z t t u t fy z s fρz s ds Y z u t ρz u t = S t µt+s ut µu t. there is an analogous estimate from below, altogether Y z t ρ z t 2M t with M t := sup S u µu. u t Consequently, Ys z ρ z s 2M s ρ z s 2M t for s t and by means of the monotonicity of f t fy z s ds t fρ z sds t fρ z s 2M t ds An analoguous estimate is valid from below and we obtain t fy z s ds t t fρ z sds 2M t fρ z t 2M t. fρ z sds 2M t fρ z t 2M t. 5 At this point we recall that under the above assumptions on the subordinator S by Donsker s invariance principle we have M t = O P t as t. iii Now we derive some upper estimates of probabilities. Given a,x R, we have for any c > PTx z β z +a z = PYt z x for all t β z +a z βz+a z = P S βz+a z z x+ fys z ds, Yt z x for all t β z +a z P S βz+a z z x+fxc+ a z + βz c z fys z ds 7

8 We now bring 5 into play. From the definition of ρ z we have, writing βy = β y, that thus because of 4 Then on the event M βz c z/8µ we have βρ z β z c z = c z, ρ z β z c z c z 2µ. ρ z β z z 2M βz c z c z 2µ c z 4µ = c z 4µ. Consequently, by means of 5 and since β z 2z/µ PTx z β z +a z P M 2z/µ > c z 8µ +P S βz+a z z x+fxc+ a z + Moreover, by definition of ρ z, βz βz z + fρ z sds = ρ z β z +µβ z = µβ z. fρ z s ds+ c z c z 4µ f. 6 4µ Therefore, if we fix ε >, let c be so large that the first right-hand probability in 6 is smaller than ε, then choose z so large that c/4µf c z 4µ ε, and also choose x > and so large that cfxc+ a ε, then we end up with PTx z β z +a z ε+p S βz+a z µβ z +2ε z. Also by the law of large numbers in probability. Therefore Moreover µβ z z, hence S βz+a z S βz µa z PT z x β z +a z 2ε+P S βz µβ z + µa+3ε z. PT z x β z +a z 2ε+P S βz µβ z + µa+4εµ /2 β z for large z. Now from assumption 2 and the central limit theorem there follows S t µt σ 2 t d L, where L denotes a standard normal random variable. Thus lim supptx z β z +a z 2ε+PL µa+4εµ /2 σ. z 8

9 Note that the choice of x dependson ε in our proof. However, since again the differences Tx z Tx z 2 are bounded in probability uniformly in z, this estimate generalizes to all x. Now letting ε we obtain T z lim supp x β z a PL µ 3/2 σ a. z z This is the first part of our claim. iv For the lower estimates we first introduce the random variable R z,x := sup{t : Y z t x} inf{t : Y z t < x} which is the length of the time interval where Y z t x is changing from positive sign to ultimately negative signnotethat thepathsof Y z arenot monotone. We claim that theserandomvariables are bounded in probability, uniformly in z and x. Indeed, with we have for t > η = η z,x because of Y z η Y z t = Y z η S t S η + Thus R z,x is bounded from above by η z,x := inf{t : Y z t < x} x and 3 t η fy z s ds x S t S η + µ 2 t η. R z,x := sup{u : S ηz,x+u S ηz,x µu/2 } These random variables are a.s. finite. Moreover, they are identically distributed, since η z,x are stopping times. This proves that the R z,x are uniformly bounded in probability. Now for the lower bounds we have for a,b R PT z x β z +a z PY z t x for all t β z +a z, R z,x b = PY z t x for all β z +a z b t β z +a z, R z,x b. For these t we have Y z t = z S t + t βz+a z b fys z ds z S β z+a z + therefore PTx z β z +a z P S βz+a z b βz+a z z x+ P S βz c z βz+a z z x+ fy z s ds, fys z ds, R z,x b fy z s ds PR z,x > b for c sufficiently large. We now bring, as in part iii, 5 into play. Proceeding analogously we obtain instead of 6 the estimate PTx z β z +a z PR z,x > b P M 2z/µ > c z 8µ βz c z +P S βz+a z z x+ fρ z sds c z c z 4µ f. 4µ 9

10 Also, since ρ z β z = and ρ z t µ/2, βz β z c z c z fρ z s ds fµs/2ds = o z. Hence, for given ε > and z sufficiently large PTx z β z +a z PR z,x > b P M 2z/µ > c z 8µ βz +P S βz+a z z ε z + fρ z s ds c z c z 4µ f. 4µ Returning to the arguments of part iii we choose b, c and then z so large that we arrive at PTx z β z +a z 2ε+P S βz+a z µβ z 2ε z and further at lim inf z PTz x β z +a z 3ε+PL µa 3εµ /2 σ. The limit ε leads to the desired lower estimate. 3 Approximating the block counting process In this section we derive a strengthening of a result in Kersting, Schweinsberg and Wakolbinger [4] on the approximation to the logarithm of the block counting processes in the dust case. To this end, let us quickly recall the Poisson point process construction of the Λ-coalescent given in [4], which is a slight variation of the construction provided by Pitman in [6]. This construction requires Λ{} =, which is fulfilled for coalescents with dust. Consider a Poisson point process Ψ on,,] [,] n with intensity dt p 2 Λdp du du n, and let Π n = {{},...,{n}} be the partition of the integers,...,n into singletons. Suppose t,p,u,...,u n is a point of Ψ, and Π n t consists of the blocks B,...,B b, ranked in order by their smallest element. Then Π n t is obtained from Π n t by merging together all of the blocks B i for which u i p into a single block. These are the only times that mergers occur. This construction is well-defined because almost surely for any fixed t <, there are only finitely many points t,p,u,...,u n of Ψ for which t t and at least two of u,...,u n are less than or equal to p. The resulting process Π n = Π n t,t is the Λ-coalescent. When t,p,u,...,u n is a point of Ψ, we say that a p-merger occurs at time t. Condition 7 allows us to approximate the number of blocks in the Λ-coalescent by a subordinator. Let φ :,,] [,] n,, ] be the function defined by φt,p,u,...,u n = t, log p. Now φψ is a Poisson point process, and we can define a pure jump subordinator St,t having the property that S = and, if t,x is a point of φψ, then St = St +x. With λ the Lévy measure of S, the formulas 9 and 2 now read µ = [,] log Λdp p p 2 and σ 2 = [,] 2 Λdp log p p 2.

11 ThissubordinatorfirstappearedintheworkofPitman[6]andwasusedtoapproximatetheblockcounting process by Gnedin et al. [2] and Möhle [5]; the benefits of a refined approximation by a drifted subordinator were discovered in [4]. We recall that the drift appears because a merging of b out of N n t lines results in a decrease by b and not by b lines, see equation 23 in [4] for an explanation of the form of the drift. The next result provides a refinement of Theorem in [4]. Proposition 6. Let [,] Λdp p <, let f be as in 2, and let Y n be the solution of 8 with z := logn. Then for any ε > there is an l < such that P sup t<τ n logn n t Y n t l, Y n τ n < l ε. Proof. From [4] we know that for given ε > there is an r < such that P sup t<τ n logn n t Y n t r ε/2. Now we consider the size n of the last jump. Letting u i,p i, i, be the points of the underlying Poisson point process with intensity measure dtλdp/p 2, the associated subordinator S has jumps of size v i = log p i at times t i. Thus for any c > we have { n logn n τ n c} = {τ n = t i and log p i logn n t i c for some i } { e c } = τ n = t i and p i N n t i for some i Given N n t this event appears at time t with rate ν n,t = p Nnt Λdp p 2. [, e c /N nt ] Using the inequalities p b = p b e pb / pb we get ν n,t e pnnt 2 e 2 Λdp pn n t Λdp. [, e c /N nt ] [, e c /N nt ] It follows [ E ] ν n,t dt E [ [,] e 2 ] pn n t I {N nt e c / p }dtλdp Lemma 4 of [4] yields the estimate with some c >, hence [ ] E N n t I {N nt e c / p }dt c e c / p p c e c [ E ] ν n,t dt c e 2 c Λ[,].

12 Therefore for c sufficiently large which implies Altogether we obtain [ E ] ν n,t dt ε/2, P n logn n τ n c = exp E [ ] ν n,t dt ε/2. P sup t<τ n logn n t Y n t r, n > logn n τ n c ε. The event in the previous formula implies Y n τ n = Y n τ n n < logn n τ n +r logn n τ n c = r+c, and the claim of the theorem follows with l = r +c. 4 Proof of the main results Proof of Theorem. Let us first assume that µ <. Then we have a coalescent with dust, and we may apply Proposition 6. Fix η >. Note that on the event that Y n τ n < l the event τ n < ηlogn/µ implies the inequality T logn l < ηlogn/µ. Thus in view of Proposition 6 there exists for any ε > an l such that Pτ n < ηlogn/µ PT logn l < ηlogn/µ+ε. Proposition 4 implies that the right-hand probability converges to as n. Letting ε we obtain lim n Pτ n < ηlogn/µ =. Also on the event sup t<τn logn n t Y n t l, the event τ n > + ηlogn/µ implies Y n t l for all t +ηlogn/µ, and consequently Pτ n > +ηlogn/µ PT logn l > +ηlogn/µ+ε. Again the right-hand probability converges to zero in view of Proposition 4, and we obtain lim Pτ n > +ηlogn/µ =. n Altogether our claim follows in the case µ <. Now assume µ =. If Λ{} >, then the coalescent comes down from infinity and τ n stays bounded in probability. The same is true if Λ{} >, thus we may assume that Λ{,} =. For given ε > define the measure Λ ε by Λ ε B := ΛB [ε, ε]. Obviously µ ε := log Λ ε dp p p 2 <. 2

13 Thus for the absorption times τ ε n of the Λ ε -coalescent we have τ ε n logn µ ε in probability as n. Now we may couple the Λ ε -coalescent in an obvious manner to the Λ-coalescent in such a way that N n t N ε nt a.s. for all t, in particular τ n τ ε n. Hence it follows that Pτ n /logn > 2/µ ε. Because of Λ{,} = we have µ ε µ = with ε, consequently for all η >. This is our claim. Pτ n /logn > η Proof of Theorem 2. Because of the condition µ < we again may apply Proposition 6. We follow the same line as in the previous proof: For ε > there exists an l such that for all a R and Pτ n < b n +a n PT logn l Pτ n > b n +a n PT logn l Now apply Proposition 5 and let ε. < b n +a n+ε > b n +a n+ε Proof of Proposition 3. i Let us first assume 6. Because of p /r minp/r, for < r < we have for α > f log r α Λdp Λdp r α Λdp r +r p r α p 2 +r α Λdp. 7 p p Also, because of p /r e p/r e p/r p/r, it follows for β > that f log r β r e r β Together with 6 these two estimates yield for α < < β Λdp p. 8 cβ /2 liminf r f log r log r limsup f log r log r cα /2. r Letting α,β we arrive at fy = c+o/ y as y and consequently as n. Now, because of logn fydy = c+o2 logn µ fy = µ + fy µµ fy 3

14 and fy = o as y, we have z dy µ fy = z µ + +o z µ 2 fydy +O 9 as z, and consequently, as claimed, κ b n = logn µ + 2c+o µ 2 logn. ii Now suppose that 5 is satisfied. Then in view of 9 with z = logn it follows that as n, or equivalently logn z fydy = 2c+o logn fydy = 2c+o z for z. This implies that fz = c + o/ z as z. For c = this claim follows because f is decreasing, which entails For c > we use the estimate zfz z fydy = o z. η z +ηz z fydy zfz z η fydy z ηz with η >. Taking the limit z and then η yields fz = c+o/ z. Now, similar as in part i we get from 7 and 8 c α liminf r With α,β we arrive at 6. log Λdp r [,r] p limsup log Λdp r c β. r [,r] p Acknowledgement. It is our pleasure to dedicate this work to Peter Jagers. References [] Herriger, Ph. and Möhle, M. 22. Conditions for excheangable coalescents to come down from infinity. ALEA 9, [2] Gnedin, A., Iksanov, A., and Marynych, A. 2. On Λ-coalescents with dust component, J. Appl. Probab. 48, [3] Goldschmidt, Ch. and Martin, J. 25. Random Recursive Trees and the Bolthausen- Sznitman Coalesent. Electron. J. Probab.,

15 [4] Kersting, G., Schweinsberg, J., and Wakolbinger, A. 27. The size of the last merger and time-reversal in Λ-coalescents. To appear in AIHP, electronic version available at [5] Möhle, M On hitting probabilities of beta coalescents and absorption times of coalescents that come down from infinity. ALEA, [6] Pitman, J Coalescents with multiple collisions. Ann. Probab. 27, [7] Sagitov, S The general coalescent with asynchronous mergers of ancestral lines. J. Appl. Probab. 36, [8] Schweinsberg, J. 2. A necessary and sufficient condition for the Λ-coalescent to come down from infinity. Electron. Comm. Probab. 5,. 5