Erratum to "Higher order scrambled digital nets achieve the optimal rate of the root mean square error for smooth integrands"

Technische Samenvatting van de Erratum

Probleemstelling
Deze nota behandelt kritieke fouten in het oorspronkelijke artikel "Higher order scrambled digital nets achieve the optimal rate of the root mean square error for smooth integrands" (J. Dick, Annals of Statistics 39, 2011). Het oorspronkelijke werk beweerde dat orde-$d$ geneste uniform gescrambelde digitale $(t, m, s)$-nets een wortelgetal van de kwadratische gemiddelde fout (RMSE) bereiken van $O(N^{-\min(\alpha, d) - 1/2 + \varepsilon})$ voor functies in ongeankerde Sobolev-ruimten met gemengde partiële afgeleiden tot orde $\alpha$. Hoewel het hoofdresultaat met betrekking tot de optimale algebraïsche snelheid geldig blijft, bevatte het oorspronkelijke artikel fundamentele gebreken in de bewijsmachinerie. Specifiek was de geclaimde equivalentie tussen een voorgestelde eindverschil-variatie en de Sobolev-norm incorrect, de covariantiestructuur die werd gebruikt in de variantie-decompositie vertrouwde op een slecht gedefinieerde inverse interlacing-afbeelding, en de variantie-grens in Stelling 10 bevatte een ontbrekend kwadraat en een onjuiste logaritmische factor voor het geval $d > \alpha$.

Methodologie en Correcties
De erratum vervangt de gebrekkige argumenten door rigoureuze correcties op de volgende sleutelgebieden:

1. Definities van Scrambling en Interlacing: De oorspronkelijke definitie van orde-$d$ scrambling maakte gebruik van een inverse interlacing-afbeelding $D_d^{-1}$, die niet goed gedefinieerd is omdat $D_d$ niet surjectief is. De correctie definieert de orde-$d$ gescrambelde punt $Y$ direct als $Y = D_d(\Pi(z))$, waarbij $\Pi$ Owen's geneste uniforme scrambling toegepast op $ds$ coördinaten is, en $D_d$ de digitale interlacing-afbeelding is. Propositie 2.1 stelt vast dat $Y$ uniform verdeeld is over $[0, 1)^s$ zonder dat een inverse afbeelding vereist is.

2. Covariantiestructuur: Een nieuwe covariantie-stelling (Propositie 2.2) is afgeleid voor de geneste uniforme scrambling. Deze propositie berekent de verwachting van producten van Walsh-functies onder scrambling met behulp van een blok-geïndexeerde functie $\chi_\ell$. Dit vervangt de vorige indexering in Lemma 6 en behandelt identieke punten zonder ad-hoc conventies.

3. Variantie-decompositie: De variantie van de estimator is opnieuw afgeleid met behulp van de gecorrigeerde covariantie. De oorspronkelijke Vergelijking (3.2) telde onjuist over alle frequenties; de gecorrigeerde versie (Vergelijking 10) telt over niet-nul frequenties. Lemma 3.1 biedt de gecorrigeerde variantie-decompositie:
   $$\text{Var}[\hat{I}_N(f)] = \sum_{\ell \in \mathbb{N}_0^{ds} \setminus \{0\}} \sigma_{d,\ell,s}^2(f) \Gamma_{d,\ell}(P)$$
   waarbij $\sigma_{d,\ell,s}^2(f)$ de blok-energie van de Walsh-coëfficiënten vertegenwoordigt en $\Gamma_{d,\ell}(P)$ de winst-coëfficiënt is die afhankelijk is van de net-structuur.

4. Intrekking van de Eindverschil-variatie: Het artikel trekt expliciet de claim in dat de eindverschil-variatie $V_\alpha$ samenvalt met de weergegeven Sobolev-norm. Een tegenvoorbeeld ($f(x) = x^2/2$) laat zien dat de variatie (een supremum over partities) niet gelijk is aan het integraal van de gekwadrateerde afgeleide. Bovendien bleek de toelaatbaarheid van incrementen in het oorspronkelijke bewijs in Appendix B gebrekkig te zijn. Bijgevolg worden alle resultaten die steunen op $V_\alpha$ ingetrokken.

5. Directe Sobolev-schatting: Om de ingetrokken eindverschil-argument te vervangen, wordt een directe Sobolev-schatting (Lemma 5.2) verstrekt. Dit lemma begrenst de blok-energie $\sigma_{d,\ell,s}(f)$ direct in termen van de ongeankerde gemengde Sobolev-norm $\|f\|_{H^{s,\alpha}}$. Het bewijs maakt gebruik van een eindige Walsh-blok-identiteit, waarbij paren cellen worden opgelost in fijne kinderen en de integraal-Minkowski-ongelijkheid wordt toegepast om de eindverschillen te begrenzen door de gemengde partiële afgeleiden.

6. Gecorrigeerde Logaritmische Factor: Het bewijs voor het geval $d > \alpha$ is herzien om correct rekening te houden met de verval van de Walsh-blokken. De oorspronkelijke stelling bevatte een onjuiste logaritmische macht. De gecorrigeerde Stelling 6.2 stelt de variantie-grens vast:
   $$\text{Var}[\hat{I}_N(f)] \le C_{b,s,\alpha,d} \|f\|_{H^{s,\alpha}}^2 b^{-(2\rho+1)M} (M+2)^{ds-1}$$
   waarbij $\rho = \min(\alpha, d)$ en $M = m-t$.

Belangrijkste Resultaten
Het primaire resultaat van de erratum is de validatie van de optimale convergentiesnelheid voor het gecorrigeerde bewijs.

• RMSE-snelheid: Voor een familie van nets met een vaste grondtal $b$, dimensie $s$, orde $d$ en gladheid $\alpha$, voldoet de wortelgetal van de kwadratische gemiddelde fout aan:
  $$\text{RMSE} = O\left(N^{-\rho - 1/2} (\log N)^{(ds-1)/2}\right) = O\left(N^{-\rho - 1/2 + \varepsilon}\right)$$
  voor elke $\varepsilon > 0$, waarbij $\rho = \min(\alpha, d)$.
• Optimaliteit: De algebraïsche exponent $-\rho - 1/2$ blijft optimaal voor de gemengde Sobolev-klasse.
• Specifieke Correcties:
  • De variantie-grens in Stelling 10 is gecorrigeerd door een kwadraat op de Sobolev-norm toe te voegen.
  • De logaritmische factor is gecorrigeerd van de oorspronkelijke stelling naar $(ds-1)/2$.
  • Typografische en indexeringsfouten met betrekking tot de digitale interlacing-afbeelding, het domein van $E_d$, en de definitie van eindverschillen zijn rechtgezet in Sectie 7.

Betekenis en Claims
De auteur stelt dat het belangrijkste resultaat betreffende gladde functies van het oorspronkelijke artikel geldig blijft: orde-$d$ geneste uniform gescrambelde digitale nets bereiken de optimale convergentiesnelheid voor gladde integranden. De numerieke experimenten en de optimale algebraïsche exponent die in het oorspronkelijke werk worden aangehaald, blijven ongewijzigd. De specifieke bewijsplaats die gebruikmaakt van de eindverschil-variatie wordt echter volledig verworpen als ongeldig. De betekenis van deze erratum ligt in het bieden van een wiskundig sluitend, direct Sobolev-gebaseerd bewijs dat de theoretische fundamenten van de convergentiesnelheden herstelt, de variantie-grenzen en logaritmische factoren corrigeert zonder de fundamentele conclusie over de efficiëntie van hogere-orde gescrambelde digitale nets te wijzigen.