97
2

The "north pole problem" and random orthogonal matrices

Abstract

This paper is motivated by the following observation. Take a 3 x 3 random (Haar distributed) orthogonal matrix Γ\Gamma, and use it to "rotate" the north pole, x0x_0 say, on the unit sphere in R3R^3. This then gives a point u=Γx0u=\Gamma x_0 that is uniformly distributed on the unit sphere. Now use the same orthogonal matrix to transform u, giving v=Γu=Γ2x0v=\Gamma u=\Gamma^2 x_0. Simulations reported in Marzetta et al (2002) suggest that v is more likely to be in the northern hemisphere than in the southern hemisphere, and, morever, that w=Γ3x0w=\Gamma^3 x_0 has higher probability of being closer to the poles ±x0\pm x_0 than the uniformly distributed point u. In this paper we prove these results, in the general setting of dimension p3p\ge 3, by deriving the exact distributions of the relevant components of u and v. The essential questions answered are the following. Let x be any fixed point on the unit sphere in RpR^p, where p3p\ge 3. What are the distributions of U2=xΓ2xU_2=x'\Gamma^2 x and U3=xΓ3xU_3=x'\Gamma^3 x? It is clear by orthogonal invariance that these distribution do not depend on x, so that we can, without loss of generality, take x to be x0=(1,0,...,0)Rpx_0=(1,0,...,0)'\in R^p. Call this the "north pole". Then x0Γkx0x_0'\Gamma^ k x_0 is the first component of the vector Γkx0\Gamma^k x_0. We derive stochastic representations for the exact distributions of U2U_2 and U3U_3 in terms of random variables with known distributions.

View on arXiv
Comments on this paper