Walsh functions, scrambled (0,m,s)-nets, and negative covariance: applying symbolic computation to quasi-Monte Carlo integration

06/11/2020
by   Jaspar Wiart, et al.
0

We investigate base b Walsh functions for which the variance of the integral estimator based on a scrambled (0,m,s)-net in base b is less than or equal to that of the Monte-Carlo estimator based on the same number of points. First we compute the Walsh decomposition for the joint probability density function of two distinct points randomly chosen from a scrambled (t,m,s)-net in base b in terms of certain counting numbers and simplify it in the special case t is zero. Using this, we obtain an expression for the covariance of the integral estimator in terms of the Walsh coefficients of the function. Finally, we prove that the covariance of the integral estimator is negative when the Walsh coefficients of the function satisfy a certain decay condition. To do this, we use creative telescoping and recurrence solving algorithms from symbolic computation to find a sign equivalent closed form expression for the covariance term.

READ FULL TEXT
research
07/11/2022

Improved bounds on the gain coefficients for digital nets in prime base

We study randomized quasi-Monte Carlo integration by scrambled nets. The...
research
06/23/2021

A note on concatenation of quasi-Monte Carlo and plain Monte Carlo rules in high dimensions

In this short note, we study a concatenation of quasi-Monte Carlo and pl...
research
03/09/2020

Toeplitz Monte Carlo

Motivated mainly by applications to partial differential equations with ...
research
06/19/2021

The nonzero gain coefficients of Sobol's sequences are always powers of two

When a plain Monte Carlo estimate on n samples has variance σ^2/n, then ...
research
11/14/2022

Regression-based Monte Carlo Integration

Monte Carlo integration is typically interpreted as an estimator of the ...
research
07/14/2020

Designed Quadrature to Approximate Integrals in Maximum Simulated Likelihood Estimation

Maximum simulated likelihood estimation of mixed multinomial logit (MMNL...

Please sign up or login with your details

Forgot password? Click here to reset