This paper is an improvement upon BayFish, our software for Bayesian inference of transcription dynamics from population snapshots of single‐molecule RNA FISH in single cells. However, Bayesian inference is computationally expensive and cannot tackle more complex models of stochastic gene expression in a reasonable amount of time.
To address this problem, we developed a Poisson Mixture Piecewise Deterministic Markov Switch Rate (PM-PDMSR) algorithm that is 1000-fold more efficient than the state-of-the-art algorithms commonly used to solve the time-dependent dynamics of complex models of stochastic gene expression. With this efficiency, we performed a full‐scale Bayesian analysis including the parameter identification, uncertainty quantification, and most importantly, an estimation of the Bayesian evidence for a given gene expression model (compared to other competing complex models). This was a fruitful collaboration with Yen Ting Lin at Los Alamos National Laboratory.
- Lin YT, Buchler NE. Exact and efficient hybrid Monte Carlo algorithm for accelerated Bayesian inference of gene expression models from snapshots of single-cell transcripts. J. Chem. Physics 151: 024106 (2019).