Search papers, labs, and topics across Lattice.
This paper computes the maximum likelihood threshold (MLT) for colored Gaussian graphical models (CGGMs), which are designed to handle high-dimensional data with small sample sizes by leveraging symmetry constraints. By developing a geometric framework to determine the minimum rank of sample covariance matrices necessary for the existence of the MLT, the authors extend existing results from uncolored to colored models. Their findings highlight the effectiveness of topological data analysis (TDA) in addressing computational challenges associated with traditional methods, revealing new insights into the likelihood geometry of CGGMs.
TDA can dramatically streamline the computation of maximum likelihood thresholds in colored Gaussian graphical models, revealing insights previously obscured by traditional methods.
Gaussian graphical models (GGMs) are essential tools for interpretable structure learning. However, in high-dimensional, small-sample regimes, the available data is often insufficient for the maximum likelihood estimator to exist. Colored Gaussian graphical models (CGGMs) mitigate this limitation by imposing symmetry constraints through graph coloring, which reduces the required sample size. This minimal number of observations needed to guarantee that the estimator exists almost surely is defined as the maximum likelihood threshold (MLT). Here, we address the computation of the MLT for CGGMs by focusing on its geometric formulation: finding the minimum rank of a sample covariance matrix such that its projection lies almost surely within the interior of the cone of sufficient statistics. We establish a unified theoretical framework, extending results from uncolored to colored models and introducing new symbolic algorithms. Furthermore, we present a computational study integrating sampling with topological data analysis (TDA) to investigate the local geometry of the cone of sufficient statistics. Our results demonstrate the potential of TDA to overcome the computational bottlenecks of traditional symbolic algebraic methods, particularly Groebner basis computations, in analyzing the likelihood geometry of CGGMs.