Sunday, 28 November 2010

pr.probability - Limit law for the number of local maxima in a square lattice of IID random variables

For i,jin1,ldots,n, let Xi,j be a real-valued random variable uniformly distributed on the interval [0,1]. The Xi,j are independent.



Let Ai,j be the indicator random variable of the event that Xi,j is a local maximum, i. e. it is the largest of the five random variables Xi,j,Xi,j+1,Xi,j1,Xi1,j,Xi+1,j. For the sake of not having to think about boundary conditions, interpret all coordinates modulo n.



Then it's clear that EAi,j=1/5, by symmetry. It's also not hard to see that:



  • E(Ai,jAi,j+1)=0, since we can't have local maxima both at (i,j) and at (i,j+1). (The case where Xi,j=Xi,j+1 can be ignored since it occurs with probability zero.


  • I believe E(Ai,jAi,j+2)=2/45 and E(Ai,jAi+1,j+1)=1/20. (In any case, these are clearly constants, and their exact values don't matter.)


  • E(Ai,jAk,l)=1/25 for all choices of i,j,k,l other than the ones I already listed and those that clearly are related to them by symmetry. That is, Ai,j and Ak,l are independent unless the rook-neighborhoods of (i,j) and (k,l) overlap.


From this we can compute the mean and variance of Mn=sumni=1sumnj=1Ai,j, the total number of maxima. Of course EM=n2/5. The variance is a bit harder to compute and I haven't actually written out the computation, but it ought to be asymptotic to c2n2 for some positive c. (It is the sum of Theta(n2) covariances each of order 1.)



Let tildeMn=(Mnn2/5)/(cn). Then tildeMn has mean 0 and variance 1. Is it true that the sequence tildeMnin=1nfty converges in distribution to the standard normal? Intuitively this seems like it should be true -- we're adding up a bunch of small, almost-independent contributions. If it's true, how can this be proven?



This problem came from last year's qualifying exam in probability at Penn; it's been making the rounds around here over the past few days but nobody seems to remember how to do it.

No comments:

Post a Comment