Saddlepoint approximations for the normalizing constant of Fisher--Bingham distributions on products of spheres and Stiefel manifolds

In an earlier paper Kume & Wood (2005) showed how the normalizing constant of the Fisher--Bingham distribution on a sphere can be approximated with high accuracy using a univariate saddlepoint density approximation. In this sequel, we extend the approach to a more general setting and derive saddlepoint approximations for the normalizing constants of multicomponent Fisher--Bingham distributions on Cartesian products of spheres, and Fisher--Bingham distributions on Stiefel manifolds. In each case, the approximation for the normalizing constant is essentially a multivariate saddlepoint density approximation for the joint distribution of a set of quadratic forms in normal variables. Both first-order and second-order saddlepoint approximations are considered. Computational algorithms, numerical results and theoretical properties of the approximations are presented. In the challenging high-dimensional settings considered in this paper the saddlepoint approximations perform very well in all examples considered. Copyright 2013, Oxford University Press.

[1]  J. Kent.  The Complex Bingham Distribution and Shape Analysis , 1994 .

[2]  Ronald W. Butler,et al.  Laplace approximation for Bessel functions of matrix argument , 2003 .

[3]  Peter D. Hoff,et al.  Simulation of the Matrix Bingham–von Mises–Fisher Distribution, With Applications to Multivariate and Relational Data , 2007, 0712.4166.

[4]  R. Muirhead.  Aspects of Multivariate Statistical Theory , 1982, Wiley Series in Probability and Statistics.

[5]  Harshinder Singh,et al.  Probabilistic model for two dependent circular variables , 2002 .

[6]  M. J. Prentice.  Orientation Statistics Without Parametric Assumptions , 1986 .

[7]  Kanti V. Mardia,et al.  Statistics of Directional Data , 1972 .

[8]  A. Wood,et al.  Saddlepoint approximations for the Bingham and Fisher–Bingham normalising constants , 2005 .

[9]  Michael R. Pierrynowski,et al.  A Directional Model for the Estimation of the Rotation Axes of the Ankle Joint , 2008 .

[10]  K. Mardia,et al.  A general correlation coefficient for directional data and related regression problems , 1980 .

[11]  K. Mardia,et al.  Maximum Likelihood Estimators for the Matrix Von Mises-Fisher and Bingham Distributions , 1979 .

[12]  Kanti V. Mardia,et al.  A multivariate von mises distribution with applications to bioinformatics , 2008 .

[13]  Karim Oualkacha,et al.  On the estimation of an average rigid body motion , 2012 .

[14]  Andrew T. A. Wood,et al.  ESTIMATION OF THE CONCENTRATION PARAMETERS OF THE FISHER MATRIX DISTRIBUTION ON 50(3) AND THE BINGHAM DISTRIBUTION ON Sq, q≥ 2 , 1993 .

[15]  D. Flinn.  Orientation Statistics , 1967, Nature.

[16]  J. Kent.  The Fisher‐Bingham Distribution on the Sphere , 1982 .

[17]  K. Mardia,et al.  The von Mises–Fisher Matrix Distribution in Orientation Statistics , 1977 .

[18]  Christopher Bingham.  An Antipodally Symmetric Distribution on the Sphere , 1974 .

[19]  K. Mardia.  Statistics of Directional Data , 1972 .

[20]  P. Jupp,et al.  Statistics of orthogonal axial frames , 2013 .