Maximum likelihood thresholds via graph rigidity

Abstract

DIB was partially supported by a Mathematical Sciences Postdoctoral Research Fellowship from the US NSF Grant DMS-1802902. SD was partially supported by the Austrian Science Fund (FWF): P31888. AN was partially supported by the Heilbronn Institute for Mathematical Research. SJG was partially supported by US NSF Grant DMS-1564473. MS was partially supported by US NSF Grant DMS-1564480 and US NSF Grant DMS-1563234.The maximum likelihood threshold (MLT) of a graph G is the minimum number of samples to almost surely guarantee existence of the maximum likelihood estimate in the corresponding Gaussian graphical model. We give a new characterization of the MLT in terms of rigidity-theoretic properties of G and use this characterization to give new combinatorial lower bounds on the MLT of any graph. We use the new lower bounds to give high-probability guarantees on the maximum likelihood thresholds of sparse Erd{ö}s-Rényi random graphs in terms of their average density. These examples show that the new lower bounds are within a polylog factor of tight, where, on the same graph families, all known lower bounds are trivial. Based on computational experiments made possible by our methods, we conjecture that the MLT of an Erd{ö}s-Rényi random graph is equal to its generic completion rank with high probability. Using structural results on rigid graphs in low dimension, we can prove the conjecture for graphs with MLT at most 4 and describe the threshold probability for the MLT to switch from 3 to 4. We also give a geometric characterization of the MLT of a graph in terms of a new "lifting" problem for frameworks that is interesting in its own right. The lifting perspective yields a new connection between the weak MLT (where the maximum likelihood estimate exists only with positive probability) and the classical Hadwiger-Nelson problem.Peer reviewe

    Similar works