We extend the Poincaré–Borel lemma to a weak approximation of a Brownian motion via simple functionals of uniform distributions on n-spheres in the Skorokhod space D([0,1]). This approach is used to simplify the proof of the self-normalized Donsker theorem in Csörgő et al. (2003). Some notes on spheres with respect to ℓp-norms are given.
Let Sn−1(d)={x∈Rn:‖x‖=d} be the (n−1)-sphere with radius d, where ‖·‖ denotes the Euclidean norm. The uniform measure on the unit sphere Sn−1:=Sn−1(1) can be characterized as
μS,n=d(X1,…,Xn)‖(X1,…,Xn)‖,
where (X1,…,Xn) is a standard n-dimensional normal random variable.
The celebrated Poincaré–Borel lemma is the classical result on the approximation of a Gaussian distribution by projections of the uniform measure on Sn−1(n) as n tends to infinity: Let n≥m and πn,m:Rn→Rm be the natural projection. The uniform measure on the sphere Sn−1(n) is given by nμS,n. Then, for every fixed m∈N,
nμS,n∘πn,m−1
converges in distribution to a standard m-dimensional normal distribution as n tends to infinity, cf. [11, Proposition 6.1]. Following the historical notes in [6, Section 6] on the earliest reference to this result by Émile Borel, we acquire the usual practice to speak about the Poincaré–Borel lemma.
Among other fields, this convergence stimulated the development of the infinite-dimensional functional analysis (cf. [12]) as well as the concentration of measure theory (cf. [10, Section 1.1]).
In particular, it inspired to consider connections of the Wiener measure and the uniform measure on an infinite-dimensional sphere [21]. Such a Donsker-type result is firstly proved in [4] by nonstandard methods. For the illustration, we make use of the notations in [7], where this result is used for statistical analysis of measures on high-dimensional unit spheres. Define the functional
Qn,2:Sn−1→C([0,1]),(x1,…,xn)↦(Qn,2(t))t∈[0,1],
such that
Qn,2(k/n):=∑i=1kxi‖(x1,…,xn)‖,
for k∈{0,…,n} and is linearly interpolated elsewhere. Then [4, Theorem 2.4] gives that the sequence of processes
μS,n∘Qn,2−1
converges weakly to a Brownian motion W:=(Wt)t∈[0,1] in the space of continuous functions C([0,1]) as n tends to infinity. The first proof without nonstandard methods in C([0,1]) and in the Skorohod space D([0,1]) is given in [17].
In this note, we present a very simple proof of the càdlàg version of this Poincaré–Borel lemma for Brownian motion. This is the content of Section 2.
Some remarks on such Donsker-type convergence results on spheres with respect to ℓp-norms are collected in Section 3.
In fact, our simple approach can be used to simplify the proof of the main result in [3] as well. This is presented in Section 4.
Poincaré–Borel lemma for Brownian motion
Suppose X1,X2,… is a sequence of i.i.d. standard normal random variables. Then (X1,…,Xn) has a standard n-dimensional normal distribution. We define the processes with càdlàg paths
Zn=(Ztn:=∑i=1⌊nt⌋Xi‖(X1,…,Xn)‖)t∈[0,1].
Thus, Zn is equivalent to μS,n∘Q‾n,2−1 for the functional
Q‾n,2:Sn−1→D([0,1]),(x1,…,xn)↦(Q‾n,2(t)=∑i=1⌊nt⌋xi‖(x1,…,xn)‖)t∈[0,1],
and therefore it is a relatively simple computation from the uniform distribution on the n-sphere. Then the following extension of the Poincaré–Borel lemma is true:
The sequence(Zn)n∈Nconverges weakly in the Skorokhod spaceD([0,1])to a standard Brownian motion W as n tends to infinity.
As the distribution of the random vector in (1) is exactly the uniform measure μS,n, the proof of the convergence of finite-dimensional distributions is in line with the classical Poincaré–Borel lemma: by the law of large numbers, 1n∑i=1nXi2→1 in probability. Hence, by the continuous mapping theorem, n/(‖(X1,…,Xn)‖)→1 in probability, and, by Donsker’s theorem and Slutsky’s theorem, we conclude the convergence of finite-dimensional distributions.
For the tightness we consider the increments of the process Zn and make use of a standard criterion. For all s≤t in [0,1], we denote
(Ztn−Zsn)2=∑⌊ns⌋<i≤⌊nt⌋Xi2∑i≤nXi2+∑⌊ns⌋<i≠j≤⌊nt⌋XiXj∑i≤nXi2=:I1t,s+I2t,s.
Due to the symmetry of the standard n-dimensional normally distributed vector (X1,…,Xn), for all pairwise different i,j,k,l, we observe
E[XiXjXkXl(∑i≤nXi2)2]=E[Xi2XjXk(∑i≤nXi2)2]=0.
Let s≤u≤t in [0,1]. Thus via (3), we conclude
E[I1t,uI2u,s]=0,E[I2t,uI1u,s]=0,E[I2t,uI2u,s]=0,
and therefore
E[(Ztn−Zun)2(Zun−Zsn)2]=E[I1t,uI1u,s].
We denote for shorthand m1:=⌊nt⌋−⌊nu⌋, m2:=⌊nu⌋−⌊ns⌋ and m3:=n−(⌊nt⌋−⌊ns⌋). Then we observe
I1t,uI1u,s=χm12χm22(χm12+χm22+χm32)2=12((χm12+χm22)2−(χm12)2−(χm22)2)(χm12+χm22+χm32)2,
for pairwise independent chi-squared random variables χm2 with m degrees of freedom. We recall that χm2χm2+χk2 is Beta(m/2,k/2)-distributed with
E[(χm2χm2+χk2)2]=(m+2m+k+2)(mm+k).
Hence a computation via (4) yields
E[I1t,uI1u,s]=m1m2(m1+m2+m3+2)(m1+m2+m3)≤(m1m1+m2+m3)(m2m1+m2+m3),
and therefore
E[(Ztn−Zun)2(Zun−Zsn)2]≤(⌊nt⌋−⌊nu⌋n)(⌊nu⌋−⌊ns⌋n)≤(⌊nt⌋−⌊ns⌋n)2.
Thus the well-known criterion [1, Theorem 15.6] (cp. Remark 1 in [15]) implies the tightness of Zn. □
(i) The heuristic connection of the Wiener measure and the uniform measure on an infinite-dimensional sphere goes back to Norbert Wiener’s study of the differential space, [21]. The first informal presentation of Theorem 1 and further historical notes can be found in [12]. The first rigorous proof is given in Section 2 of [4]. However, the authors make use of nonstandard analysis and the functional Qn,2. To the best of our knowledge, the first proof of Theorem 1 is [17]. In contrast, our proof is based on the pretty decoupling in the tightness argument. Moreover, this approach is extended in Section 4 to a simpler proof of Theorem 1 in [3].
(ii) According to the historical comments in [20, Section 2.2], the Poincaré–Borel lemma could be also attributed to Maxwell and Mehler.
In this section, we consider uniform measures on ℓpn-spheres and prove that the limit in Theorem 1 is the only case such that a simple Q‾-type pathwise functional leads to a nontrivial limit (Theorem 5).
Furthermore, we present random variables living on ℓpn-spheres, with a similar characterization for a fractional Brownian motion (Theorem 6).
Concerning the ℓpn norm ‖x‖p=(∑i=1n|xi|p)1/p for p∈[1,∞) and defining the ℓpn unit sphere
Spn−1:={x∈Rn:‖x‖p=1},
the uniform measure μS,n,p on Spn−1 is characterized similarly to the uniform measure on the Euclidean unit sphere by independent results in [18, Lemma 1] and [16, Lemma 3.1]:
SupposeX,X1,X2,…is a sequence of i.i.d. random variables with densityf(x)=exp(−|x|p/p)2p1/pΓ(1+1/p).ThenμS,n,p=d(X1,…,Xn)‖(X1,…,Xn)‖p.
(i) We notice that the uniform measure on the ℓpn-sphere equals the surface measure only in the cases p∈{1,2,∞}, see e.g. [16, Section 3] or the interesting study of the total variation distance of these measures for p≥1 in [14].
(ii) In particular, we have a counterpart of the classical Poincaré–Borel lemma for finite-dimensional distributions: For every fixed m∈N,
n1/pμS,n,p∘πn,m−1
converges in distribution to the random vector (X1,…,Xm) as n tends to infinity. This follows immediately from E[|X|p]=1 and the law of large numbers, cf. [11, Proposition 6.1] or the finite-dimensional convergence in Theorem 1.
Similarly to the characterization of the central limit theorem, cp. [9, Theo- rem 4.23], but in contrast to the convergence of the projection on a finite number of coordinates in Remark 4, we have a uniqueness result for the processes constructed according to the Q‾-type pathwise functionals.
In the following we denote the convergence in distribution by →d and the almost sure convergence by →a.s..
Supposep≥1and denoteQ‾n,p:(x1,…,xn)↦(∑i=1⌊nt⌋xi‖(x1,…,xn)‖p)t∈[0,1].Then, in the Skorokhod spaceD([0,1]), as n tends to infinity:μS,n,p∘Q‾n,p−1→a.s.0,p<2,→dW,p=2,is divergent,p>2.2.\end{array}\right.\]]]>
The strong law of large numbers [9, Theorem 4.23] implies that n1/p/‖(X1,…,Xn)‖p→1 almost surely for all p≥1. Moreover, for p<2, it gives as well that 1n1/p∑i=1⌊nt⌋Xi→0 almost surely for all t∈[0,1]. Thanks to Proposition 3, we have
μS,n,p∘Q‾n,p−1=dn1/p‖(X1,…,Xn)‖p(n−1/p∑i=1⌊n·⌋Xi).
Thus we conclude via n−1/p=n−1/2n(p−2)/2p, Donsker’s theorem and Slutsky’s theorem. □
However, the ℓpn spheres can be involved in another convergence result. The fractional Brownian motion BH=(BtH)t≥0 with Hurst parameter H∈(0,1) is a centered Gaussian process with the covariance E[BtHBsH]=12(t2H+s2H−|t−s|2H). We refer to [13] for further information on this generalization of the Brownian motion beyond semimartingales. In particular, there is the following random walk approximation ([19, Theorem 2.1] or [13, Lemma 1.15.9]): Let {Xi}i≥1 be a stationary Gaussian sequence with E[Xi]=0 and correlations
∑i,j=1nE[XiXj]∼n2HL(n),
as n tends to infinity for a slowly varying function L. Then 1n2HL(n)∑i=1⌊nt⌋Xi converges weakly in the Skorohod space D([0,1]) towards a fractional Brownian motion with Hurst parameter H. For simplification let Xi=BiH−Bi−1H, i∈N, be the correlated increments of the fractional Brownian motion BH. The stationarity and the ergodic theorem imply, for p>00$]]> and the constant cH:=E[|B1H|1/H], that
(‖(X1,…,Xn)‖p/nH)p=n−Hp∑i=1n|Xi|p→a.s.0,p>1/H,cHp=1/H,+∞,p<1/H,1/H,\\{} c_{H}\hspace{1em}& p=1/H,\\{} +\infty ,\hspace{1em}& p<1/H,\end{array}\right.\]]]>
(see e.g. [13, Eq. (1.18.3)]). With this at hand, we obtain a similar uniqueness result:
LetXi=BiH−Bi−1H,i∈N, be the increments of a fractional Brownian motionBH. Then, in the Skorokhod spaceD([0,1]), as n tends to infinity:Q‾n,p(X1,…,Xn)=(∑i=1⌊nt⌋Xi‖(X1,…,Xn)‖p)t∈[0,1]→a.s.0,p<1/H,→dBH/cHH,p=1/H,is divergent,p>1/H.1/H.\end{array}\right.\]]]>
Taqqu’s limit theorem implies, for all H∈(0,1),
(n−H∑i=1⌊nt⌋Xi)t∈[0,1]→dBH
in the Skorokhod space D([0,1]). Then, thanks to (5), we conclude as in Theorem 5. □
Due to the different correlations between the random variables Xi in Theorem 6, there is no symmetric and trivial sequence of measures μˆS,n,p on the ℓpn-spheres and some simple Q‾n,p-type pathwise functionals, which represent the distributions of Q‾n,p(X1,…,Xn). However, it would be interesting, whether some uniform or surface measures on geometric objects in combination with simple Q‾p-type pathwise functionals allow similar Donsker-type theorems for fractional Brownian motion or other Gaussian processes?
The self-normalized Donsker theorem
Suppose X,X1,X2,… is a sequence of i.i.d. nondegenerate random variables and we denote for all n∈N,
Sn:=∑i=1nXi,Vn2:=∑i=1nXi2.
Limit theorems for self-normalized sums Sn/Vn play an important role in statistics, see e.g. [8], and have been extensively studied during the last decades, cf. the monograph on self-normalizes processes [5].
In [3], the following invariance principle for self-normalized sum processes is established.
(Theorem 1 in [<xref ref-type="bibr" rid="j_vmsta82_ref_003">3</xref>]).
Assume the notations above and denoteZtn:=S⌊nt⌋/Vn.Then the following assertions, with n tending to infinity, are equivalent:
E[X]=0and X is in the domain of attraction of the normal law (i.e. there exists a sequence(bn)n≥1withSn/bn→dN(0,1)).
For allt0∈(0,1],Zt0n→dN(0,t0).
(Ztn)t∈[0,1]converges weakly to(Wt)t∈[0,1]on(D([0,1]),ρ), where ρ denotes the uniform topology.
On an appropriate joint probability space, the following is valid:supt∈[0,1]|Ztn−W(nt)/n|=oP(1).
The equivalence of (a) and (b) is the celebrated result [8, Theorem 3]. Since the implications (d)⇒(c)⇒(b) are trivial, the proof in [3] is completed by showing (a)⇒(d).
Thanks to a tightness argument as in the proof of Theorem 1, we obtain a simpler alternative for the proof.
As stated in the remark, we already know that (d)⇒(c)⇒(b)⇔(a). We denote
(Ztn)t∈[0,1] converges weakly to (Wt)t∈[0,1] on the Skorokhod space D([0,1]).
By the continuity of the paths of the Brownian motion and [1, Section 18], we obtain the equivalence (c)⇔(c0). We denote by d0 the Skorokhod metric on D([0,1]) which makes it a Polish space. The Skorokhod–Dudley Theorem [9, Theorem 4.30] and (c0) imply
d0((Ztn)t∈[0,1],(Wt)t∈[0,1])→0,
almost surely on an appropriate probability space. Since the uniform topology is finer than the Skorokhod topology ([1, Section 18]), we conclude assertion (d). Thus it remains to prove (a)⇒(c0). Firstly we consider finite-dimensional distributions. Due to [8, Lemma 3.2], the sequence (bn)n∈N with Sn/bn→dN(0,1) fulfills Vn/bn→1 in probability and bn=nL(n) for some slowly varying at infinity function L. The continuous mapping theorem implies bn/Vn→1 in probability. Take arbitrary N∈N, a1,…,aN∈R and t1,…,tN∈[0,1]. Without loss of generality, we assume t1<⋯<tN and denote t0:=0 and tN+1:=1. Then, by the independence of the random variables S⌊nti⌋−S⌊nti−1⌋, i=1,…,N+1, for every fixed n∈N, Lévy’s continuity theorem and the normality of the random vector (Y1,…,YN+1), we obtain
(S⌊nt1⌋−S⌊nt0⌋(⌊nt1⌋),…,S⌊ntN+1⌋−S⌊ntN⌋(⌊ntN+1⌋−⌊ntN⌋))→d(Y1,…,YN+1),
as n tends to infinity. As the sequence (bn)n∈N is regularly varying with exponent 1/2, it is easily seen that
b⌊nti⌋−⌊nti−1⌋bn→ti−ti−1.
Via the continuous mapping theorem, we conclude
∑iaiS⌊nti⌋bn=∑i=1N+1(∑j≤iaj)(b⌊nti⌋−⌊nti−1⌋)bn(S⌊nti⌋−S⌊nti−1⌋b⌊nti⌋−⌊nti−1⌋)→d∑i=1N+1(∑j≤iaj)ti−ti−1Yi=d∑i=1N+1aiWti.
Slutsky’s theorem implies
∑i=1N+1aiZtin=(bnVn)(∑iaiS⌊nti⌋bn)→d∑i=1N+1aiWti,
what means the convergence of finite-dimensional distributions.
The tightness follows again by the criterion [1, Theorem 15.6]. By the identical distribution, for all m≤n, we have
E[(∑i≤mXi2∑i≤nXi2)2]=E[mX14(∑i≤nXi2)2]+E[m(m−1)X12X22(∑i≤nXi2)2].
Thanks to the value 1 on the left hand side in (6) for m=n, we conclude
0≤E[X12X22(∑i≤nXi2)2]≤1n(n−1).
In contrast to (3), for possibly nonsymmetric random variables, the Cauchy–Schwarz inequality and [8, (3.10)] yields a constant cX<∞ such that for every r∈{2,3,4},
maxi,j,k,l≤n|{i,j,k,l}|=rE[|XiXjXkXl|(∑i≤nXi2)2]≤cXn−r.
Applying the estimates in (7) on the terms in (2) gives that
maxi,j∈{1,2}E[Iit,uIju,s]≤cX(⌊nt⌋−⌊ns⌋n)2.
Hence, we obtain
E[(Ztn−Zun)2(Zun−Zsn)2]=E[(I1t,u+I2t,u)(I1u,s+I2u,s)]≤4cX(⌊nt⌋−⌊ns⌋n)2,
and the proof concludes as in Theorem 1. □
(i) By the same reasoning, we obtain Theorem 5 for the sequence of i.i.d. variables X,X1,X2,… such that Theorem 8(a) is fulfilled.
(ii) In [2], a similar counterpart of Theorem 8 for α-stable Lévy processes is established. An interesting question would be on a uniqueness result similar to Theorem 5.
Acknowledgments
The author thanks the anonymous referees for careful reading and helpful comments.
ReferencesBillingsley, P.: Convergence of Probability Measures. Wiley, New York (1968). MR0233396Csörgő, M., Hu, Z.: Weak convergence of self-normalized partial sums processes. In: Asymptotic Laws and Methods in Stochastics. Fields Inst. Commun., 76, pp. 3–15. Springer (2015). MR3409822. doi:10.1007/978-1-4939-3076-0_1Csörgő, M., Szyszkowicz, B., Wang, Q.: Donsker’s theorem for self-normalized partial sums processes. Ann. Probab.31 (3), 1228–1240 (2003). MR1988470. doi:10.1214/aop/1055425777Cutland, N., Ng, S.-A.: The wiener sphere and wiener measure. Ann. Probab.21 (1), 1–13 (1993). MR1207212de la Peña, V.H., Lai, T.L., Shao, Q.-M.: Self-normalized Processes. Springer, Berlin (2009). MR2488094. doi:10.1007/978-3-540-85636-8Diaconis, P., Freedman, D.: A dozen de finetti-style results in search of a theory. Ann. Inst. Henri Poincaré Probab. Stat.23 (2), 397–423 (1987). MR0898502Dryden, I.L.: Statistical analysis on high-dimensional spheres and shape spaces. Ann. Stat.33 (4), 1643–1665 (2005). MR2166558. doi:10.1214/009053605000000264Giné, E., Götze, F., Mason, D.M.: When is the student t-statistic asymptotically standard normal?Ann. Probab.25, 1514–1531 (1997). MR1457629. doi:10.1214/aop/1024404523Kallenberg, O.: Foundations of Modern Probability. Second Edition. Probability and Its Applications.Springer, New York (2002). MR1876169. doi:10.1007/978-1-4757-4015-8Ledoux, M., Talagrand, M.: Probability in Banach Spaces. Springer, Berlin (1991).Lifshits, M.: Lectures on Gaussian Processes. Springer, New York (2012). MR3024389. doi:10.1007/978-3-642-24939-6McKean, H.P.: Geometry of differential space. Ann. Probab.1 (2), 197–206 (1973). MR0353471Mishura, Y.: Stochastic Calculus for Fractional Brownian Motion and Related Processes. Springer, Berlin (2008). MR2378138. doi:10.1007/978-3-540-75873-0Naor, A., Romik, D.: Projecting the surface measure of the sphere of ℓpn.Ann. Inst. Henri Poincaré Probab. Stat.39 (2), 241–261 (2003). MR1962135. doi:10.1016/S0246-0203(02)00008-0Parczewski, P.: A fractional donsker theorem. Stoch. Anal. Appl.32 (2), 328–347 (2014). MR3177074. doi:10.1080/07362994.2013.866521Rachev, S.T., Rüschendorf, L.: Approximate independence of distributions on spheres and their stability properties. Ann. Probab.19 (3), 1311–1337 (1991). MR1112418Račkauskas, A., Suquet, C.: Convergence of self-normalized partial sums processes in c[0,1] and d[0,1]. Publications IRMA de Lille 53 - VI (2000).Schechtman, G., Zinn, J.: On the volume of the intersection of two lpn balls.Proc. Am. Math. Soc.110 (1), 217–224 (1990). MR1015684. doi:10.2307/2048262Taqqu, M.S.: Weak convergence to fractional brownian motion and to the rosenblatt process.Z. Wahrscheinlichkeitstheor. Verw. Geb.31, 287–302 (1975). MR0400329. doi:10.1007/BF00532868Vershik, A.M.: Does a lebesgue measure in an infinite-dimensional space exist?Tr. Mat. Inst. Steklova259, 256–281 (2007). MR2433687. doi:10.1134/S0081543807040153Wiener, N.: Differential space. J. Math. Phys.2, 132–174 (1923).