Fundamental limits of distributed covariance matrix estimation via a conditional strong data processing inequality

TL;DR

This paper establishes fundamental limits and nearly optimal protocols for distributed covariance matrix estimation under communication constraints, introducing a novel conditional SDPI to analyze the problem's information-theoretic bounds.

Contribution

It introduces the Conditional Strong Data Processing Inequality (C-SDPI), a new tool for analyzing distributed estimation limits, and derives tight lower bounds and matching protocols without assuming infinite samples or Gaussianity.

Findings

Derived nearly tight minimax lower bounds for covariance estimation

Proposed a nearly optimal estimation protocol matching bounds

Showed interaction reduces communication costs significantly

Abstract

Estimating high-dimensional covariance matrices is a key task across many fields. This paper explores the theoretical limits of distributed covariance estimation in a feature-split setting, where communication between agents is constrained. Specifically, we study a scenario in which multiple agents each observe different components of i.i.d. samples drawn from a sub-Gaussian random vector. A central server seeks to estimate the complete covariance matrix using a limited number of bits communicated by each agent. We obtain a nearly tight minimax lower bound for covariance matrix estimation under operator norm and Frobenius norm. Our main technical tool is a novel generalization of the strong data processing inequality (SDPI), termed the Conditional Strong Data Processing Inequality (C-SDPI) coefficient, introduced in this work. The C-SDPI coefficient shares key properties such as tensorization with the conventional SDPI. Crucially, it quantifies the average contraction in a state-dependent channel and can be significantly lower than the worst-case SDPI coefficient over the state input. Utilizing the doubling trick of Geng-Nair and an operator Jensen inequality, we compute this coefficient for Gaussian mixture channels. We then employ it to establish minimax lower bounds on estimation error, capturing the trade-offs among sample size, communication cost, and data dimensionality. Building on this, we present a nearly optimal estimation protocol whose sample and communication requirements match the lower bounds up to logarithmic factors. Unlike much of the existing literature, our framework does not assume infinite samples or Gaussian distributions, making it broadly applicable. Finally, we extend our analysis to interactive protocols, showing interaction can significantly reduce communication requirements compared to non-interactive schemes.