Robustness of Lognormal Confidence Regions for Means of Symmetric Positive Matrices

• Activity Number: 339 - Model-Fitting, Likelihood-Based Inference, and Their Applications
• Type: Contributed
• Date/Time: Tuesday, August 1, 2017 : 10:30 AM to 12:20 PM
• Sponsor: IMS
• Author(s): Benoit Ahanda* and Daniel E Osborne and Leif Ellingson
• Companies: Texas Tech University and Florida A&M University and Texas Tech University

Abstract

Symmetric positive definite (SPD) matrices arise in a wide range of applications including diffusion tensor imaging (DTI), cosmic background radiation, and as covariance matrices. A complication when working with such data is that the space of SPD matrices is a manifold, so traditional statistical methods may not be directly applied. However, there are nonparametric procedures based on resampling for statistical inference for such data, but these can be slow and computationally tedious. Schwartzman (2015) introduced a lognormal distribution on the space of SPD matrices, providing a convenient framework for parametric inference on this space. Our goal is to check how robust confidence regions based on this distributional assumption are to a lack of lognormality. The methods are illustrated in a simulation study by examining the coverage probability of various mixtures of distributions.