Interpretable statistics for complex modelling: quantile and topological learning

As the complexity of our data increased exponentially in the last decades, so has our need for interpretable features. This thesis revolves around two paradigms to approach this quest for insights. In the first part we focus on parametric models, where the problem of interpretability can be seen as a “parametrization selection”. We introduce a quantile-centric parametrization and we show the advantages of our proposal in the context of regression, where it allows to bridge the gap between classical generalized linear (mixed) models and increasingly popular quantile methods. The second part of the thesis, concerned with topological learning, tackles the problem from a non-parametric perspective. As topology can be thought of as a way of characterizing data in terms of their connectivity structure, it allows to represent complex and possibly high dimensional through few features, such as the number of connected components, loops and voids. We illustrate how the emerging branch of statistics devoted to recovering topological structures in the data, Topological Data Analysis, can be exploited both for exploratory and inferential purposes with a special emphasis on kernels that preserve the topological information in the data. Finally, we show with an application how these two approaches can borrow strength from one another in the identification and description of brain activity through fMRI data from the ABIDE project

Polar mesoscale cyclones in the northeast Atlantic: Comparing climatologies from ERA-40 and satellite imagery

Polar mesoscale cyclones over the subarctic are thought to be an important component of the coupled atmosphere–ocean climate system. However, the relatively small scale of these features presents some concern as to their representation in the meteorological reanalysis datasets that are commonly used to drive ocean models. Here polar mesocyclones are detected in the 40-Year European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis dataset (ERA-40) in mean sea level pressure and 500-hPa geopotential height, using an automated cyclone detection algorithm. The results are compared to polar mesocyclones detected in satellite imagery over the northeast Atlantic, for the period October 1993–September 1995. Similar trends in monthly cyclone numbers and a similar spatial distribution are found. However, there is a bias in the size of cyclones detected in the reanalysis. Up to 80% of cyclones larger than 500 km are detected in MSL pressure, but this hit rate decreases, approximately linearly, to ∼40% for 250-km-scale cyclones and to ∼20% for 100-km-scale cyclones. Consequently a substantial component of the associated air–sea fluxes may be missing from the reanalysis, presenting a serious shortcoming when using such reanalysis data for ocean modeling simulations. Eight maxima in cyclone density are apparent in the mean sea level pressure, clustered around synoptic observing stations in the northeast Atlantic. They are likely spurious, and a result of unidentified shortcomings in the ERA-40 data assimilation procedure

Multi-Material Mesh Representation of Anatomical Structures for Deep Brain Stimulation Planning

The Dual Contouring algorithm (DC) is a grid-based process used to generate surface meshes from volumetric data. However, DC is unable to guarantee 2-manifold and watertight meshes due to the fact that it produces only one vertex for each grid cube. We present a modified Dual Contouring algorithm that is capable of overcoming this limitation. The proposed method decomposes an ambiguous grid cube into a set of tetrahedral cells and uses novel polygon generation rules that produce 2-manifold and watertight surface meshes with good-quality triangles. These meshes, being watertight and 2-manifold, are geometrically correct, and therefore can be used to initialize tetrahedral meshes. The 2-manifold DC method has been extended into the multi-material domain. Due to its multi-material nature, multi-material surface meshes will contain non-manifold elements along material interfaces or shared boundaries. The proposed multi-material DC algorithm can (1) generate multi-material surface meshes where each material sub-mesh is a 2-manifold and watertight mesh, (2) preserve the non-manifold elements along the material interfaces, and (3) ensure that the material interface or shared boundary between materials is consistent. The proposed method is used to generate multi-material surface meshes of deep brain anatomical structures from a digital atlas of the basal ganglia and thalamus. Although deep brain anatomical structures can be labeled as functionally separate, they are in fact continuous tracts of soft tissue in close proximity to each other. The multi-material meshes generated by the proposed DC algorithm can accurately represent the closely-packed deep brain structures as a single mesh consisting of multiple material sub-meshes. Each sub-mesh represents a distinct functional structure of the brain. Printed and/or digital atlases are important tools for medical research and surgical intervention. While these atlases can provide guidance in identifying anatomical structures, they do not take into account the wide variations in the shape and size of anatomical structures that occur from patient to patient. Accurate, patient-specific representations are especially important for surgical interventions like deep brain stimulation, where even small inaccuracies can result in dangerous complications. The last part of this research effort extends the discrete deformable 2-simplex mesh into the multi-material domain where geometry-based internal forces and image-based external forces are used in the deformation process. This multi-material deformable framework is used to segment anatomical structures of the deep brain region from Magnetic Resonance (MR) data

Constructing Desirable Scalar Fields for Morse Analysis on Meshes

Morse theory is a powerful mathematical tool that uses the local differential properties of a manifold to make conclusions about global topological aspects of the manifold. Morse theory has been proven to be a very useful tool in computer graphics, geometric data processing and understanding. This work is divided into two parts. The first part is concerned with constructing geometry and symmetry aware scalar functions on a triangulated $2$-manifold. To effectively apply Morse theory to discrete manifolds, one needs to design scalar functions on them with certain properties such as respecting the symmetry and the geometry of the surface and having the critical points of the scalar function coincide with feature or symmetry points on the surface. In this work, we study multiple methods that were suggested in the literature to construct such functions such as isometry invariant scalar functions, Poisson fields and discrete conformal factors. We suggest multiple refinements to each family of these functions and we propose new methods to construct geometry and symmetry-aware scalar functions using mainly the theory of the Laplace-Beltrami operator. Our proposed methods are general and can be applied in many areas such as parametrization, shape analysis, symmetry detection and segmentation. In the second part of the thesis we utilize Morse theory to give topologically consistent segmentation algorithms