Summary
The authors provide a framework that allows training of GFlowNet models with streaming data by checkpointing a previously trained GFlowNet model and proposing the streaming balance condition. The problem as well as the approach is well motivated and theoretically sound and this work would be a valuable contribution to the research community. The authors conduct experiments on synthetic tasks and phylogenetic inference, and show improved scalability. Additionally, there is some theory that describes and bounds the approximation error based on errors on the checkpointed GFlowNet model and current estimation problem. Overall, I think the work is quite interesting but it would be nice to consider another experiment since currently most of the experiments are based on synthetic data.
Strengths
- The authors provide useful theoretical contributions that bound the approximation errors into two terms: error from solving the current distribution matching problem and error coming from suboptimal solution to the previous distribution matching problem.
- The problem of training sampling methods that can incorporate streaming data is quite important and relevant to the field, especially for Bayesian Inference methodologies.
- The authors do provide a good set of experiments, and show the benefits of their proposed approach over the standard training of GFlowNets.
Weaknesses
- I found Definition 2 to be quite confusing to read and it is still not clear to me how the equivalence holds. Could the authors clarify what the KL divergence is between, why is it valid and how does it circumvent the problem of parameterizing a partition function?
- The original setup proposed by the authors follows from a trivial extension of the GFlowNets framework which begs to question the novelty of the work.
- It is unclear how their formulation solves the problem of permutation invariance / iid treatment of observations. In particular, given two sets $D_1$ and $D_2$, we know that the posterior distribution $p(x | D_1, D_2)$ can be equivalently written as
$$
p(x | D_1, D_2) = \frac{p(D_1, D_2 | x) p(x)}{p(D_1, D_2)} = \frac{p(D_1 | x) p(D_2 | x) p(x)}{p(D_1, D_2)} = \frac{p(D_1 | x) p(x|D_2)}{p(D_1 | D_2)} = \frac{p(D_2 | x) p(x|D_1)}{p(D_2 | D_1)}
$$
Then, in particular, how is it maintained that given any split of the data $D$ into $D_1$ and $D_2$, the final GFN trained leads to the same solution. If the GFN learns the right solution, then it will follow simply, but it would be nice if the model could satisfy this permutation invariance inherently. Could the authors provide some ablations and control experiments to show to what extent is this constraint satisfied over $t$?
- The authors should add an experiment on discovering causal graphs conditioned on observational samples, which would strengthen their paper considerably.
- The authors should consider the following related work, which are either connected to GFlowNets, sampling from an unnormalized density, or modeling Bayesian posterior inference:
Cranmer, Kyle, Johann Brehmer, and Gilles Louppe. "The frontier of simulation-based inference." Proceedings of the National Academy of Sciences 117.48 (2020): 30055-30062.
Mittal, Sarthak, et al. "Exploring Exchangeable Dataset Amortization for Bayesian Posterior Inference." ICML 2023 Workshop on Structured Probabilistic Inference {\&} Generative Modeling. 2023.
Zhang, Qinsheng, and Yongxin Chen. "Path integral sampler: a stochastic control approach for sampling." arXiv preprint arXiv:2111.15141 (2021).
Richter, Lorenz, et al. "VarGrad: a low-variance gradient estimator for variational inference." Advances in Neural Information Processing Systems 33 (2020): 13481-13492.
Sendera, Marcin, et al. "On diffusion models for amortized inference: Benchmarking and improving stochastic control and sampling." arXiv preprint arXiv:2402.05098 (2024).
Berner, Julius, Lorenz Richter, and Karen Ullrich. "An optimal control perspective on diffusion-based generative modeling." arXiv preprint arXiv:2211.01364 (2022).
Akhound-Sadegh, Tara, et al. "Iterated denoising energy matching for sampling from Boltzmann densities." arXiv preprint arXiv:2402.06121 (2024).
Richter, Lorenz, Julius Berner, and Guan-Horng Liu. "Improved sampling via learned diffusions." arXiv preprint arXiv:2307.01198 (2023).
Questions
- It is not clear what do the authors mean by $\mathbf{A}$ in their notation at the start of Section 2.
- The authors seem to have a typo in the equation under Line 142, essentially the right hand equation should have $p_F^{t+1}$.
- Can the authors provide a visualization of the Set generation task? If I understand correctly, the authors first sample $d$ points randomly and the corresponding $f_i(j)$ randomly as well where $i = 1, ..., K+1$ and $j = 1, ..., d$. Then they define $R_i$ over each set, and a tempered version of $R_i$ with temperature $\alpha$. Is that it?