We study properties of Fisher distribution (von Mises-Fisher distribution,
matrix Langevin distribution) on the rotation group SO(3). In particular we
apply the holonomic gradient descent, introduced by Nakayama et al. (2011), and
a method of series expansion for evaluating the normalizing constant of the
distribution and for computing the maximum likelihood estimate. The rotation
group can be identified with the Stiefel manifold of two orthonormal vectors.
Therefore from the viewpoint of statistical modeling, it is of interest to
compare Fisher distributions on these manifolds.
We give a new algorithm to find local maximum and minimum of a holonomic
function and apply it for the Fisher-Bingham integral on the sphere $S^n$,
which is used in the directional statistics. The method utilizes the theory and
algorithms of holonomic systems.