2025-11-25T13:52:17.039270

On weak convergence of Gaussian conditional distributions

Lumpp, Drton
Weak convergence of joint distributions generally does not imply convergence of conditional distributions. In particular, conditional distributions need not converge when joint Gaussian distributions converge to a singular Gaussian limit. Algebraically, this is due to the fact that at singular covariance matrices, Schur complements are not continuous functions of the matrix entries. Our results lay out special conditions under which convergence of Gaussian conditional distributions nevertheless occurs, and we exemplify how this allows one to reason about conditional independence in a new class of graphical models.
academic

Über schwache Konvergenz von bedingten Gaußschen Verteilungen

Grundinformationen

  • Paper-ID: 2510.12412
  • Titel: On weak convergence of Gaussian conditional distributions
  • Autoren: Sarah Lumpp, Mathias Drton
  • Klassifikation: math.ST stat.TH
  • Veröffentlichungsdatum: 14. Oktober 2025 (arXiv-Preprint)
  • Paper-Link: https://arxiv.org/abs/2510.12412

Zusammenfassung

Die schwache Konvergenz von Gelenkverteilungen impliziert in der Regel nicht die Konvergenz von bedingten Verteilungen. Insbesondere wenn gemeinsame Gaußsche Verteilungen gegen eine singuläre Gaußsche Grenzverteilung konvergieren, konvergieren die bedingten Verteilungen nicht notwendigerweise. Aus algebraischer Perspektive liegt dies daran, dass das Schur-Komplement bei singulären Kovarianzmatrizen keine stetige Funktion der Matrixelemente ist. Die Ergebnisse dieser Arbeit erläutern spezielle Bedingungen, unter denen Gaußsche bedingte Verteilungen dennoch konvergieren, und zeigen durch Beispiele, wie dies es ermöglicht, in einer neuen Klasse von Graphenmodellen über bedingte Unabhängigkeit zu argumentieren.

Forschungshintergrund und Motivation

Kernproblem

Das Kernproblem dieser Untersuchung ist: Wenn eine Folge von Gaußschen Verteilungen X(m)N(0,Σ(m))X^{(m)} \sim N(0, Σ^{(m)}) schwach gegen eine singuläre Gaußsche Verteilung XN(0,Σ)X \sim N(0, Σ) konvergiert, konvergieren dann die bedingten Verteilungen N(0,ΣRS(m))N(0, Σ_{R|S}^{(m)}) gegen N(0,ΣRS)N(0, Σ_{R|S})?

Bedeutung des Problems

  1. Theoretische Relevanz: Die Konvergenz bedingter Verteilungen ist ein grundlegendes Problem der Wahrscheinlichkeitstheorie und von großer Bedeutung für das Verständnis des Grenzverhaltens von Verteilungen
  2. Anwendungswert: In Graphenmodellen hängt die Analyse von bedingter Unabhängigkeit von den Eigenschaften bedingter Verteilungen ab, besonders in neuartigen Graphenmodellen, die durch kontinuierliche Lyapunov-Gleichungen definiert werden

Einschränkungen bestehender Methoden

  1. Unstetigkeit im allgemeinen Fall: Die Pseudoinverse ()+(·)^+ ist bei singulären Matrizen unstetig, was dazu führt, dass die bedingte Kovarianzmatrix ΣRS=ΣR,RΣR,S(ΣS,S)+ΣS,RΣ_{R|S} = Σ_{R,R} - Σ_{R,S}(Σ_{S,S})^+Σ_{S,R} möglicherweise nicht konvergiert
  2. Einschränkungen in Spezialfällen: Bekannte Ergebnisse gelten hauptsächlich für Matrixfolgen konstanten Ranges, während in praktischen Anwendungen häufig Matrizen vollen Ranges gegen singuläre Matrizen konvergieren

Forschungsmotivation

Inspiriert durch die Arbeiten von Boege et al. (2025) zur bedingten Unabhängigkeit betrachten die Autoren den anspruchsvolleren Fall: Jede Σ(m)Σ^{(m)} ist invertierbar, während die Grenz-Kovarianzmatrix ΣΣ singulär ist.

