Si una matriz de covarianza inversa es escasa, ¿qué puedo decir sobre la matriz de covarianza?

8

¿Cómo afecta la condición de dispersión en una matriz de covarianza inversa a la matriz de covarianza real?

Lepisma
fuente
Desafortunadamente, no es así
Yair Daon
Tiene que haber algo. Por ejemplo, la matriz de identidad es escasa, y también es inversa.
Aksakal

Respuestas:

4

Como ya comentó Yair, no existe una condición de escasez específica de la matriz de covarianza inversa que afecte a la matriz de covarianza real o viceversa. Cualquier otra cosa que no sean patrones triviales de matriz de dispersión (es decir, diagonal) no garantiza que se reflejarán tanto en una matriz particular como en su inverso. Incluso las matrices tridiagonales pueden tener inversas no dispersas fácilmente.

Para casos particulares donde la escasez de la matriz se produce en bloques, es posible que pueda obtener algunos resultados derivados del algoritmo pseudoinverso de la matriz de bloques que establece que:

[UNAsiCre]-1=[(UNA-sire-1C)-1-UNA-1si(re-CUNA-1si)-1-re-1C(UNA-sire-1C)-1(re-CUNA-1si)-1]

pero eso es probablemente sobre eso (puramente anecdótico, he tratado de imponer patrones de dispersión a través de la descomposición de Cholesky de una matriz PSD pero fallé en mi incursión de prueba y error). También es posible que desee considerar investigar el algoritmo de Cuthill-McKee (CM) si espera que alguna característica de adyacencia se refleje en la matriz de covarianza. El algoritmo CM permuta una matriz dispersa que tiene un patrón de dispersión simétrica en una forma de matriz de banda con un ancho de banda pequeño, esto podría ayudar a preservar cierta dispersión hacia las entradas fuera de la diagonal de la matriz inversa, pero eso no está garantizado. (La aplicación de CM, si es razonable, puede ser muy útil para aplicaciones particulares (por ejemplo, en rutinas de suavizado 2D) y puede acelerar significativamente sus cálculos.

usεr11852
fuente
2
(+1) Debido a que las matrices de diagonal de bloque de cualquier forma dada son un anillo, sus inversas (siempre que existan) tienen la misma estructura de diagonal de bloque y conservan gran parte del patrón de dispersión. Como un ejemplo extremo, las matrices diagonales son en diagonal de bloque, ejemplificando así el caso señalado por @Aksakal. Lo más lejos que se puede ir en esta dirección es conjugar matrices de bloques en diagonal por matrices de permutación (que obviamente conserva todas las entradas cero y no cero, pero solo las mueve).
whuber
1
(+1) Gracias por este comentario (respuesta corta realmente). Es realmente perspicaz. Definitivamente lo consideraré en el futuro.
usεr11852