跳到正文
研踪PaperTrail
← 研踪 PaperTrail

arXiv:2605.18316 · v2

Dynamic Elliptical Graph Factor Models via Riemannian Optimization with Geodesic Temporal Regularization

首次提交 2026-05-18 20:34
版本更新 2026-11-09 21:12
雷达观测 2026-11-10 15:22
arXiv ↗PDF ↗

摘要

Inferring time-varying graphs from high-dimensional nodal observations is challenging when windows contain fewer samples than nodes, observations are heavy-tailed, and the dependence drifts multiplicatively over time. We propose the dynamic elliptical graph factor model (\textsc{Degfm}), combining: an elliptical likelihood that down-weights outliers; a low-rank-plus-diagonal (LRaD) parameterization of each precision matrix, shown to be the posterior precision of a field conditioned on a few aggregate linear measurements; and a temporal regularizer given by the squared affine-invariant distance between consecutive precision matrices, capturing multiplicative evolution and invariant to node rescaling and to the choice between precision and covariance. Estimation on a product of quotient manifolds is solved by a safeguarded Riemannian conjugate-gradient method with closed-form temporal gradients via generalized eigendecompositions. We prove that (i) LRaD is locally identifiable modulo orthogonal transformations under a tangent-space transversality (coherence) condition; (ii) with rank-preserving Armijo backtracking, the solver is well defined and converges to critical points with $\min_{k<K}\|\grad F(θ^{(k)})\|^2=\mathcal{O}(1/K)$; (iii) a localized estimator obeys a finite-sample oracle inequality whose error scales with the intrinsic dimension pr+p−r(r−1)/2pr+p-r(r-1)/2 times an effective temporal dimension decreasing from 11 to 1/T1/T as coupling strengthens, quantifying the variance reduction and bias of temporal coupling, and extending to misspecified non-LRaD models; (iv) for Student-tt observations, the local restricted strong convexity is proved rather than assumed, with curvature of order ν/(p+ν)ν/(p+ν) shown unimprovable; and (v) partial-correlation graphs are recovered exactly by thresholding. Experiments on synthetic data, S\&P~500 returns, and HCP fMRI data validate the approach.

AI 论文解读

基于摘要 · 非全文审阅

此版本尚未生成 AI 解读。解读由后台任务生成并按论文版本保存,登录后共享已有结果。

版本差异

本站首次收录为 v2,暂无更早版本可对照。

观测到的版本

相邻记录可能跳过中间版本;首次发现 v2 及以上也会标为版本更新。

v2 · 2026-11-09 21:12 · 最新观测版本

Dynamic Elliptical Graph Factor Models via Riemannian Optimization with Geodesic Temporal Regularization

Inferring time-varying graphs from high-dimensional nodal observations is challenging when windows contain fewer samples than nodes, observations are heavy-tailed, and the dependence drifts multiplicatively over time. We propose the dynamic elliptical graph factor model (\textsc{Degfm}), combining: an elliptical likelihood that down-weights outliers; a low-rank-plus-diagonal (LRaD) parameterization of each precision matrix, shown to be the posterior precision of a field conditioned on a few aggregate linear measurements; and a temporal regularizer given by the squared affine-invariant distance between consecutive precision matrices, capturing multiplicative evolution and invariant to node rescaling and to the choice between precision and covariance. Estimation on a product of quotient manifolds is solved by a safeguarded Riemannian conjugate-gradient method with closed-form temporal gradients via generalized eigendecompositions. We prove that (i) LRaD is locally identifiable modulo orthogonal transformations under a tangent-space transversality (coherence) condition; (ii) with rank-preserving Armijo backtracking, the solver is well defined and converges to critical points with $\min_{k<K}\|\grad F(θ^{(k)})\|^2=\mathcal{O}(1/K)$; (iii) a localized estimator obeys a finite-sample oracle inequality whose error scales with the intrinsic dimension pr+p−r(r−1)/2pr+p-r(r-1)/2 times an effective temporal dimension decreasing from 11 to 1/T1/T as coupling strengthens, quantifying the variance reduction and bias of temporal coupling, and extending to misspecified non-LRaD models; (iv) for Student-tt observations, the local restricted strong convexity is proved rather than assumed, with curvature of order ν/(p+ν)ν/(p+ν) shown unimprovable; and (v) partial-correlation graphs are recovered exactly by thresholding. Experiments on synthetic data, S\&P~500 returns, and HCP fMRI data validate the approach.

查看此版本解读 →查看此版本原文 ↗