A population Monte Carlo scheme with transformed weights and its application to stochastic kinetic models
From MaRDI portal
Publication:5962748
Abstract: This paper addresses the problem of Monte Carlo approximation of posterior probability distributions. In particular, we have considered a recently proposed technique known as population Monte Carlo (PMC), which is based on an iterative importance sampling approach. An important drawback of this methodology is the degeneracy of the importance weights when the dimension of either the observations or the variables of interest is high. To alleviate this difficulty, we propose a novel method that performs a nonlinear transformation on the importance weights. This operation reduces the weight variation, hence it avoids their degeneracy and increases the efficiency of the importance sampling scheme, specially when drawing from a proposal functions which are poorly adapted to the true posterior. For the sake of illustration, we have applied the proposed algorithm to the estimation of the parameters of a Gaussian mixture model. This is a very simple problem that enables us to clearly show and discuss the main features of the proposed technique. As a practical application, we have also considered the popular (and challenging) problem of estimating the rate parameters of stochastic kinetic models (SKM). SKMs are highly multivariate systems that model molecular interactions in biological and chemical problems. We introduce a particularization of the proposed algorithm to SKMs and present numerical results.
Recommendations
- Bayesian computation methods for inference in stochastic kinetic models
- Efficient \(\mathrm{SMC}^2\) schemes for stochastic kinetic models
- Population Monte Carlo algorithm in high dimensions
- On variance stabilisation in population Monte Carlo by double Rao-Blackwellisation
- Particle Markov chain Monte Carlo for efficient numerical simulation
Cites work
- scientific article; zbMATH DE number 2117879 (Why is no real title available?)
- A Particle Filtering Scheme for Processing Time Series Corrupted by Outliers
- Advanced Lectures on Machine Learning
- Bayesian Inference in Econometric Models Using Monte Carlo Integration
- Curse-of-dimensionality revisited: Collapse of the particle filter in very large scale systems
- Fundamentals of stochastic filtering
- Inference for Lévy-driven stochastic volatility models via adaptive sequential Monte Carlo
- Moment closure based parameter inference of stochastic kinetic models
- On the stability of sequential Monte Carlo methods in high dimensions
- Particle Markov Chain Monte Carlo Methods
- Probability Inequalities for Sums of Bounded Random Variables
- Probability with Martingales
- Sequential Imputations and Bayesian Missing Data Problems
- Sequential Monte Carlo Methods in Practice
- Sequential Monte Carlo Samplers
Cited in
(7)- Nested particle filters for online parameter estimation in discrete-time state-space Markov models
- Parameter inference for a stochastic kinetic model of expanded polyglutamine proteins
- Accelerating inference for stochastic kinetic models
- Implicitly adaptive importance sampling
- Direct statistical inference for finite Markov jump processes via the matrix exponential
- Correlated pseudo-marginal schemes for time-discretised stochastic kinetic models
- Bayesian computation methods for inference in stochastic kinetic models
This page was built for publication: A population Monte Carlo scheme with transformed weights and its application to stochastic kinetic models
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5962748)
