Search papers, labs, and topics across Lattice.
This paper introduces a methodology for non-parametrically recovering the causal drift mechanism of sparse multivariate stochastic systems from steady-state observations, particularly in contexts like gene expression analysis where data collection is limited. By assuming a time-homogeneous diffusion process has reached equilibrium, the authors demonstrate that the drift function can be identified under a weak non-explosion criterion. The proposed non-parametric kernel estimator is shown to be consistent, and the authors provide a cross-validation scheme for hyperparameter tuning, with implications for generative diffusion models and low-frequency data sampling.
Non-parametric identification of causal mechanisms from steady-state observations could revolutionize data analysis in fields constrained by experimental limitations.
We consider sparse multivariate stochastic systems that evolve in continuous time according to a causal mechanism and present methodology to recover the system's time-infinitesimal transition mechanism from mere cross-sectional data. This observational paradigm is motivated by applications such as gene expression analysis, where destructive experimental techniques may only allow recording data once over a cell's lifetime. Precisely, we assume the system follows a time-homogeneous diffusion process that has reached an equilibrium distribution at observation time. Further, we assume the causal mechanism is fully described by the diffusion drift, is acyclic, and its causal structure graph is known. In this setting, we prove that the full causal mechanism, i.e., the drift function, can be non-parametrically identified under a weak non-explosion criterion. We derive a non-parametric kernel estimator for this challenging inverse problem and prove its consistency. Moreover, we propose a cross-validation scheme for hyperparameter tuning, illustrate the behavior of our estimator in simulations, and we discuss connections with irreversible generative diffusion models and low-frequency sampled data.