Decentralized Inference for Spatial Data Using Low-Rank Models
Jianwei Shi, Sameh Abdulah, Ying Sun, Marc G. Genton
TL;DR
This work tackles scalable decentralized parameter inference for spatial low-rank models by replacing the non-decomposable log-likelihood with a variational ELBO objective. The ELBO decomposes into local and common components, enabling decentralized optimization via a block coordinate descent framework that leverages dynamic consensus averaging and multi-consensus. The authors establish local convexity of the negative ELBO near the true parameters and prove consistency and asymptotic normality of the estimators, along with convergence guarantees for the algorithm. Through extensive simulations and a Total Precipitable Water data application, the method demonstrates robustness to data partitioning, network structure, and varying ranks, while delivering competitive predictive performance and substantial computational advantages in distributed settings.
Abstract
Advancements in information technology have enabled the creation of massive spatial datasets, driving the need for scalable and efficient computational methodologies. While offering viable solutions, centralized frameworks are limited by vulnerabilities such as single-point failures and communication bottlenecks. This paper presents a decentralized framework tailored for parameter inference in spatial low-rank models to address these challenges. A key obstacle arises from the spatial dependence among observations, which prevents the log-likelihood from being expressed as a summation-a critical requirement for decentralized optimization approaches. To overcome this challenge, we propose a novel objective function leveraging the evidence lower bound, which facilitates the use of decentralized optimization techniques. Our approach employs a block descent method integrated with multi-consensus and dynamic consensus averaging for effective parameter optimization. We prove the convexity of the new objective function in the vicinity of the true parameters, ensuring the convergence of the proposed method. Additionally, we present the first theoretical results establishing the consistency and asymptotic normality of the estimator within the context of spatial low-rank models. Extensive simulations and real-world data experiments corroborate these theoretical findings, showcasing the robustness and scalability of the framework.
