arxiv: v1 [math.pr] 29 Dec 2015

Transcription

1 Royen s proof of the Gaussian correlation inequality arxiv:5.8776v [math.pr] 9 Dec 5 Rafa l Lata la and Dariusz Matlak Abstract We present in detail Thomas Royen s proof of the Gaussian correlation inequality which states that µk L µkµl for any centered Gaussian measure µ on R d and symmetric convex sets K,L in R d. Introduction Theaim of this note is to present in aself contained way thebeautiful proofof the Gaussian correlation inequality, due to Thomas Royen [7]. Although the method is rather simple and elementary, we found the original paper not too easy to follow. One of the reasons behind it is that in [7] the correlation inequality was established for more general class of probability measures. Moreover, the author assumed that the reader is familiar with properties of certain distributions and may justify some calculations by herself/himself. We decided to reorganize a bit Royen s proof, restrict it only to the Gaussian case and add some missing details. We hope that this way a wider readership may appreciate the remarkable result of Royen. The statement of the Gaussian correlation inequality is as follows. Theorem. For any closed symmetric sets K,L in R d and any centered Gaussian measure µ on R d we have µk L µkµl. For d = the result was proved by Pitt [5]. In the case when one of the sets K,L is a symmetric strip which corresponds to min{n,n } = in Theorem below inequality was established independently by Khatri [3] and Šidák [9]. Hargé [] generalized the Khatri-Šidak result to the case when one of the sets is a symmetric ellipsoid. Some other partial results may be found in papers of Borell [] and Schechtman, Schlumprecht and Zinn [8]. Up to our best knowledge Thomas Royen was the first to present a complete proof of the Gaussian correlation inequality. Some other recent attempts may be found in [4] and [6], however both papers are very long and difficult to check. The first version of [4], placed on the arxiv before Royen s paper, contained a fundamental mistake Lemma 6.3 there was wrong.

2 Since any symmetric closed set is a countable intersection of symmetric strips, it is enough to show in the case when K = {x R d : i n x,v i t i } and L = {x R d : n + i n +n x,v i t i }, where v i are vectors in R d and t i nonnegative numbers. If we set n = n +n, X i := v i,g, where G is the Gaussian random vector distributed according to µ, we obtain the following equivalent form of Theorem. Theorem. Let n = n +n and X be an n-dimensional centered Gaussian vector. Then for any t,...,t n >, P X t,..., X n t n P X t,..., X n t n P X n + t n +,..., X n t n. Remark 3. i The standard approximation argument shows that the Gaussian correlation inequality holds for centered Gaussian measures on separable Banach spaces. ii Thomas Royen established Theorem for more general class of random vectors X such that X = X,...,X n has an n-variate gamma distribution see [7] for details. Notation. By N, C we denote the centered Gaussian measure with the covariance matrix C. We write M n m for a set of n m matrices and A for the determinant of a square matrix A. For a matrix A = a ij i,j n and J [n];= {,...,n} by A J we denote the square matrix a ij i,j J and by J the cardinality of J. Proof of Theorem Without loss of generality we may and will assume that the covariance matrix C of X is nondegenerate i.e. strictly positively defined. We may write C as C C C =, C C where C ij is the n i n j matrix. Let C τc Cτ := τc C, τ. Set Z i τ := X iτ, i n, where Xτ N,Cτ. We may restate the assertion as PZ s,...,z n s n PZ s,...,z n s n, where s = t i. Therefore it is enough to show that the function τ PZ τ s,...,z n τ s n is nondecreasing on [,].

