We consider continuous time Markovian processes where populations of individual agents interact stochastically according to kinetic rules. Despite the increasing prominence of such models in fields ranging from biology to smart cities, Bayesian inference for such systems remains challenging, as these are continuous time, discrete state systems with potentially infinite state-space. Here we propose a novel efficient algorithm for joint state/parameter posterior sampling in population Markov Jump processes. We introduce a class of pseudo-marginal sampling algorithms based on a random truncation method which enables a principled treatment of infinite state spaces. Extensive evaluation on a number of benchmark models shows that this approach achieves considerable savings compared to state of the art methods, retaining accuracy and fast convergence. We also present results on a synthetic biology data set showing the potential for practical usefulness of our work.

Unbiased Bayesian inference for population Markov jump processes via random truncations / Georgoulas, A.; Hillston, J.; Sanguinetti, G.. - In: STATISTICS AND COMPUTING. - ISSN 0960-3174. - 27:(2017), pp. 991-1002. [10.1007/s11222-016-9667-9]

Unbiased Bayesian inference for population Markov jump processes via random truncations

Sanguinetti, G.
2017-01-01

Abstract

We consider continuous time Markovian processes where populations of individual agents interact stochastically according to kinetic rules. Despite the increasing prominence of such models in fields ranging from biology to smart cities, Bayesian inference for such systems remains challenging, as these are continuous time, discrete state systems with potentially infinite state-space. Here we propose a novel efficient algorithm for joint state/parameter posterior sampling in population Markov Jump processes. We introduce a class of pseudo-marginal sampling algorithms based on a random truncation method which enables a principled treatment of infinite state spaces. Extensive evaluation on a number of benchmark models shows that this approach achieves considerable savings compared to state of the art methods, retaining accuracy and fast convergence. We also present results on a synthetic biology data set showing the potential for practical usefulness of our work.
2017
27
991
1002
https://arxiv.org/abs/1509.08327
Georgoulas, A.; Hillston, J.; Sanguinetti, G.
File in questo prodotto:
File Dimensione Formato  
Georgoulas2017_Article_UnbiasedBayesianInferenceForPo.pdf

non disponibili

Tipologia: Versione Editoriale (PDF)
Licenza: Non specificato
Dimensione 807.37 kB
Formato Adobe PDF
807.37 kB Adobe PDF   Visualizza/Apri   Richiedi una copia

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/20.500.11767/117262
Citazioni
  • ???jsp.display-item.citation.pmc??? 3
  • Scopus 18
  • ???jsp.display-item.citation.isi??? 13
social impact