Kernbeiträge

  1. Konvergenzbedingungen etablieren: Ausreichende Bedingungen für die Konvergenz Gaußscher bedingter Verteilungen werden vorgestellt, auch im singulären Grenzfall
  2. Verallgemeinerung des Matrixdeterminanten-Lemmas: Eine erweiterte Version des Matrixdeterminanten-Lemmas wird entwickelt, um das asymptotische Verhalten singulärer Matrizen zu analysieren
  3. Anwendung auf Graphenmodelle: Die Anwendung der theoretischen Ergebnisse auf die Analyse bedingter Unabhängigkeit in Graphenmodellen, die durch kontinuierliche Lyapunov-Gleichungen definiert werden, wird demonstriert
  4. Spezialfall Toeplitz-Matrizen: Ein intuitiver alternativer Beweis wird für spezielle Toeplitz-Matrizen bereitgestellt

Methodische Details

Aufgabendefinition

Gegeben ist eine Folge von Gaußschen Zufallsvektoren X(m)N(0,Σ(m))X^{(m)} \sim N(0, Σ^{(m)}), die schwach gegen XN(0,Σ)X \sim N(0, Σ) konvergiert, wobei Σ(m)Σ^{(m)} invertierbar, aber ΣΣ singulär ist. Für disjunkte Teilmengen S,R[p]S, R ⊆ [p] ist zu bestimmen, ob die bedingte Kovarianzmatrix ΣRS(m)Σ_{R|S}^{(m)} gegen ΣRSΣ_{R|S} konvergiert.

Zentrale theoretische Ergebnisse

Hauptsatz (Satz 1.2)

Sei Σ(m)Σ^{(m)} eine Folge invertierbarer symmetrischer p×pp×p-Matrizen mit singulärer Grenzmatrix ΣΣ und Entwicklung: Σ(m)=Σ+1mΣ(1,)+o(1m)Σ^{(m)} = Σ + \frac{1}{m}Σ^{(1,∞)} + o\left(\frac{1}{m}\right)

Für S[p]S ⊂ [p] mit k=S<pk = |S| < p, R=[p]SR = [p] \setminus S, konvergiert die bedingte Kovarianzmatrix ΣRS(m)Σ_{R|S}^{(m)} gegen ΣRSΣ_{R|S}, wenn die Matrizen ΣΣ und Σ(1,)Σ^{(1,∞)} folgende Bedingungen erfüllen:

(i) f(U,ΣS,S(1,))0f(U, Σ_{S,S}^{(1,∞)}) ≠ 0, wobei UU aus der symmetrischen Rangzerlegung ΣS,S=UBUTΣ_{S,S} = UBU^T stammt
(ii) im(ΣS,R)im(ΣS,S)\text{im}(Σ_{S,R}) ⊆ \text{im}(Σ_{S,S})

Definition der Schlüsselfunktion

Die Funktion f:Rk×r×Rk×kRf: \mathbb{R}^{k×r} × \mathbb{R}^{k×k} → \mathbb{R} ist definiert als: f(U,G)=(1)s([kr])α[k]α=kr(1)s(α)det(G[k],αU)det(U[k]α,[r])f(U,G) = (-1)^{s([k-r])} \sum_{\substack{α⊆[k] \\ |α|=k-r}} (-1)^{s(α)} \det\begin{pmatrix} G_{[k],α} & U \end{pmatrix} \det(U_{[k]\setminus α,[r]})

Technische Innovationen

1. Verallgemeinerung des Matrixdeterminanten-Lemmas

Lemma 3.1: Für Matrizen U,VRk×rU, V ∈ \mathbb{R}^{k×r}, GRk×kG ∈ \mathbb{R}^{k×k}, BRr×rB ∈ \mathbb{R}^{r×r} mit r<kr < k: det(UBVT+1mG)=1mkrdet(B)fasy(U,V,G)+o(1mkr)\det\left(UBV^T + \frac{1}{m}G\right) = \frac{1}{m^{k-r}} \cdot \det(B) \cdot f_{\text{asy}}(U, V, G) + o\left(\frac{1}{m^{k-r}}\right)

