904 research outputs found
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
Higher order scrambled digital nets achieve the optimal rate of the root mean square error for smooth integrands
We study a random sampling technique to approximate integrals
by averaging the function
at some sampling points. We focus on cases where the integrand is smooth, which
is a problem which occurs in statistics. The convergence rate of the
approximation error depends on the smoothness of the function and the
sampling technique. For instance, Monte Carlo (MC) sampling yields a
convergence of the root mean square error (RMSE) of order (where
is the number of samples) for functions with finite variance. Randomized
QMC (RQMC), a combination of MC and quasi-Monte Carlo (QMC), achieves a RMSE of
order under the stronger assumption that the integrand
has bounded variation. A combination of RQMC with local antithetic sampling
achieves a convergence of the RMSE of order (where
is the dimension) for functions with mixed partial derivatives up to
order two. Additional smoothness of the integrand does not improve the rate of
convergence of these algorithms in general. On the other hand, it is known that
without additional smoothness of the integrand it is not possible to improve
the convergence rate. This paper introduces a new RQMC algorithm, for which we
prove that it achieves a convergence of the root mean square error (RMSE) of
order provided the integrand satisfies the strong
assumption that it has square integrable partial mixed derivatives up to order
in each variable. Known lower bounds on the RMSE show that this rate
of convergence cannot be improved in general for integrands with this
smoothness. We provide numerical examples for which the RMSE converges
approximately with order and , in accordance with the
theoretical upper bound.Comment: Published in at http://dx.doi.org/10.1214/11-AOS880 the Annals of
Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical
Statistics (http://www.imstat.org
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
Efficient calculation of the worst-case error and (fast) component-by-component construction of higher order polynomial lattice rules
We show how to obtain a fast component-by-component construction algorithm
for higher order polynomial lattice rules. Such rules are useful for
multivariate quadrature of high-dimensional smooth functions over the unit cube
as they achieve the near optimal order of convergence. The main problem
addressed in this paper is to find an efficient way of computing the worst-case
error. A general algorithm is presented and explicit expressions for base~2 are
given. To obtain an efficient component-by-component construction algorithm we
exploit the structure of the underlying cyclic group.
We compare our new higher order multivariate quadrature rules to existing
quadrature rules based on higher order digital nets by computing their
worst-case error. These numerical results show that the higher order polynomial
lattice rules improve upon the known constructions of quasi-Monte Carlo rules
based on higher order digital nets
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
- …