On the discrepancy of jittered sampling

From MaRDI portal



Abstract: We study the discrepancy of jittered sampling sets: such a set mathcalPsubset[0,1]d is generated for fixed minmathbbN by partitioning [0,1]d into md axis aligned cubes of equal measure and placing a random point inside each of the N=md cubes. We prove that, for N sufficiently large, frac{1}{10}frac{d}{N^{frac{1}{2} + frac{1}{2d}}} leq mathbb{E} D_N^*(mathcal{P}) leq frac{sqrt{d} (log{N})^{frac{1}{2}}}{N^{frac{1}{2} + frac{1}{2d}}}, where the upper bound with an unspecified constant Cd was proven earlier by Beck. Our proof makes crucial use of the sharp Dvoretzky-Kiefer-Wolfowitz inequality and a suitably taylored Bernstein inequality; we have reasons to believe that the upper bound has the sharp scaling in N. Additional heuristics suggest that jittered sampling should be able to improve known bounds on the inverse of the star-discrepancy in the regime Ngtrsimdd. We also prove a partition principle showing that every partition of [0,1]d combined with a jittered sampling construction gives rise to a set whose expected squared L2−discrepancy is smaller than that of purely random points.




Cites work









This page was built for publication: On the discrepancy of jittered sampling

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5963459)