2. Anwendung der Cramerschen Regel

Jedes Element der bedingten Kovarianzmatrix wird als: vmTAm1wm=det(Am[j])det(Am)v_m^T A_m^{-1} w_m = \frac{\det(A_m[j])}{\det(A_m)} ausgedrückt und das asymptotische Verhalten von Zähler und Nenner wird analysiert, um die Konvergenz zu beweisen.

3. Rauschbehandlung

Lemma 3.2 behandelt die Auswirkungen zusätzlicher Rauschterme D(m)=o(1/m)D^{(m)} = o(1/m) und stellt sicher, dass der Hauptasymptotikterm nicht beeinträchtigt wird.

Experimentelle Einrichtung

Anwendungsszenario: Kontinuierliche Lyapunov-Gleichung

Betrachtet werden Gaußsche Verteilungen, die die kontinuierliche Lyapunov-Gleichung erfüllen: MΣ+ΣMT+2Ip=0MΣ + ΣM^T + 2I_p = 0 wobei MM eine stabile Matrix ist, deren Nullmuster die Graphenstruktur kodiert.

Konkrete Beispiele

Beispiel 2.1: Trek T4:1234T_4: 1 ← 2 → 3 → 4

Konstruktion einer stabilen Matrix:

-1 & 1 & 0 & 0 \\ 0 & -1 & 0 & 0 \\ 0 & 1 & -1 & 0 \\ 0 & 0 & 1 & -1 \end{pmatrix}$$ Erhaltene Lösung: $$Σ_4 = \begin{pmatrix} 3/2 & 1/2 & 1/2 & 3/8 \\ 1/2 & 1 & 1/2 & 1/4 \\ 1/2 & 1/2 & 3/2 & 7/8 \\ 3/8 & 1/4 & 7/8 & 15/8 \end{pmatrix}$$ Verifikation: $(Σ_4)_{14|23} = 1/10 ≠ 0$. ## Experimentelle Ergebnisse ### Anwendung auf Toeplitz-Matrizen Für $n$-dimensionale symmetrische Toeplitz-Matrizen: $$T_n = (-|i-j|)_{i,j=1,\ldots,n}$$ **Lemma 4.1** beweist für spezifische Matrizen $U$ und $G$: $$f(U,G) = 2^{k+1-|\{1,n_l\}|-|\{1,n_r\}|} ≠ 0$$ ### Verifikation der Hauptergebnisse **Proposition 4.2** (Proposition 4.5 von Boege et al. 2025): Für jeden Trek beliebiger Länge $T: 1 ← ··· → p$ existiert eine Kovarianzmatrix $Σ$, die die kontinuierliche Lyapunov-Gleichung erfüllt, so dass $X_1 \not\perp X_p | X_2, \ldots, X_{p-1}$. Der Beweis wird durch die Konstruktion einer Matrixfolge $M^{(m)}$ erbracht, so dass die entsprechenden Lösungen die Bedingungen von Satz 1.2 erfüllen und somit die bedingte Unabhängigkeit nicht besteht. ## Verwandte Arbeiten ### Theoretische Grundlagen 1. **Stewart (1969), Ben-Israel and Greville (2003)**: Stetigkeit der Pseudoinversen für Matrixfolgen konstanten Ranges 2. **Lauritzen (1996, 2024)**: Allgemeine Theorie der Konvergenz bedingter Verteilungen 3. **Bérczi et al. (2023)**: Diese Arbeit verallgemeinert deren Ergebnisse aus Abschnitt 7 ### Anwendungshintergrund 1. **Dettling et al. (2023), Varando and Hansen (2020)**: Kontinuierliche Lyapunov-Graphenmodelle 2. **Boege et al. (2025)**: Charakterisierung bedingter Unabhängigkeit in Graphenmodellen ## Schlussfolgerungen und Diskussion ### Hauptschlussfolgerungen 1. Unter bestimmten Bedingungen konvergieren Gaußsche bedingte Verteilungen auch dann, wenn die Grenz-Kovarianzmatrix singulär ist 2. Bedingung (i) stellt sicher, dass der Nenner des Schur-Komplements nicht gegen Null geht, Bedingung (ii) gewährleistet die Kompatibilität der Bildräume 3. Die theoretischen Ergebnisse werden erfolgreich auf die Analyse bedingter Unabhängigkeit in kontinuierlichen Lyapunov-Graphenmodellen angewendet ### Einschränkungen 1. **Bedingungsrestriktivität**: Es müssen spezifische Matrixstruktur-Bedingungen erfüllt sein 2. **Rechenkomplexität**: Die Berechnung der Funktion $f$ beinhaltet Mehrfachsummen und kann für große Matrizen rechnerisch schwierig sein 3. **Anwendungsbereich**: Hauptsächlich auf Gaußsche Verteilungen ausgerichtet; die Erweiterbarkeit auf andere Verteilungsfamilien ist ungeklärt ### Zukünftige Richtungen 1. **Bedingungslokkerung**: Suche nach allgemeineren Konvergenzbedingungen 2. **Rechnerische Optimierung**: Entwicklung effizienterer Algorithmen zur Bedingungsverifikation 3. **Nicht-Gaußsche Erweiterung**: Verallgemeinerung der Ergebnisse auf andere Verteilungsfamilien ## Tiefgreifende Bewertung ### Stärken 1. **Theoretische Strenge**: Bietet einen vollständigen mathematischen Beweisrahmen, von der Verallgemeinerung des Matrixdeterminanten-Lemmas bis zum Hauptsatz 2. **Praktischer Wert**: Löst praktische Probleme in Graphenmodellen, insbesondere die Anwendung auf kontinuierliche Lyapunov-Gleichungen 3. **Methodische Innovation**: Geschickte Kombination von Matrixanalyse, Wahrscheinlichkeitstheorie und Graphentheorie 4. **Klare Darstellung**: Die Papierstruktur ist klar, mit reichhaltigen Beispielen, die das Verständnis erleichtern ### Mängel 1. **Bedingungskomplexität**: Die Verifikation der Satzvoraussetzungen kann schwierig sein, besonders die Berechnung der Funktion $f$ 2. **Anwendungsbegrenzungen**: Hauptsächlich auf spezifische Matrixstrukturen anwendbar; die Allgemeinheit ist verbesserungsbedürftig 3. **Numerische Stabilität**: Bei Matrizen nahe der Singularität können numerische Stabilitätsprobleme auftreten ### Einflussfähigkeit 1. **Theoretischer Beitrag**: Bietet neue theoretische Werkzeuge für die Konvergenz bedingter Verteilungen unter singulären Grenzen 2. **Anwendungsperspektiven**: Hat potenzielle Anwendungen in Graphenmodellen, Zeitreihenanalyse und anderen Bereichen 3. **Methodische Inspiration**: Bietet neue Ansätze für die Behandlung ähnlicher Matrixgrenzwertprobleme ### Anwendungsszenarien 1. **Graphenmodellanalyse**: Besonders für Modelle, die durch kontinuierliche Lyapunov-Gleichungen definiert werden 2. **Zeitreihen**: Multivariate Zeitreihen mit spezifischen Kovarianzstrukturen 3. **Statistische Inferenz**: Statistische Probleme, die die Analyse des Grenzverhaltens bedingter Verteilungen erfordern ## Literaturverzeichnis Das Papier zitiert mehrere wichtige verwandte Arbeiten, darunter: - Ben-Israel and Greville (2003): Theorie der verallgemeinerten Inversen - Lauritzen (1996, 2024): Graphenmodelle und Theorie bedingter Verteilungen - Boege et al. (2025): Anwendung bedingter Unabhängigkeit in Diffusionsprozessen - Stewart (1969): Stetigkeitstheorie der verallgemeinerten Inversen --- Dieses Papier leistet wichtige Beiträge auf theoretischer und anwendungsorientierter Ebene und bietet neue mathematische Werkzeuge zum Verständnis des Verhaltens bedingter Verteilungen unter singulären Grenzen mit erfolgreicher Anwendung auf die Analyse bedingter Unabhängigkeit in Graphenmodellen.