3 Let fx,τ denote the density of the random vector Zτ and K = [,s ] [,s n ]. We have τ PZ τ s,...,z n τ s n = fx,τdx = τ K K τ fx,τdx, where the last equation follows by Lemma 6 applied to λ =... = λ n =. Therefore it is enough to show that K τfx,τ. To this end we will compute the Laplace transform of τfx,τ. By Lemma 6, applied to K = [, n, we have for any λ...,λ n, τ fx,τdx = fx,τdx. τ [, n e However by Lemma 4 we have fx,τdx = Eexp [, n e where Λ = diagλ,...,λ n. Formula below yields I +ΛCτ = + [, n e λ k Xk τ = I +ΛCτ /, k= ΛCτ J = + Cτ J j Jλ j. Fix J [n]. Then J = J J, where J := [n ] J, J := J \ [n ] and Cτ J = CJ τc J J. If J τc J J C = or J = then Cτ J = C J, otherwise by 3 we get J Cτ J = C J C J I J τ C / C J J CJ C J J C / J J = C J C J τ µ J,J i, where µ J,J i, i J denote the eigenvalues of C / J C J J C J C J J C / J by 4 they belong to [,]. Thus for any J [n] and τ [,] we have Therefore a J τ := τ Cτ J. τ I +ΛCτ / = I +ΛCτ 3/ = I +ΛCτ 3/ a J τ j J J τ Cτ J Λ J λ j. 3

4 We have thus shown that τ fx,τdx = [, n e a Jτ I +ΛCτ 3/ λ j. j J Let h τ := h 3,Cτ be the density function on, n defined by 5. By Lemmas 8 and 7 iii we know that This shows that thus I +ΛCτ 3/ j = j Jλ τ fx,τ =, n e J x J h τ. a Jτ J x J h τ x. Finally recall that a J τ and observe that by Lemma 7 ii, K I lim h τ x for i / I [n], x i + x I J h τ xdx = x h τ t J,x J cdx J c, J j J c[,t j] where J c = [n]\j and y = t J,x J c if y i = t i for i J and y i = x i for i J c. 3 Auxiliary Lemmas Lemma 4. Let X be an n dimensional centered Gaussian vector with the covariance matrix C. Then for any λ,...,λ n we have Eexp λ i Xi = I n +ΛC /, where Λ := diagλ,...,λ n. Proof. Let A be a symmetric positively defined matrix. Then A = UDU T for some U On and D = diagd,d,...,d n. Hence exp Ax,x dx = exp Dx,x dx = R n R n n k= π d k = π n/ D / = π n/ A /. 4

5 Therefore for a canonical Gaussian vector Y N,I n and a symmetric matrix B such that B < I n we have Eexp BY,Y = π R n/ exp n I n B x,x dx = n/ I / n B = I n B /. We may represent X N,C as X AY with Y N,I n and C = AA T. Thus Eexp λ i Xi = Eexp ΛX,X = Eexp ΛAY,AY = Eexp A T ΛAY,Y = I n +A T ΛA / = I n +ΛC /, where to get the last equality we used the fact that I n +A A = I n +A A for A,A n M n n. Lemma 5. i For any matrix A M n n, I n +A = + A J. A A ii Suppose that n = n +n and A M n n has the block representation A =, A A where A ij M ni n j and A, A are invertible. Then A = A A I n A / A A. 3 Moreover, if A is symmetric and positively defined then A A / A / A A A A / I n. 4 Proof. i This formula may be verified in several ways e.g. by induction on n or by using the Leibniz formula for the determinant. ii We have A A A A = A / A / I n A / A A / A / A A / I n A / A / and I n A / A A / A / A A / I n = = I n A / A A A A / A / A A / I n A / A A I n A A /. 5

6 To show the last part of the statement notice that A / A A A A / = B T B, where B := A / A A /. If A is positively defined then for any t R, x R n and y R n we have t A x,x + t A x,y + A y,y. This implies A x,y A x,x A y,y. Replacing x by A / x and y by A / y we get Bx,y x y. Choosing y = Bx we get B T Bx,x x, i.e. B T B I n. Lemma 6. Let fx,τ be the density of the random vector Zτ defined above. Then for any Borel set K in [, n and any λ,...,λ n, e τ fx,τdx = e fx,τdx. τ K Proof. The matrix C is nondegenerate, therefore matrices C and C are nondegerate and Cτ is nondegenerate for any τ [,]. Random vector Xτ N,Cτ has the density Cτ / π n/ exp Cτ x,x. Standard calculation shows that Zτ has the density fx,τ = Cτ / 4π n/ x x n K ε {,} n e Cτ where for ε {,} n and x, n we set x ε := ε i xi i. The function τ Cτ / is smooth on [,], in particular sup Cτ / + sup τ [,] τ [,] τ Cτ / =: M <. Since Cτ = τc+ τc we have τcτ = C C and x ε, x ε, nx, e Cτ xε, xε = Cτ C CCτ x τ ε, x ε e Cτ x ε, x ε. The continuity of the function τ Cτ gives and Cτ x ε, x ε a x ε, x ε = a x i Cτ C CCτ x ε, x ε b x ε, x ε = b for some a >, b <. Hence for x, n sup τ fx,τ gx := Mπ n/ +b x x n τ [,] x i x i e a xi. Since gx L, n and e the statement easily follows by the Lebesgue dominated convergence theorem. 6

7 Let for α >, g α x,y := e x y k= x k+α y k Γk+α k! x >,y. For µ,α,...,α n > and a random vector Y = Y,...,Y n such that PY i = we set [ n ] h α,...,α n,µ,yx,...,x n := E µ g xi α i µ,y i, x,...,x n >. Lemma 7. Let µ > and Y be a random n-dimensional vector with nonnegative coordinates. For α = α,...,α n, n set h α := f α,...,α n,µ,y. i For any α, n, h α and, n h α xdx =. ii If α, n and α i > then lim xi +h α x =, x i h α x exists and iii If α, n then for any J [n], Moreover for λ,...,λ n, e J, n x i h α x = h α ei h α. x J h α xdx = i J J x J h α x exists and belongs to L, n. λ i, n e Proof. i Obviously h α [, ]. We have for any y and α >, x µ g α µ,y dx = g α x,ydx =. Hence by the Fubini theorem,, n h α xdx = E k j= h α xdx. µ g xi α i µ,y i dx i =. ii It is well known that Γx is decreasing on,x ] and increasing on [x,, where < x < and Γx > /. Therefore for k =,... and α >, Γk + α Γk = k! and g α x,y e x k= x k+α Γk+α x α e x +x α x k k! e x = x α e x +x. k= 7

8 This implies that for α > and < a < b <, g α x,y Cα,a,b < for x a,b and y. Moreover, h α x µ n n αi xi + x i. µ µ In particular lim xi +h α x = if α i. Observe that for α >, x g α = g α g α. Standard application of the Lebegue dominated convergence theorem concludes the proof of part ii. iii By ii we get J h α = J i J δ i f x α i J δ ie i L, n. J δ {,} J Moreover lim xj + J x J h α x = for j / J. We finish the proof by induction on J using integration by parts. Let C be a strictly positively defined symmetric n n matrix. Then there exists µ > such that C µi n is positively defined, so C = µi n + AA T for some A M n n. Let g l j j n,l k be i.i.d. N, r.v s. Set and Y i = µ k l= j,j n g l j gl j a i,j a i,j = k g l µ l= j= j a i,j, i n h k,c := hk,...,k,µ,y. 5 Lemma 8. For any λ,...,λ n we have h k,c x = I n +ΛC k, where Λ = diagλ,...,λ n., n e Proof. We have for any α,µ > and λ,y x µ e λx g α µ,y dx = e y y k k!γk +α k= = e y y k k= e λ+ µ x x k+α dx µ k+α k!+µλ k+α = +µλ α e µλ +µλ y. 8

9 By the Fubini theorem we have, n e h k,c xdx = E n = I n +µλ k Ee e λ ix i µ g xi k/ µ,y i dx i µλ i +µλ i Y i. Observe that Y i = k l= Xl i, where X l := X l i i n are independent N, µ AAT. Therefore by Lemma 4 we have h k,c xdx = I n +µλ k I +µλi +µλ µ AAT k = In +C k., n e References [] C. Borell, A Gaussian correlation inequality for certain bodies in R n, Math. Ann , [] G. Hargé, A particular case of correlation inequality for the Gaussian measure, Ann. Probab , [3] C. G. Khatri, On certain inequalities for normal distributions and their applications to simultaneous confidence bounds, Ann. Math. Statist , [4] Y. Memarian, The Gaussian Correlation Conjecture Proof, arxiv: [5] L. D. Pitt, A Gaussian correlation inequality for symmetric convex sets, Ann. Probab , [6] G. Qingyang, The Gaussian Correlation Inequality for Symmetric Convex Sets, arxiv:.676. [7] T. Royen, A simple proof of the Gaussian correlation conjecture extended to multivariate gamma distributions, Far East J. Theor. Stat. 48 4, [8] G. Schechtman, T. Schlumprecht, J. Zinn, On the Gaussian measure of the intersection, Ann. Probab , [9] Z. Šidák, Rectangular confidence regions for the means of multivariate normal distributions, J. Amer. Statist. Assoc , Institute of Mathematics University of Warsaw Banacha -97 Warszawa, Poland rlatala@mimuw.edu.pl, ddmatlak@gmail.com 9