2025-11-10T03:05:00.236892

Stein's method for Fréchet approximation: a regularly varying functions approach

Paul, Guillaume, Yvik
We develop a variant of Stein's method of comparison of generators to bound the Kolmogorov, total variation, and Wasserstein-1 distances between distributions on the real line. Our discrepancy is expressed in terms of the ratio of reverse hazard rates; it therefore remains tractable even when density derivatives are intractable. Our main application concerns the approximation of normalized extremes by Fréchet laws. In this setting, the new discrepancy provides a quantitative measure of distributional proximity in terms of the average regular variation at infinity of the underlying cumulative distribution function. We illustrate the approach through explicit computations for maxima of Pareto, Cauchy, and Burr~XII distributions.
academic

Steins Methode für Fréchet-Approximation: ein Ansatz mit regelmäßig variierenden Funktionen

Grundinformationen

  • Papier-ID: 2510.14016
  • Titel: Stein's method for Fréchet approximation: a regularly varying functions approach
  • Autoren: Paul Mansanarez, Guillaume Poly, Yvik Swan
  • Klassifizierung: math.PR (Wahrscheinlichkeitstheorie)
  • Veröffentlichungsdatum: 17. Oktober 2025
  • Papierlink: https://arxiv.org/abs/2510.14016

Zusammenfassung

In diesem Artikel wird eine Variante der Steinschen Methode entwickelt, um Erzeugeroperatoren zu vergleichen und die Kolmogorov-Distanz, die Totalvariationsdistanz und die Wasserstein-1-Distanz zwischen Verteilungen auf der reellen Linie zu begrenzen. Die Diskrepanz der Methode wird durch das Verhältnis der umgekehrten Hazard-Raten ausgedrückt und bleibt daher auch in Fällen handhabbar, in denen Dichtederivate schwer zu handhaben sind. Die Hauptanwendung konzentriert sich auf die Approximation normalisierter Extremwerte durch das Fréchet-Gesetz. In dieser Einstellung bietet die neue Diskrepanz ein quantitatives Maß für die Verteilungsnähe, das auf der durchschnittlichen regelmäßigen Variation der zugrunde liegenden kumulativen Verteilungsfunktion im Unendlichen basiert. Die Methode wird durch explizite Berechnungen für die Maxima von Pareto-, Cauchy- und Burr-XII-Verteilungen illustriert.

Forschungshintergrund und Motivation

Problemhintergrund

  1. Quantifizierungsproblem der Approximation in der Extremwerttheorie: In der Extremwerttheorie ist die Fréchet-Verteilung eine von drei verallgemeinerten Extremwertverteilungen (GEV) und ist die einzige mögliche Grenzverteilung für das Maximum unabhängig identisch verteilter Zufallsvariablen nach angemessener Normalisierung. Die Quantifizierung der Konvergenzgeschwindigkeit und -genauigkeit ist ein wichtiges theoretisches und praktisches Problem.
  2. Einschränkungen bestehender Methoden:
    • Traditionelle Steinsche Methoden beruhen auf dem Vergleich von Score-Funktionen (Ableitungen des Logarithmus der Wahrscheinlichkeitsdichte), die in der GEV-Einstellung schwer zu handhaben sind
    • Bestehende Konvergenzratenergebnisse sind technisch anspruchsvoll und schwer anzuwenden
    • Es fehlen quantifizierbare Werkzeuge, um die Verteilungsnähe mit der Theorie der regelmäßigen Variation zu verbinden
  3. Forschungsmotivation:
    • Entwicklung einer neuen Steinschen Methodenvariante, die die Berechnung von Dichtederivaten vermeidet
    • Bereitstellung handhabbarerer Grenzen in der GEV-Approximation
    • Etablierung einer quantitativen Verbindung zwischen Verteilungsnähe und regelmäßiger Variation

Kernbeiträge

  1. Neue Steinsche Diskrepanzmaße: Vorschlag einer Steinschen Diskrepanz basierend auf dem Verhältnis der umgekehrten Hazard-Raten, die die Rechenkomplexität von Dichtederivaten vermeidet
  2. Einheitliche Grenzen für drei Distanzen: Bereitstellung einheitlicher Obergrenzen für die Kolmogorov-Distanz, Totalvariationsdistanz und Wasserstein-1-Distanz
  3. Quantitative Charakterisierung der regelmäßigen Variation: Verbindung der Steinschen Diskrepanz mit der regelmäßigen Variation der kumulativen Verteilungsfunktion im Unendlichen und Bereitstellung einer quantitativen Konvergenzcharakterisierung
  4. Konkrete Anwendungsbeispiele: Demonstration der praktischen Anwendbarkeit der Methode durch explizite Berechnungen für Pareto-, Cauchy- und Burr-XII-Verteilungen

Methodische Erläuterung

Aufgabendefinition

