Mirror Likelihood for Modeling Manifold Data with Application on Correlation Time Series
Abstract
Manifold-valued data arise in a wide range of scientific applications. Existing statistical methods typically follow either an intrinsic approach, which characterizes the manifold itself via quantities such as geodesic distance or the Fréchet mean, or an extrinsic approach, which embeds the manifold into an ambient Euclidean space, performs unconstrained statistical estimation, and projects the result back to the manifold. While the latter is conceptually and computationally simple, likelihood-based extrinsic solutions remain scarce, due to the lack of a general recipe for choosing an embedding map and for characterizing its impact on the transformed distribution. In this article, we tackle the problem using an augmented likelihood approach: we view the manifold-valued data as the Fenchel mirror map of a latent variable, which is nearly unconstrained, and model the latent variable with a canonical distribution, such as Gaussian. The key transform, the Fenchel mirror map, is a generalization of the Euclidean projection and can effectively map into the interior of the manifold. Taking a Bayesian approach, the proposed framework enjoys a generative model interpretation and an augmented likelihood free from intractable normalization constants, and it supports principled uncertainty quantification and flexible model-based extensions. We study the theoretical and computational properties of the framework and demonstrate its utility in analyzing neuroimaging data of high-dimensional correlation time series, where each correlation matrix is rank deficient.
Related
- Online adaptive non-intrusive model reduction via manifold interpolation and subspace updates: application to FSI convergence acceleration
- Extension of Likelihood Ratio Analysis Method to Skorokhod M1 Topology (with Application to Poissonian Smooth Change-Point Model)
- Reciprocal-Manifold Annealed KKT Flows for Constrained Optimization: Application to the Nonconvex AC Optimal Power Flow
- Nonparametric Correlation Estimator via Solving Cubic Equations and its Application to Brain Functional Connectivity Analysis
- Modeling Dynamic Correlation Matrices with Shrinkage Priors
BibTeX
@misc{hasenauer2026mirror,
title = {{Mirror Likelihood for Modeling Manifold Data with Application on Correlation Time Series}},
author = {Luis Hasenauer and Arkaprava Roy and Leo Duan},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.36448},
eprint = {2610.36448},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.36448},
}
RIS
TY - UNPB
TI - Mirror Likelihood for Modeling Manifold Data with Application on Correlation Time Series
AU - Luis Hasenauer
AU - Arkaprava Roy
AU - Leo Duan
PY - 2026
DA - 2026/10/27
PB - arXiv
DO - 10.48550/arxiv.2610.36448
UR - https://arxiv.org/abs/2610.36448
AB - Manifold-valued data arise in a wide range of scientific applications. Existing statistical methods typically follow either an intrinsic approach, which characterizes the manifold itself via quantities such as geodesic distance or the Fréchet mean, or an extrinsic approach, which embeds the manifold into an ambient Euclidean space, performs unconstrained statistical estimation, and projects the result back to the manifold. While the latter is conceptually and computationally simple, likelihood-based extrinsic solutions remain scarce, due to the lack of a general recipe for choosing an embedding map and for characterizing its impact on the transformed distribution. In this article, we tackle the problem using an augmented likelihood approach: we view the manifold-valued data as the Fenchel mirror map of a latent variable, which is nearly unconstrained, and model the latent variable with a canonical distribution, such as Gaussian. The key transform, the Fenchel mirror map, is a generalization of the Euclidean projection and can effectively map into the interior of the manifold. Taking a Bayesian approach, the proposed framework enjoys a generative model interpretation and an augmented likelihood free from intractable normalization constants, and it supports principled uncertainty quantification and flexible model-based extensions. We study the theoretical and computational properties of the framework and demonstrate its utility in analyzing neuroimaging data of high-dimensional correlation time series, where each correlation matrix is rank deficient.
ER -