Hostname: page-component-848d4c4894-v5vhk Total loading time: 0 Render date: 2024-06-15T11:34:59.866Z Has data issue: false hasContentIssue false

Simulation from quasi-stationary distributions on reducible state spaces

Published online by Cambridge University Press:  08 September 2017

A. Griffin*
University of Warwick
P. A. Jenkins*
University of Warwick
G. O. Roberts*
University of Warwick
S. E. F. Spencer*
University of Warwick
* Postal address: Department of Statistics, University of Warwick, Coventry CV4 7AL, UK.
* Postal address: Department of Statistics, University of Warwick, Coventry CV4 7AL, UK.
* Postal address: Department of Statistics, University of Warwick, Coventry CV4 7AL, UK.
* Postal address: Department of Statistics, University of Warwick, Coventry CV4 7AL, UK.


Quasi-stationary distributions (QSDs) arise from stochastic processes that exhibit transient equilibrium behaviour on the way to absorption. QSDs are often mathematically intractable and even drawing samples from them is not straightforward. In this paper the framework of sequential Monte Carlo samplers is utilised to simulate QSDs and several novel resampling techniques are proposed to accommodate models with reducible state spaces, with particular focus on preserving particle diversity on discrete spaces. Finally, an approach is considered to estimate eigenvalues associated with QSDs, such as the decay parameter.

Research Article
Copyright © Applied Probability Trust 2017 

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)


[1] Blanchet, J., Glynn, P. and Zheng, S. (2016). Analysis of a stochastic approximation algorithm for computing quasi-stationary distributions. Adv. Appl. Prob. 48, 792811. Google Scholar
[2] Clancy, D. and Pollett, P. K. (2003). A note on quasi-stationary distributions of birth–death processes and the SIS logistic epidemic. J. Appl. Prob. 40, 821825. CrossRefGoogle Scholar
[3] Darroch, J. N. and Seneta, E. (1967). On quasi-stationary distributions in absorbing continuous-time finite Markov chains. J. Appl. Prob. 4, 192196. CrossRefGoogle Scholar
[4] Del Moral, P. (2004). Feynman-Kac Formulae: Genealogical and Interacting Particle Systems with Applications. Springer, New York. Google Scholar
[5] Del Moral, P., Doucet, A. and Jasra, A. (2006). Sequential Monte Carlo samplers. J. R. Statist. Soc. B 68, 411436. Google Scholar
[6] Douc, R. and Moulines, E. (2008). Limit theorems for weighted samples with applications to sequential Monte Carlo methods. Ann. Statist. 36, 23442376. Google Scholar
[7] Doucet, A., de Freitas, N. and Gordon, N. (eds) (2001). Sequential Monte Carlo Methods in Practice. Springer, New York. Google Scholar
[8] Etheridge, A. (2011). Some Mathematical Models from Population Genetics. Springer, Heidelberg. Google Scholar
[9] Griffin, A. (2016). Quasi-stationary distributions for epidemic models: simulation and characterisation. Doctoral thesis, University of Warwick. Google Scholar
[10] Groisman, P. and Jonckheere, M. (2013). Simulation of quasi-stationary distributions on countable spaces. Markov Process. Relat. Fields 19, 521542. Google Scholar
[11] Hairer, M., Stuart, A. M. and Vollmer, S. J. (2014). Spectral gaps for a Metropolis–Hastings algorithm in infinite dimensions. Ann. Appl. Prob. 24, 24552490. Google Scholar
[12] Jenkins, P. A. (2012). Stopping-time resampling and population genetic inference under coalescent models. Stat. Appl. Genet. Molec. Biol. 11, 9. Google Scholar
[13] Kong, A., Liu, J. S. and Wong, W. H. (1994). Sequential imputations and Bayesian missing data problems. J. Amer. Statist. Assoc. 89, 278288. Google Scholar
[14] Lambert, A. (2008). Population dynamics and random genealogies. Stoch. Models 24, 45163. Google Scholar
[15] Lin, M., Chen, R. and Liu, J. S. (2013). Lookahead strategies for sequential Monte Carlo. Statist. Sci. 28, 6994. CrossRefGoogle Scholar
[16] Liu, J. S. and Chen, R. (1995). Blind deconvolution via sequential imputations. J. Amer. Statist. Assoc. 90, 567576. Google Scholar
[17] Liu, J. S. and Chen, R. (1998). Sequential Monte Carlo methods for dynamic systems. J. Amer. Statist. Assoc. 93, 10321044. Google Scholar
[18] Méléard, S. and Villemonais, D. (2012). Quasi-stationary distributions and population processes. Prob. Surveys 9, 340410. Google Scholar
[19] Nåsell, I. (1999). On the quasi-stationary distribution of the stochastic logistic epidemic. Math. Biosci. 156, 2140. Google Scholar
[20] Neal, P. (2014). Endemic behaviour of SIS epidemics with general infectious period distributions. Adv. Appl. Prob. 46, 241255. Google Scholar
[21] Pollock, M., Fearnhead, P., Johansen, A. M. and Roberts, G. O. (2016). The scalable Langevin exact algorithm: Bayesian inference for big data. Preprint. Available at Google Scholar
[22] Van Doorn, E. A. (2015). Representations for the decay parameter of a birth–death process based on the Courant–Fischer theorem. J. Appl. Prob. 52, 278289. Google Scholar
[23] Van Doorn, E. A. and Pollett, P. K. (2008). Survival in a quasi-death process. Linear Algebra Appl. 429, 776791. Google Scholar
[24] Van Doorn, E. A. and Pollett, P. K. (2013). Quasi-stationary distributions for discrete-state models. Europ. J. Operat. Res. 230, 114. Google Scholar
[25] Wright, S. (1931). Evolution in Mendelian populations. Genetics 16, 97159. CrossRefGoogle ScholarPubMed