Optimal Transport for Stationary Markov Chains via Policy Iteration
arXiv:2006.07998
Abstract
We study the optimal transport problem for pairs of stationary finite-state Markov chains, with an emphasis on the computation of optimal transition couplings. Transition couplings are a constrained family of transport plans that capture the dynamics of Markov chains. Solutions of the optimal transition coupling (OTC) problem correspond to alignments of the two chains that minimize long-term average cost. We establish a connection between the OTC problem and Markov decision processes, and show that solutions of the OTC problem can be obtained via an adaptation of policy iteration. For settings with large state spaces, we develop a fast approximate algorithm based on an entropy-regularized version of the OTC problem, and provide bounds on its per-iteration complexity. We establish a stability result for both the regularized and unregularized algorithms, from which a statistical consistency result follows as a corollary. We validate our theoretical results empirically through a simulation study, demonstrating that the approximate algorithm exhibits faster overall runtime with low error. Finally, we extend the setting and application of our methods to hidden Markov models, and illustrate the potential use of the proposed algorithms in practice with an application to computer-generated music.
References in corpus (9)
- Learning with a Wasserstein Loss
- Near-linear time approximation algorithms for optimal transport via Sinkhorn iteration
- Differential Properties of Sinkhorn Approximation for Learning with Wasserstein Distance
- On Efficient Optimal Transport: An Analysis of Greedy and Accelerated Mirror Descent Algorithms
- Statistical Optimal Transport via Factored Couplings
- Distilled Wasserstein Learning for Word Embedding and Topic Modeling
- Fast Algorithms for Computational Optimal Transport and Wasserstein Barycenter
- Spatio-Temporal Alignments: Optimal transport through space and time
- Classical Music Composition Using State Space Models