M-estimation for location and regression parameters in group models: A case study using Stiefel manifolds

摘要

We discuss here a general approach to the calculation of the asymptotic covariance of M-estimates for location parameters in statistical group models when an invariant objective function is used. The calculation reduces to standard tools in group representation theory and the calculation of some constants. Only the constants depend upon the precise forms of the density or of the objective function. If the group is sufficiently large this represents a major simplification in the computation of the asymptotic covariance. Following the approach of Chang and Tsai we define a regression model for group models and derive the asymptotic distribution of estimates in the regression model from the corresponding distribution theory for the location model. The location model is not, in general, a subcase of the regression model. We illustrate these techniques using Stiefel manifolds. The Stiefel manifold V-p,V-m is the collection of p x m matrices X which satisfy the condition (XX)-X-T = I-m where m less than or equal to p. Under the assumption that X has a distribution proportional to exp(Tr((FX)-X-T)), for some p x m matrix F, Downs (1972) gives approximations to maximum likelihood estimation of F. In this paper, we consider a somewhat different location problem: under the assumption that X has a distribution of the form f(Tr(theta X-T(0))) for some 0(0) is an element of V-p,V-m, we calculate the asymptotic distribution of M-estimates which minimize an objective function of the form Sigma (i) rho (Tr(0(T)X(i))). The assumptions on the form of the density and the objective function can be relaxed to a more general invariant form. In this case, the calculation of the asymptotic distribution of <()over cap> reduces to the calculation of four constants and we present consistent estimators for these constants. Prentice (1989) introduced a regression model for Stiefel manifolds. In the Prentice model, u(1), u(2), . . ., u(n) is an element of V-p,V-m are fixed, V-1, V-2,..., V-n is an element of V-p,V-m are independent random so that the distribution of Vi depends only upon Tr(V(i)(T)A(2)u(i)A(1)(T)) for unknown (A(1), A(2)) is an element of SO(m) x SO(p). We discuss here M-estimation of A(1) and A(2) under general invariance conditions for both the density and the objective function. Using a well-studied example on vector cardiograms we discuss the physical interpretation of the invariance assumption as well as of the parameters (A(1),A(2)) in the Prentice regression model. In particular, A(1) represents a rotation of the u's to the V's in a coordinate system relative to the u's and A(2) represents a rotation of the u's to the V's in a coordinate system fixed to the external world. [References: 21]