Discrepancy of stratified samples from partitions of the unit cube

被引:0
|
作者
Markus Kiderlen
Florian Pausinger
机构
[1] Aarhus University,
[2] Queen’s University Belfast,undefined
来源
Monatshefte für Mathematik | 2021年 / 195卷
关键词
Jittered sampling; Stratified sampling; -discrepancy; Sets of positive reach; Primary 11K38; 60C05; Secondary 05A18; 60D99;
D O I
暂无
中图分类号
学科分类号
摘要
We extend the notion of jittered sampling to arbitrary partitions and study the discrepancy of the related point sets. Let Ω=(Ω1,…,ΩN)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\Omega }}=(\Omega _1,\ldots ,\Omega _N)$$\end{document} be a partition of [0,1]d\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[0,1]^d$$\end{document} and let the ith point in P\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\mathcal {P}}}$$\end{document} be chosen uniformly in the ith set of the partition (and stochastically independent of the other points), i=1,…,N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i=1,\ldots ,N$$\end{document}. For the study of such sets we introduce the concept of a uniformly distributed triangular array and compare this notion to related notions in the literature. We prove that the expected Lp\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal {L}}}_p}$$\end{document}-discrepancy, ELp(PΩ)p\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\mathbb {E}}}{{{\mathcal {L}}}_p}({{\mathcal {P}}}_{\varvec{\Omega }})^p$$\end{document}, of a point set PΩ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\mathcal {P}}}_{\varvec{\Omega }}$$\end{document} generated from any equivolume partition Ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\Omega }}$$\end{document} is always strictly smaller than the expected Lp\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal {L}}}_p}$$\end{document}-discrepancy of a set of N uniform random samples for p>1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p>1$$\end{document}. For fixed N we consider classes of stratified samples based on equivolume partitions of the unit cube into convex sets or into sets with a uniform positive lower bound on their reach. It is shown that these classes contain at least one minimizer of the expected Lp\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal {L}}}_p}$$\end{document}-discrepancy. We illustrate our results with explicit constructions for small N. In addition, we present a family of partitions that seems to improve the expected discrepancy of Monte Carlo sampling by a factor of 2 for every N.
引用
收藏
页码:267 / 306
页数:39
相关论文
共 50 条
  • [31] On entropies for random partitions of the unit segment
    Bieniek, Milena
    Szynal, Dominik
    KYBERNETIKA, 2008, 44 (01) : 75 - 94
  • [32] On some subvarieties of IMTn-algebras and the partitions of the m-cube
    Juan M. Cornejo
    Laura A. Rueda
    Algebra universalis, 2015, 73 : 23 - 42
  • [33] On some subvarieties of IMTn-algebras and the partitions of the m-cube
    Cornejo, Juan M.
    Rueda, Laura A.
    ALGEBRA UNIVERSALIS, 2015, 73 (01) : 23 - 42
  • [34] A chain of controllable partitions of unity on the cube and the approximation of Holder continuous functions
    Richter, C
    ILLINOIS JOURNAL OF MATHEMATICS, 1999, 43 (01) : 159 - 191
  • [35] ROLE OF MASTER POLYTOPES IN UNIT CUBE
    HAMMER, PL
    JOHNSON, EL
    PELED, UN
    NOTICES OF THE AMERICAN MATHEMATICAL SOCIETY, 1975, 22 (01): : A255 - A255
  • [36] ROLE OF MASTER POLYTOPES IN UNIT CUBE
    HAMMER, PL
    JOHNSON, EL
    PELED, UN
    SIAM JOURNAL ON APPLIED MATHEMATICS, 1977, 32 (04) : 711 - 716
  • [37] EULERIAN-NUMBERS AND THE UNIT CUBE
    HENSLEY, D
    FIBONACCI QUARTERLY, 1982, 20 (04): : 344 - 348
  • [38] UNIT CUBE - ITS POSITION IN SPACE
    WRIGHT, IL
    TEXAS JOURNAL OF SCIENCE, 1975, 26 (1-2): : 67 - 91
  • [39] ROBUST EMPIRICAL BAYES ESTIMATION OF VARIANCES FROM STRATIFIED SAMPLES
    GHOSH, M
    LAHIRI, P
    SANKHYA-THE INDIAN JOURNAL OF STATISTICS SERIES B, 1987, 49 : 78 - 89
  • [40] ONLINE COVERING THE UNIT CUBE BY CUBES
    JANUSZEWSKI, J
    LASSAK, M
    DISCRETE & COMPUTATIONAL GEOMETRY, 1994, 12 (04) : 433 - 438