Abstract
We derive and study sequential quasi Monte Carlo (SQMC), a class of algorithms
obtained by introducing QMC point sets in particle filtering. SQMC is related to, and may be seen as an extension of, the array-RQMC algorithm of L’Ecuyer and his colleagues. The complexity of SQMC is O(N logN) , where N is the number of simulations at each iteration, and its error rate is smaller than the Monte Carlo rate O_P(N^{1/2}). The only requirement to implement SQMC algorithms is the ability to write the simulation of particle x_t^n given x_{t-1}^n as a deterministic
function of x_{t-1}^n and a fixed number of uniform variates. We show that SQMC is amenable to the same extensions as standard SMC, such as forward smoothing, backward smoothing and unbiased likelihood evaluation. In particular, SQMC may replace SMC within a particle Markov chain Monte Carlo algorithm. We establish several convergence results. We provide numerical
evidence that SQMC may significantly outperform SMC in practical scenarios
obtained by introducing QMC point sets in particle filtering. SQMC is related to, and may be seen as an extension of, the array-RQMC algorithm of L’Ecuyer and his colleagues. The complexity of SQMC is O(N logN) , where N is the number of simulations at each iteration, and its error rate is smaller than the Monte Carlo rate O_P(N^{1/2}). The only requirement to implement SQMC algorithms is the ability to write the simulation of particle x_t^n given x_{t-1}^n as a deterministic
function of x_{t-1}^n and a fixed number of uniform variates. We show that SQMC is amenable to the same extensions as standard SMC, such as forward smoothing, backward smoothing and unbiased likelihood evaluation. In particular, SQMC may replace SMC within a particle Markov chain Monte Carlo algorithm. We establish several convergence results. We provide numerical
evidence that SQMC may significantly outperform SMC in practical scenarios
Original language | English |
---|---|
Pages (from-to) | 509-579 |
Number of pages | 71 |
Journal | Journal of the Royal Statistical Society: Series B |
Volume | 77 |
Issue number | 3 |
Early online date | 12 May 2015 |
DOIs | |
Publication status | Published - 1 Jun 2015 |
Keywords
- Array-randomized quasi Monte Carlo
- Low discrepancy
- Particle filtering
- quasi-Monte Carlo
- Randomized quasi Monte Carlo
- Sequential Monte Carlo
Fingerprint
Dive into the research topics of 'Sequential quasi Monte Carlo'. Together they form a unique fingerprint.Profiles
-
Dr Mathieu Gerber
- School of Mathematics - Senior Lecturer in Statistical Science
- Statistical Science
Person: Academic , Member