45 research outputs found
Application of Sequential Quasi-Monte Carlo to Autonomous Positioning
Sequential Monte Carlo algorithms (also known as particle filters) are
popular methods to approximate filtering (and related) distributions of
state-space models. However, they converge at the slow rate, which
may be an issue in real-time data-intensive scenarios. We give a brief outline
of SQMC (Sequential Quasi-Monte Carlo), a variant of SMC based on
low-discrepancy point sets proposed by Gerber and Chopin (2015), which
converges at a faster rate, and we illustrate the greater performance of SQMC
on autonomous positioning problems.Comment: 5 pages, 4 figure
Improved bounds on the gain coefficients for digital nets in prime power base
We study randomized quasi-Monte Carlo integration by scrambled nets. The
scrambled net quadrature has long gained its popularity because it is an
unbiased estimator of the true integral, allows for a practical error
estimation, achieves a high order decay of the variance for smooth functions,
and works even for -functions with any . The variance of the
scrambled net quadrature for -functions can be evaluated through the set
of the so-called gain coefficients.
In this paper, based on the system of Walsh functions and the concept of dual
nets, we provide improved upper bounds on the gain coefficients for digital
nets in general prime power base. Our results explain the known bound by Owen
(1997) for Faure sequences, the recently improved bound by Pan and Owen (2021)
for digital nets in base 2 (including Sobol' sequences as a special case), and
their finding that all the nonzero gain coefficients for digital nets in base 2
must be powers of two, all in a unified way.Comment: minor revision, 14 page
The Discrepancy and Gain Coefficients of Scrambled Digital Nets
AbstractDigital sequences and nets are among the most popular kinds of low discrepancy sequences and sets and are often used for quasi-Monte Carlo quadrature rules. Several years ago Owen proposed a method of scrambling digital sequences and recently Faure and Tezuka have proposed another method. This article considers the discrepancy of digital nets under these scramblings. The first main result of this article is a formula for the discrepancy of a scrambled digital (λ, t, m, s)-net in base b with n=λbm points that requires only O(n) operations to evaluate. The second main result is exact formulas for the gain coefficients of a digital (t, m, s)-net in terms of its generator matrices. The gain coefficients, as defined by Owen, determine both the worst-case and random-case analyses of quadrature error
On Integration Methods Based on Scrambled Nets of Arbitrary Size
We consider the problem of evaluating for a function . In situations where
can be approximated by an estimate of the form
, with a point set in
, it is now well known that the Monte Carlo
convergence rate can be improved by taking for the first
points, , of a scrambled
-sequence in base . In this paper we derive a bound for the
variance of scrambled net quadrature rules which is of order
without any restriction on . As a corollary, this bound allows us to provide
simple conditions to get, for any pattern of , an integration error of size
for functions that depend on the quadrature size . Notably,
we establish that sequential quasi-Monte Carlo (M. Gerber and N. Chopin, 2015,
\emph{J. R. Statist. Soc. B, to appear.}) reaches the
convergence rate for any values of . In a numerical study, we show that for
scrambled net quadrature rules we can relax the constraint on without any
loss of efficiency when the integrand is a discontinuous function
while, for sequential quasi-Monte Carlo, taking may only
provide moderate gains.Comment: 27 pages, 2 figures (final version, to appear in The Journal of
Complexity
Local antithetic sampling with scrambled nets
We consider the problem of computing an approximation to the integral
. Monte Carlo (MC) sampling typically attains a root
mean squared error (RMSE) of from independent random function
evaluations. By contrast, quasi-Monte Carlo (QMC) sampling using carefully
equispaced evaluation points can attain the rate for
any and randomized QMC (RQMC) can attain the RMSE
, both under mild conditions on . Classical
variance reduction methods for MC can be adapted to QMC. Published results
combining QMC with importance sampling and with control variates have found
worthwhile improvements, but no change in the error rate. This paper extends
the classical variance reduction method of antithetic sampling and combines it
with RQMC. One such method is shown to bring a modest improvement in the RMSE
rate, attaining for any , for
smooth enough .Comment: Published in at http://dx.doi.org/10.1214/07-AOS548 the Annals of
Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical
Statistics (http://www.imstat.org