Gegeben seien zwei kumulative Verteilungsfunktionen P und Q, die Annahme 0 erfüllen. Das Ziel besteht darin, die Wahrscheinlichkeitsdistanzen zwischen ihnen zu begrenzen:

  • Kolmogorov-Distanz: Kol(P,Q)\text{Kol}(P,Q)
  • Totalvariationsdistanz: TV(P,Q)\text{TV}(P,Q)
  • Wasserstein-1-Distanz: Wass(P,Q)\text{Wass}(P,Q)

Kernannnahmen

Annahme 0: Es existieren eine Anfangsmasse f00f_0 \geq 0, ein linker Endpunkt cF[,+)c_F \in [-\infty,+\infty) und eine Dichtefunktion f:RR+f: \mathbb{R} \to \mathbb{R}_+, so dass die kumulative Verteilungsfunktion gegeben ist durch: F(x)={0,x<cFf0+cFxf(u)du,xcFF(x) = \begin{cases} 0, & x < c_F \\ f_0 + \int_{c_F}^x f(u) du, & x \geq c_F \end{cases}

Definition der Steinschen Diskrepanz

Für Verteilungen P und Q, die Annahme 0 erfüllen, wird die umgekehrte Hazard-Rate definiert als: rF(x)=f(x)F(x)I(cF,+)(x)=(logF)(x)I(cF,+)(x)r_F(x) = \frac{f(x)}{F(x)}I_{(c_F,+\infty)}(x) = (\log F)'(x)I_{(c_F,+\infty)}(x)

Hauptsteinsche Diskrepanz: Δ(QP):=cQ+1rP(x)rQ(x)q(x)dx\Delta(Q | P) := \int_{c_Q}^{+\infty} \left|1 - \frac{r_P(x)}{r_Q(x)}\right| q(x) dx

Gewichtete Steinsche Diskrepanz (für Wasserstein-Distanz): Δw(QP):=cQ+x1rP(x)rQ(x)q(x)dx\Delta_w(Q | P) := \int_{c_Q}^{+\infty} x \left|1 - \frac{r_P(x)}{r_Q(x)}\right| q(x) dx

Haupttheoretische Ergebnisse

Satz 1.1 (Grenzen für Kolmogorov- und Totalvariationsdistanz): Kol(P,Q)Δ(QP)+q0\text{Kol}(P,Q) \leq \Delta(Q | P) + q_0TV(P,Q)2Δ(QP)+q0\text{TV}(P,Q) \leq 2\Delta(Q | P) + q_0

Satz 1.2 (Grenze für Wasserstein-1-Distanz): Wass(P,Q)2μΔ(QP)+3Δw(QP)\text{Wass}(P,Q) \leq 2\mu\Delta(Q | P) + 3\Delta_w(Q | P)

Drei-Schritte-Rahmen der Steinschen Methode

  1. Konstruktion des Stein-Operators: Definition des Reparametrisierungsoperators Apϕ:=(Pˉpϕ+ϕ)I(cP,+)A^p\phi := \left(\frac{\bar{P}}{p}\phi' + \phi\right)I_{(c_P,+\infty)}
  2. Lösen der Stein-Gleichung: Für eine gegebene Testfunktion hh wird die ODE gelöst Pˉpϕ+ϕ=hP(h)\frac{\bar{P}}{p}\phi' + \phi = h - P(h)
  3. Transferprinzip: Nutzung der Operatordifferenzen zwischen verschiedenen Verteilungen zur Etablierung von Distanzgrenzen

Experimentelle Einrichtung

Theoretische Verifikation

Das Papier verifiziert die Wirksamkeit der Methode durch die folgenden Verteilungen:

  1. Pareto-Verteilung: F(x)=(1xα)I(1,)(x)F(x) = (1-x^{-\alpha})I_{(1,\infty)}(x)
  2. Cauchy-Verteilung: Standardmäßige Cauchy-Verteilung
  3. Burr-XII-Verteilung: F(x)=1(1+xα)τF(x) = 1-(1+x^\alpha)^{-\tau}, x>0x > 0
  4. Log-Logistische Verteilung: F(x)=xα1+xαI(0,)(x)F(x) = \frac{x^\alpha}{1+x^\alpha}I_{(0,\infty)}(x)

Bewertungsmetriken

  • Genaue Konstanten der Konvergenzraten
  • Numerische Berechnungsverifikation
  • Vergleich mit bekannten theoretischen Ergebnissen

Experimentelle Ergebnisse

Ergebnisse für die Pareto-Verteilung

Für die Pareto-Verteilung wird eine genaue Konvergenzrate erhalten: Δ(FnΦα)=1n+1\Delta(F_n | \Phi_\alpha) = \frac{1}{n+1} Dies stimmt vollständig mit den bekannten Konvergenzraten der Kolmogorov- und Totalvariationsdistanzen überein.

Ergebnisse für die Cauchy-Verteilung

Für die Cauchy-Verteilung wird durch Integralsubstitution erhalten: Δ(Φ1Fn)=E[1R(En)]\Delta(\Phi_1 | F_n) = \mathbb{E}[|1-R(E_n)|] wobei EnExp(n)E_n \sim \text{Exp}(n), und es wird eine Konvergenzrate von O(1/n)O(1/n) nachgewiesen.

Ergebnisse für die Burr-XII-Verteilung

Numerische Berechnungen zeigen eine Konvergenzrate von O(n1/τ)O(n^{-1/\tau}), die mit theoretischen Erwartungen übereinstimmt.

Vergleich zwischen Fréchet-Verteilungen

Wenn P und Q beide Fréchet-Verteilungen sind, wird ein expliziter Ausdruck für die Diskrepanz gegeben, der Gamma-Funktionen und unvollständige Gamma-Funktionen beinhaltet.

Verwandte Arbeiten

Steinsche Methode in der Extremwerttheorie

  • Feidt (2013): Steinsche Methode für multivariate Extremwerte
  • Bartholmé & Swan (2013): Vergleichsstrategie für Erzeugeroperatoren
  • Costacèque-Cecchi (2024): Rahmen basierend auf Markov-Halbgruppen
  • Kusumoto & Takeuchi (2020): Konvergenzraten durch Stein-Gleichungen

Theorie der regelmäßigen Variation

  • Karamata-Theorie: Grundlegende Theorie regelmäßig variierender Funktionen
  • Anziehung in Extremwertbereichen: Charakterisierung von Konvergenzbedingungen durch regelmäßige Variation

Schlussfolgerungen und Diskussion

Hauptschlussfolgerungen

  1. Methodische Wirksamkeit: Die neue Steinsche Diskrepanz bietet wirksame Grenzen für Verteilungsdistanzen
  2. Rechnerische Vorteile: Die Vermeidung von Dichtederivaten macht die Methode bei komplexen Verteilungen leichter anwendbar
  3. Theoretische Verbindung: Etablierung einer quantitativen Verbindung zwischen Verteilungsnähe und regelmäßiger Variation

Einschränkungen

  1. Annahmebedingungen: Erfordernis der Erfüllung von Annahme 0 begrenzt den Anwendungsbereich
  2. Trägerbedingungen: Anforderungen an die relative Position der Verteilungsträger
  3. Technische Bedingungen: Einige Ergebnisse erfordern zusätzliche Beschränktheitsbedingungen

Zukünftige Richtungen

  1. Verallgemeinerung auf multivariate Fälle: Erweiterung auf multivariate Extremwertverteilungen
  2. Andere GEV-Verteilungen: Anwendung auf Gumbel- und Weibull-Verteilungen
  3. Analyse höherer Ordnungsterme: Untersuchung feinerer asymptotischer Entwicklungen

Tiefgreifende Bewertung

Stärken

  1. Theoretische Innovation: Vorschlag einer neuen auf umgekehrten Hazard-Raten basierenden Steinschen Diskrepanz, die technische Schwierigkeiten traditioneller Methoden vermeidet
  2. Einheitlicher Rahmen: Einheitliche Behandlung von drei wichtigen Wahrscheinlichkeitsdistanzen
  3. Praktische Durchführbarkeit: Demonstration der rechnerischen Machbarkeit der Methode durch konkrete Beispiele
  4. Theoretische Tiefe: Tiefe Kombination der Steinschen Methode mit der Theorie der regelmäßigen Variation

Mängel

  1. Annahmebeschränkungen: Annahme 0 stellt starke Anforderungen an die Verteilungsform und kann einige praktische Anwendungen ausschließen
  2. Straffheit der Grenzen: Unzureichende Diskussion der Straffheit der erhaltenen Grenzen
  3. Numerische Experimente: Mangel an umfangreichen numerischen Experimenten zur Verifikation der praktischen Leistung der Methode

Einflussfähigkeit

  1. Theoretischer Beitrag: Bereitstellung neuer Werkzeuge für die quantifizierte Approximation in der Extremwerttheorie
  2. Methodologischer Wert: Die neue Variante der Steinschen Methode kann andere Anwendungen inspirieren
  3. Praktischer Wert: Potenzielle Anwendungen in Risikomanagement, Versicherung und anderen Bereichen

Anwendungsszenarien

  • Verteilungsapproximation in der Extremwertstatistik
  • Tail-Risk-Bewertung im Risikomanagement
  • Modellierung extremer Ereignisse in der Versicherungsmathematik
  • VaR- und ES-Schätzung im Finanzwesen

Literaturverzeichnis

Das Papier zitiert klassische Literatur aus der Extremwerttheorie, der Theorie der regelmäßigen Variation und der Steinschen Methode, einschließlich:

  • Bingham, Goldie & Teugels (1987): Regular Variation
  • de Haan & Ferreira (2006): Extreme Value Theory
  • Resnick (2008): Extreme Values, Regular Variation, and Point Processes
  • Ley, Reinert & Swan (2017): Stein's Method for Comparison of Univariate Distributions

Dieses Papier stellt theoretisch eine wichtige Innovation dar und wendet die Steinsche Methode erfolgreich auf das Problem der Extremwertverteilungsapproximation an, wobei es neue Analysewerkzeuge für dieses Gebiet bereitstellt. Die praktische Anwendbarkeit der Methode wird durch Berechnungen mehrerer klassischer Verteilungen verifiziert und zeigt gute Anwendungsaussichten.