A rare event approach to high-dimensional Approximate Bayesian computation
Autor: | Richard G. Everitt, Dennis Prangle, Theodore Kypraios |
---|---|
Jazyk: | angličtina |
Rok vydání: | 2018 |
Předmět: |
0301 basic medicine
Statistics and Probability Computer science Monte Carlo method Slice sampling Latent variable Machine learning computer.software_genre 01 natural sciences Theoretical Computer Science Hybrid Monte Carlo 010104 statistics & probability 03 medical and health sciences symbols.namesake 0101 mathematics ABC Markov chain Monte Carlo Sequential Monte Carlo Slice sampling Infectious disease modelling business.industry Markov chain Monte Carlo Statistics::Computation Approximate inference 030104 developmental biology Computational Theory and Mathematics symbols Artificial intelligence Statistics Probability and Uncertainty Approximate Bayesian computation business Particle filter computer Algorithm |
ISSN: | 1573-1375 0960-3174 |
Popis: | Approximate Bayesian computation (ABC) methods permit approximate inference for intractable likelihoods when it is possible to simulate from the model. However they perform poorly for high dimensional data, and in practice must usually be used in conjunction with dimension reduction methods, resulting in a loss of accuracy which is hard to quantify or control. We propose a new ABC method for high dimensional data based on rare event methods which we refer to as RE-ABC. This uses a latent variable representation of the model. For a given parameter value, we estimate the probability of the rare event that the latent variables correspond to data roughly consistent with the observations. This is performed using sequential Monte Carlo and slice sampling to systematically search the space of latent variables. In contrast standard ABC can be viewed as using a more naive Monte Carlo estimate. We use our rare event probability estimator as a likelihood estimate within the pseudo-marginal Metropolis-Hastings algorithm for parameter inference.\ud We provide asymptotics showing that RE-ABC has a lower computational cost for high dimensional data than standard ABC methods. We also illustrate our approach empirically, on a Gaussian distribution and an application in infectious disease modelling. |
Databáze: | OpenAIRE |
Externí odkaz: |