Sobre la convergencia débil de distribuciones condicionales gaussianas

Información Básica

• ID del Artículo: 2510.12412
• Título: On weak convergence of Gaussian conditional distributions
• Autores: Sarah Lumpp, Mathias Drton
• Clasificación: math.ST stat.TH
• Fecha de Publicación: 14 de octubre de 2025 (preimpresión arXiv)
• Enlace del Artículo: https://arxiv.org/abs/2510.12412

Resumen

La convergencia débil de distribuciones conjuntas generalmente no implica la convergencia de distribuciones condicionales. En particular, cuando distribuciones gaussianas conjuntas convergen a un límite gaussiano singular, las distribuciones condicionales no necesariamente convergen. Desde una perspectiva algebraica, esto ocurre porque el complemento de Schur no es una función continua de los elementos de la matriz en matrices de covarianza singulares. Los resultados de este artículo elucidan las condiciones especiales bajo las cuales las distribuciones condicionales gaussianas aún convergen, e ilustran cómo esto permite razonar sobre independencia condicional en una nueva clase de modelos gráficos.

Contexto de Investigación y Motivación

Problema Central

El problema central que aborda esta investigación es: cuando una secuencia de distribuciones gaussianas $X^{(m)} \sim N(0, Σ^{(m)})$ converge débilmente a una distribución gaussiana singular $X \sim N(0, Σ)$, ¿convergen las distribuciones condicionales $N(0, Σ_{R|S}^{(m)})$ a $N(0, Σ_{R|S})$?

Importancia del Problema

1. Significado Teórico: La convergencia de distribuciones condicionales es un problema fundamental en teoría de probabilidades, con importancia crucial para comprender el comportamiento límite de distribuciones
2. Valor Aplicado: En modelos gráficos, el análisis de independencia condicional depende de las propiedades de distribuciones condicionales, especialmente en nuevos modelos gráficos definidos por ecuaciones de Lyapunov continuas

Limitaciones de Métodos Existentes

1. Discontinuidad en Casos Generales: La operación de pseudoinversa $(·)^+$ es discontinua en matrices singulares, lo que causa que la matriz de covarianza condicional $Σ_{R|S} = Σ_{R,R} - Σ_{R,S}(Σ_{S,S})^+Σ_{S,R}$ pueda no converger
2. Restricciones en Casos Especiales: Los resultados conocidos se aplican principalmente a secuencias de matrices de rango constante, mientras que en aplicaciones prácticas frecuentemente se encuentran matrices de rango completo convergiendo a matrices singulares

Motivación de la Investigación

Inspirados por el trabajo de Boege et al. (2025) sobre independencia condicional, los autores consideran el caso más desafiante: cada $Σ^{(m)}$ es de rango completo, mientras que la matriz de covarianza límite $Σ$ es singular.

Contribuciones Principales

1. Establecimiento de Condiciones de Convergencia: Se proponen condiciones suficientes para la convergencia de distribuciones condicionales gaussianas, incluso en casos de límites singulares
2. Generalización del Lema del Determinante de Matrices: Se desarrolla una versión generalizada del lema del determinante de matrices para analizar el comportamiento asintótico de matrices singulares
3. Aplicación a Modelos Gráficos: Se demuestra cómo aplicar resultados teóricos al análisis de independencia condicional en modelos gráficos definidos por ecuaciones de Lyapunov continuas
4. Caso Especial de Matrices Toeplitz: Se proporciona una prueba alternativa intuitiva para matrices Toeplitz especiales

Explicación Detallada de Métodos

Definición de la Tarea

Dada una secuencia de vectores aleatorios gaussianos $X^{(m)} \sim N(0, Σ^{(m)})$ que converge débilmente a $X \sim N(0, Σ)$, donde $Σ^{(m)}$ es invertible pero $Σ$ es singular. Para subconjuntos disjuntos $S, R ⊆ [p]$, determinar si la matriz de covarianza condicional $Σ_{R|S}^{(m)}$ converge a $Σ_{R|S}$.

Resultados Teóricos Principales

Teorema Principal (Teorema 1.2)

Sea $Σ^{(m)}$ una secuencia de matrices simétricas $p×p$ invertibles con límite singular $Σ$, satisfaciendo la expansión: $Σ^{(m)} = Σ + \frac{1}{m}Σ^{(1,∞)} + o\left(\frac{1}{m}\right)$

Para $S ⊂ [p]$ con $k = |S| < p$, $R = [p] \setminus S$, si las matrices $Σ$ y $Σ^{(1,∞)}$ satisfacen:

(i) $f(U, Σ_{S,S}^{(1,∞)}) ≠ 0$, donde $U$ proviene de la descomposición de rango simétrico $Σ_{S,S} = UBU^T$
(ii) $\text{im}(Σ_{S,R}) ⊆ \text{im}(Σ_{S,S})$

entonces la matriz de covarianza condicional $Σ_{R|S}^{(m)}$ converge a $Σ_{R|S}$.

Definición de Función Clave

La función $f: \mathbb{R}^{k×r} × \mathbb{R}^{k×k} → \mathbb{R}$ se define como: $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]})$

Puntos de Innovación Técnica

1. Generalización del Lema del Determinante de Matrices

Lema 3.1: Para matrices $U, V ∈ \mathbb{R}^{k×r}$, $G ∈ \mathbb{R}^{k×k}$, $B ∈ \mathbb{R}^{r×r}$, donde $r < k$: $\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. Aplicación de la Regla de Cramer

Mediante la representación de cada elemento de la matriz de covarianza condicional como: $v_m^T A_m^{-1} w_m = \frac{\det(A_m[j])}{\det(A_m)}$ y analizando el comportamiento asintótico del numerador y denominador para probar la convergencia.

3. Manejo de Ruido

Lema 3.2 aborda el efecto de términos de ruido adicionales $D^{(m)} = o(1/m)$, asegurando que el término asintótico principal no sea perturbado.

Configuración Experimental

Escenario de Aplicación: Ecuación de Lyapunov Continua

Se considera distribuciones gaussianas que satisfacen la ecuación de Lyapunov continua: $MΣ + ΣM^T + 2I_p = 0$ donde $M$ es una matriz estable cuyo patrón de ceros codifica la estructura gráfica.

Ejemplos Concretos

Ejemplo 2.1: Trek $T_4: 1 ← 2 → 3 → 4$

Se construye la matriz estable: