Learning stochastic dynamics from snapshots through regularized unbalanced optimal transport
About
Reconstructing dynamics using samples from sparsely time-resolved snapshots is an important problem in both natural sciences and machine learning. Here, we introduce a new deep learning approach for solving regularized unbalanced optimal transport (RUOT) and inferring continuous unbalanced stochastic dynamics from observed snapshots. Based on the RUOT form, our method models these dynamics without requiring prior knowledge of growth and death processes or additional information, allowing them to be learned directly from data. Theoretically, we explore the connections between the RUOT and Schr\"odinger bridge problem and discuss the key challenges and potential solutions. The effectiveness of our method is demonstrated with a synthetic gene regulatory network, high-dimensional Gaussian Mixture Model, and single-cell RNA-seq data from blood development. Compared with other methods, our approach accurately identifies growth and transition patterns, eliminates false transitions, and constructs the Waddington developmental landscape. Our code is available at: https://github.com/zhenyiizhang/DeepRUOT.
Related benchmarks
| Task | Dataset | Result | Rank | |
|---|---|---|---|---|
| Trajectory Interpolation | EB 50D (held-out time points) | Mean W16.847 | 24 | |
| Reconstruction of discrete branching dynamics | CITE 50D (t=3) | W19.378 | 24 | |
| Reconstruction of discrete branching dynamics | 5D EB t=1 | W10.386 | 24 | |
| Trajectory Inference | EB dataset 5D (test) | W1 (t=1)0.386 | 23 | |
| Trajectory reconstruction | Gaussian Mixtures 1000D | W1 Distance3.785 | 18 | |
| Single-cell interpolation | EB 5D | EMD0.776 | 15 | |
| Discrete Branching Dynamics Reconstruction | Dyngen (t=2) | W1 Distance0.481 | 15 | |
| Discrete Branching Dynamics Reconstruction | Dyngen t=3 | W1 Distance0.87 | 15 | |
| Discrete Branching Dynamics Reconstruction | Dyngen t=4 | W1 Distance0.688 | 15 | |
| Discrete Branching Dynamics Reconstruction | Dyngen t=1 | W1 Distance0.454 | 15 |