The uniform distribution on a Stiefel manifold is a matrix-variate distribution that plays an important role in multivariate statistics. There one often encounters integrals over the orthogonal group or over the Stiefel manifold with respect to an invariant measure. For example, this distribution arises in the study of the functional determinant under transformations involving orthogonal or semi-orthogonal matrices. The uniform distribution on the Stiefel manifold corresponds to the normalized Haar measure on the Stiefel manifold. A random matrix uniformly distributed on the Stiefel manifold is invariant under the two-sided group action of the product O ( p ) × O ( n ) {\displaystyle O(p)\times O(n)} of orthogonal groups, i.e. X ∼ V 1 X V 2 {\displaystyle X\sim V_{1}XV_{2}} for all V 1 ∈ O ( p ) {\displaystyle V_{1}\in O(p)} and V 2 ∈ O ( n ) {\displaystyle V_{2}\in O(n)} .
Uniform Distribution on a Stiefel Manifold
Introduction Let V p , n := V n ( R p ) {\displaystyle V_{p,n}:=V_{n}(\mathbb {R} ^{p})} be the Stiefel manifold, i.e., the set of all orthonormal n {\displaystyle n} -frames in R p {\displaystyle \mathbb {R} ^{p}} for n ≤ p {\displaystyle n\leq p} . This manifold can also be represented as the matrix set
V p , n = { X ∈ R p × n : X ′ X = I n } {\displaystyle V_{p,n}=\{X\in \mathbb {R} ^{p\times n}\colon X'X=I_{n}\}} . The Stiefel manifold is homeomorphic to the quotient space of the orthogonal groups
V p , n ≅ O ( p ) / O ( p − n ) . {\displaystyle V_{p,n}\cong O(p)/O(p-n).} These two can be identified, and in the case p = n {\displaystyle p=n} we obtain the full orthogonal group. The Stiefel manifold inherits the left group action
X ↦ V X , V ∈ O ( p ) . {\displaystyle X\mapsto VX,\quad V\in O(p).}
Here, O ( p − n ) {\displaystyle O(p-n)} is a compact, closed Lie subgroup of O ( p ) {\displaystyle O(p)} . By Haar's theorem there exists a Haar measure on O ( p ) {\displaystyle O(p)} which induces an invariant measure on the quotient space O ( p ) / O ( p − n ) {\displaystyle O(p)/O(p-n)} .
Derivation of the Haar Measure on the Stiefel Manifold Let X ∈ O ( p ) {\displaystyle X\in O(p)} . Differentiating X ′ X = I p {\displaystyle X'X=I_{p}} yields: d X ′ X + X ′ d X = 0. {\displaystyle dX'X+X'dX=0.} Let x 1 , … , x p {\displaystyle x_{1},\dots ,x_{p}} be the columns of X = ( x 1 , … , x p ) {\displaystyle X=(x_{1},\dots ,x_{p})} . The exterior product of the superdiagonal elements defines a differential form
… excerpt ends here. Continue reading the full article.
