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 | |
|---|---|---|---|---|
| Population Growth Modeling | 50D CITE t=2 | W1 Score27.908 | 12 | |
| Population Growth Modeling | 50D CITE t=3 | W132.95 | 12 | |
| Population Growth Modeling | 50D CITE t=1 | W128.245 | 12 | |
| Trajectory Inference | 100D EB dataset (t=4) | W1 Error12.406 | 11 | |
| Trajectory Inference | 100D EB dataset t=2 | W111.103 | 11 | |
| Trajectory Inference | 100D EB dataset (t=3) | W1 Error11.529 | 11 | |
| Cell population dynamics prediction | 50D Mouse (t=1) | W1 Score6.052 | 11 | |
| Trajectory Inference | 100D EB dataset t=1 | W110.256 | 11 | |
| Trajectory Inference | EB 50D t=4 | W1 Score9.733 | 11 | |
| Trajectory Inference | Dyngen (t=2) | W10.481 | 11 |