movMF: An R Package for Fitting Mixtures of von Mises-Fisher Distributions

Finite mixtures of von Mises-Fisher distributions allow to apply model-based clustering methods to data which is of standardized length, i.e., all data points lie on the unit sphere. The R package movMF contains functionality to draw samples from finite mixtures of von Mises-Fisher distributions and to fit these models using the expectation-maximization algorithm for maximum likelihood estimation. Special features are the possibility to use sparse matrix representations for the input data, different variants of the expectation-maximization algorithm, different methods for determining the concentration parameters in the M-step and to impose constraints on the concentration parameters over the components. In this paper we describe the main fitting function of the package and illustrate its application. In addition we compare the clustering performance of finite mixtures of von Mises-Fisher distributions to spherical k-means. We also discuss the resolution of several numerical issues which occur for estimating the concentration parameters and for determining the normalizing constant of the von Mises-Fisher distribution

Computational Representation of White Matter Fiber Orientations

We present a new methodology based on directional data clustering to represent white matter fiber orientations in magnetic resonance analyses for high angular resolution diffusion imaging. A probabilistic methodology is proposed for estimating intravoxel principal fiber directions, based on clustering directional data arising from orientation distribution function (ODF) profiles. ODF reconstructions are used to estimate intravoxel fiber directions using mixtures of von Mises-Fisher distributions. The method focuses on clustering data on the unit sphere, where complexity arises from representing ODF profiles as directional data. The proposed method is validated on synthetic simulations, as well as on a real data experiment. Based on experiments, we show that by clustering profile data using mixtures of von Mises-Fisher distributions it is possible to estimate multiple fiber configurations in a more robust manner than currently used approaches, without recourse to regularization or sharpening procedures. The method holds promise to support robust tractographic methodologies and to build realistic models of white matter tracts in the human brain

Recent advances in directional statistics

Mainstream statistical methodology is generally applicable to data observed in Euclidean space. There are, however, numerous contexts of considerable scientific interest in which the natural supports for the data under consideration are Riemannian manifolds like the unit circle, torus, sphere and their extensions. Typically, such data can be represented using one or more directions, and directional statistics is the branch of statistics that deals with their analysis. In this paper we provide a review of the many recent developments in the field since the publication of Mardia and Jupp (1999), still the most comprehensive text on directional statistics. Many of those developments have been stimulated by interesting applications in fields as diverse as astronomy, medicine, genetics, neurology, aeronautics, acoustics, image analysis, text mining, environmetrics, and machine learning. We begin by considering developments for the exploratory analysis of directional data before progressing to distributional models, general approaches to inference, hypothesis testing, regression, nonparametric curve estimation, methods for dimension reduction, classification and clustering, and the modelling of time series, spatial and spatio-temporal data. An overview of currently available software for analysing directional data is also provided, and potential future developments discussed.Comment: 61 page