On dropping the first Sobol' point

08/18/2020
by   Art B. Owen, et al.
0

Quasi-Monte Carlo (QMC) points are a substitute for plain Monte Carlo (MC) points that greatly improve integration accuracy under mild assumptions on the problem. Because QMC can give errors that are o(1/n) as n→∞, changing even one point can change the estimate by an amount much larger than the error would have been and worsen the convergence rate. As a result, certain practices that fit quite naturally and intuitively with MC points are very detrimental to QMC performance. These include thinning, burn-in, and taking sample sizes such as powers of 10, other than the ones for which the QMC points were designed. This article looks at the effects of a common practice in which one skips the first point of a Sobol' sequence. The retained points ordinarily fail to be a digital net and when scrambling is applied, skipping over the first point can increase the numerical error by a factor proportional to √(n) where n is the number of function evaluations used.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
11/05/2019

Quasi-Monte Carlo sampling for machine-learning partial differential equations

Solving partial differential equations in high dimensions by deep neural...
research
08/24/2011

Using Supervised Learning to Improve Monte Carlo Integral Estimation

Monte Carlo (MC) techniques are often used to estimate integrals of a mu...
research
10/22/2019

A quasi-Monte Carlo Method for an Optimal Control Problem Under Uncertainty

We study an optimal control problem under uncertainty, where the target ...
research
06/09/2017

A randomized Halton algorithm in R

Randomized quasi-Monte Carlo (RQMC) sampling can bring orders of magnitu...
research
09/21/2023

Quasi-Monte Carlo for 3D Sliced Wasserstein

Monte Carlo (MC) approximation has been used as the standard computation...
research
08/02/2019

Why Simple Quadrature is just as good as Monte Carlo

We motive and calculate Newton-Cotes quadrature integration variance and...
research
03/04/2019

Statistical approach to detection of signals by Monte Carlo singular spectrum analysis: Multiple testing

The statistical approach to detection of a signal in noisy series is con...

Please sign up or login with your details

Forgot password? Click here to reset