Articulo de referencia

función de estructura de Kolmogorov

En 1973, Andrey Kolmogorov propuso un enfoque no probabilístico para la estadística y la selección de modelos . Sea cada dato una cadena binaria finita y un modelo un conjunto f...

En 1973, Andrey Kolmogorov propuso un enfoque no probabilístico para la estadística y la selección de modelos . Sea cada dato una cadena binaria finita y un modelo un conjunto finito de cadenas binarias. Consideremos clases de modelos que consisten en modelos con una complejidad de Kolmogorov máxima dada . La función de estructura de Kolmogorov de una cadena de datos individual expresa la relación entre la restricción del nivel de complejidad en una clase de modelos y la cardinalidad logarítmica mínima de un modelo en la clase que contiene los datos. La función de estructura determina todas las propiedades estocásticas de la cadena de datos individual: para cada clase de modelo restringida, determina el modelo individual que mejor se ajusta a la clase, independientemente de si el modelo verdadero está en la clase de modelos considerada o no. En el caso clásico, hablamos de un conjunto de datos con una distribución de probabilidad , y las propiedades son las de las esperanzas. En cambio, aquí tratamos con cadenas de datos individuales y las propiedades de la cadena individual en la que nos centramos. En este contexto, una propiedad se cumple con certeza en lugar de con alta probabilidad, como en el caso clásico. La función de estructura de Kolmogorov cuantifica con precisión la bondad de ajuste de un modelo individual con respecto a los datos individuales.

La función de estructura de Kolmogorov se utiliza en la teoría de la información algorítmica , también conocida como teoría de la complejidad de Kolmogorov, para describir la estructura de una cadena mediante el uso de modelos de complejidad creciente.

Definición de Kolmogorov

Kolmogorov (izquierda) habla sobre la función de estructura (ver dibujo en la pizarra) en ( Tallin , 1973).

La función de estructura fue propuesta originalmente por Kolmogorov en 1973 en un simposio soviético de Teoría de la Información en Tallin, pero estos resultados no fueron publicados [ 1 ] pág.  182. Sin embargo, los resultados fueron anunciados en [ 2 ] en 1974, el único registro escrito por el propio Kolmogorov. Una de sus últimas declaraciones científicas es (traducida del ruso original por L.A. Levin):

A cada objeto constructivo le corresponde una función.Φincógnita(k){\displaystyle \Phi _{x}(k)}de un número natural k: el logaritmo de la cardinalidad mínima de los conjuntos que contienen x y que permiten definiciones de complejidad como máximo k. Si el elemento x mismo permite una definición simple, entonces la funciónΦ{\displaystyle \Phi }cae a 0 incluso para k pequeño. Al carecer de tal definición, el elemento es "aleatorio" en un sentido negativo. Pero es "probabilísticamente aleatorio" positivamente solo cuando la funciónΦ{\displaystyle \Phi }habiendo tomado el valorΦ0{\displaystyle \Phi _{0}}a una escala relativamente pequeñak=k0{\displaystyle k=k_{0}}, luego cambia aproximadamente comoΦ(k)=Φ0(kk0){\displaystyle \Phi (k)=\Phi _{0}-(k-k_{0})}.

Kolmogorov , anuncio citado anteriormente

Definición contemporánea

Se discute en Cover y Thomas. [ 1 ] Se estudia extensamente en Vereshchagin y Vitányi [ 3 ] donde también se resuelven las propiedades principales. La función de estructura de Kolmogorov se puede escribir como hincógnita(α)=minS{registro|S|:incógnitaS,K(S)α}{\displaystyle h_{x}(\alpha )=\min _{S}\{\log |S|:x\in S,K(S)\leq \alpha \}} dóndeincógnita{\displaystyle x}es una cadena binaria de longitudnorte{\displaystyle n}conincógnitaS{\displaystyle x\in S}dóndeS{\displaystyle S} es un modelo contemplado (conjunto de cadenas de longitud n) paraincógnita{\displaystyle x},K(S){\displaystyle K(S)}es la complejidad de Kolmogorov deS{\displaystyle S}yα{\displaystyle \alpha } es un valor entero no negativo que limita la complejidad de lo contemplado.S{\displaystyle S}'s. Claramente, esta función no es creciente y alcanzaregistro|{incógnita}|=0{\displaystyle \log \left|\{x\}\right|=0}paraα=K(incógnita)+do{\displaystyle \alpha =K(x)+c}dóndedo{\displaystyle c}es el número de bits necesarios para cambiarincógnita{\displaystyle x}en{incógnita}{\displaystyle \{x\}}yK(incógnita){\displaystyle K(x)}es la complejidad de Kolmogorov deincógnita{\displaystyle x}.

La estadística suficiente algorítmica

Definimos un conjuntoS{\displaystyle S}que contieneincógnita{\displaystyle x}de tal manera que K(S)+K(incógnita|S)=K(incógnita)+O(1).{\displaystyle K(S)+K(x|S)=K(x)+O(1).} La funciónhincógnita(α){\displaystyle h_{x}(\alpha )}nunca disminuye más que una constante independiente fija por debajo de la diagonal llamada línea de suficiencia L definida por L(α)+α=K(incógnita).{\displaystyle L(\alpha )+\alpha =K(x).} Se aproxima a una distancia constante mediante la gráfica dehincógnita{\displaystyle h_{x}}por ciertos argumentos (por ejemplo, paraα=K(incógnita)+do{\displaystyle \alpha =K(x)+c}). Para estosα{\displaystyle \alpha }tenemosα+hincógnita(α)=K(incógnita)+O(1){\displaystyle \alpha +h_{x}(\alpha )=K(x)+O(1)}y el modelo asociadoS{\displaystyle S}(testigo dehincógnita(α){\displaystyle h_{x}(\alpha )}) se denomina un conjunto óptimo paraincógnita{\displaystyle x}y su descripción deK(S)α{\displaystyle K(S)\leq \alpha }bits es, por lo tanto, una estadística suficiente algorítmica . Por convención, escribimos «algorítmica» en lugar de «complejidad de Kolmogorov». Las principales propiedades de una estadística suficiente algorítmica son las siguientes: SiS{\displaystyle S}es una estadística suficiente algorítmica paraincógnita{\displaystyle x}, entonces K(S)+registro|S|=K(incógnita)+O(1).{\displaystyle K(S)+\log |S|=K(x)+O(1).} Es decir, la descripción en dos partes deincógnita{\displaystyle x}utilizando el modeloS{\displaystyle S}y como código de datos a modelo el índice deincógnita{\displaystyle x}en la enumeración deS{\displaystyle S}enregistro|S|{\displaystyle \log |S|}bits, es tan conciso como el código de una parte más corto deincógnita{\displaystyle x}enK(incógnita){\displaystyle K(x)}bits. Esto se puede ver fácilmente de la siguiente manera: K(incógnita)K(incógnita,S)+O(1)K(S)+K(incógnita|S)+O(1)K(S)+registro|S|+O(1)K(incógnita)+O(1),{\displaystyle {\begin{aligned}K(x)&\leq K(x,S)+O(1)\\&\leq K(S)+K(x|S)+O(1)\\&\leq K(S)+\log |S|+O(1)\\&\leq K(x)+O(1),\end{aligned}}}

Funciones de estructura
Funciones de estructurahincógnita(α),βincógnita(α),λincógnita(α){\displaystyle h_{x}(\alpha ),\beta _{x}(\alpha ),\lambda _{x}(\alpha )}y estadística mínima suficiente.

Utilizando desigualdades sencillas y la propiedad de suficiencia, encontramos queK(incógnita|S)=registro|S|+O(1){\displaystyle K(x|S)=\log |S|+O(1)}. (Por ejemplo, dadoSincógnita{\displaystyle S\ni x}, podemos describirincógnita{\displaystyle x}autolimitante (puedes determinar su final) enregistro|S|+O(1){\displaystyle \log |S|+O(1)}bits.) Por lo tanto, la deficiencia de aleatoriedadregistro|S|K(incógnita|S){\displaystyle \log |S|-K(x|S)}deincógnita{\displaystyle x}enS{\displaystyle S}es una constante, lo que significa queincógnita{\displaystyle x}es un elemento típico (aleatorio) de S. Sin embargo, puede haber modelosS{\displaystyle S}que contieneincógnita{\displaystyle x}que no son estadísticas suficientes. Una estadística suficiente algorítmicaS{\displaystyle S}paraincógnita{\displaystyle x}tiene la propiedad adicional, aparte de ser un modelo de mejor ajuste, de queK(incógnita,S)=K(incógnita)+O(1){\displaystyle K(x,S)=K(x)+O(1)}y por lo tanto por la simetría de complejidad de Kolmogorov de la información (la información sobreincógnita{\displaystyle x}enS{\displaystyle S}es aproximadamente lo mismo que la información sobreS{\displaystyle S}en x) tenemosK(S|incógnita)=O(1){\displaystyle K(S|x^{*})=O(1)}: la estadística suficiente algorítmicaS{\displaystyle S}es un modelo de mejor ajuste que está casi completamente determinado porincógnita{\displaystyle x}. (incógnita{\displaystyle x^{*}}es el programa más corto paraincógnita{\displaystyle x}.) La estadística suficiente algorítmica asociada con el menor de talesα{\displaystyle \alpha }se denomina estadística mínima suficiente algorítmica .

Con respecto a la imagen: La función de estructura MDLλincógnita(α){\displaystyle \lambda _{x}(\alpha )}Se explica a continuación. La función de estructura de bondad de ajusteβincógnita(α){\displaystyle \beta _{x}(\alpha )}es la menor deficiencia de aleatoriedad (ver arriba) de cualquier modeloSincógnita{\displaystyle S\ni x}paraincógnita{\displaystyle x}de tal manera queK(S)α{\displaystyle K(S)\leq \alpha }Esta función de estructura proporciona la bondad de ajuste de un modelo.S{\displaystyle S}(que contiene x) para la cadena x. Cuando es bajo, el modelo se ajusta bien, y cuando es alto, el modelo no se ajusta bien. Siβincógnita(α)=0{\displaystyle \beta _{x}(\alpha )=0}para algunosα{\displaystyle \alpha }Entonces existe un modelo típico.Sincógnita{\displaystyle S\ni x}paraincógnita{\displaystyle x}de tal manera queK(S)α{\displaystyle K(S)\leq \alpha }yincógnita{\displaystyle x}es típico (aleatorio) para S. Es decir,S{\displaystyle S}es el modelo que mejor se ajusta a x. Para más detalles, consulte [ 1 ] y especialmente [ 3 ] y [ 4 ] .

Selección de propiedades

Dentro de las restricciones de que la gráfica desciende en un ángulo de al menos 45 grados, que comienza en n y termina aproximadamente enK(incógnita){\displaystyle K(x)}, cada gráfico (hasta unO(registronorte){\displaystyle O(\log n)}El término aditivo en el argumento y el valor se realiza mediante la función de estructura de algunos datos x y viceversa. Donde el gráfico toca la diagonal primero, el argumento (complejidad) es el de la estadística suficiente mínima. Es incalculable determinar este lugar. Véase [ 3 ] .

Propiedad principal

Está demostrado que en cada nivelα{\displaystyle \alpha }de complejidad la función de estructura nos permite seleccionar el mejor modeloS{\displaystyle S}para la cadena individual x dentro de una tira deO(registronorte){\displaystyle O(\log n)}con certeza, no con gran probabilidad. [ 3 ]

La variante MDL

La función de longitud de descripción mínima (MDL): La longitud del código mínimo de dos partes para x que consta del costo del modelo K(S) y la longitud del índice de x en S, en la clase de modelos de conjuntos de complejidad de Kolmogorov máxima dadaα{\displaystyle \alpha }, la complejidad de S está limitada superiormente porα{\displaystyle \alpha }, viene dado por la función MDL o el estimador MDL restringido:

λincógnita(α)=minS{Λ(S):Sincógnita,K(S)α},{\displaystyle \lambda _{x}(\alpha )=\min _{S}\{\Lambda (S):S\ni x,\;K(S)\leq \alpha \},} dóndeΛ(S)=registro|S|+K(S)K(incógnita)O(1){\displaystyle \Lambda (S)=\log |S|+K(S)\geq K(x)-O(1)}es la longitud total del código de dos partes de x con la ayuda del modelo S.

Propiedad principal

Está demostrado que en cada nivelα{\displaystyle \alpha }de complejidad la función de estructura nos permite seleccionar el mejor modelo S para la cadena individual x dentro de una tira deO(registronorte){\displaystyle O(\log n)}con certeza, no con gran probabilidad. [ 3 ]

Aplicación en estadística

Las matemáticas desarrolladas anteriormente fueron tomadas como fundamento del MDL por su inventor Jorma Rissanen . [ 5 ]

Modelos de probabilidad

Para cada distribución de probabilidad computablePAG{\displaystyle P}Se puede demostrar [ 6 ] que registroPAG(incógnita)=registro|S|+O(registronorte).{\displaystyle -\log P(x)=\log |S|+O(\log n).} Por ejemplo, siPAG{\displaystyle P}es alguna distribución computable en el conjuntoS{\displaystyle S}de cadenas de longitudnorte{\displaystyle n}, luego cada unoincógnitaS{\displaystyle x\in S}tiene probabilidadPAG(incógnita)=exp(O(registronorte))/|S|=norteO(1)/|S|{\displaystyle P(x)=\exp(O(\log n))/|S|=n^{O(1)}/|S|}La función de estructura de Kolmogorov se convierte en hincógnita(α)=minPAG{registroPAG(incógnita):PAG(incógnita)>0,K(PAG)α}{\displaystyle h'_{x}(\alpha )=\min _{P}\left\{-\log P(x):P(x)>0,K(P)\leq \alpha \right\}} donde x es una cadena binaria de longitud n conregistroPAG(incógnita)>0{\displaystyle -\log P(x)>0}dóndePAG{\displaystyle P} es un modelo contemplado (probabilidad computable denorte{\displaystyle n}-cadenas de longitud) paraincógnita{\displaystyle x},K(PAG){\displaystyle K(P)}es la complejidad de Kolmogorov dePAG{\displaystyle P}yα{\displaystyle \alpha } es un valor entero que delimita la complejidad de lo contemplado.PAG{\displaystyle P}'s. Claramente, esta función no es creciente y alcanzaregistro|{incógnita}|=0{\displaystyle \log |\{x\}|=0}paraα=K(incógnita)+do{\displaystyle \alpha =K(x)+c}donde c es el número de bits necesarios para cambiarincógnita{\displaystyle x}en{incógnita}{\displaystyle \{x\}}yK(incógnita){\displaystyle K(x)}es la complejidad de Kolmogorov deincógnita{\displaystyle x}. Entonceshincógnita(α)=hincógnita(α)+O(registronorte){\displaystyle h'_{x}(\alpha )=h_{x}(\alpha )+O(\log n)}Para cada nivel de complejidadα{\displaystyle \alpha }la funciónhincógnita(α){\displaystyle h'_{x}(\alpha )} es la versión de complejidad de Kolmogorov de la máxima verosimilitud (ML).

Propiedad principal

Está demostrado que en cada nivelα{\displaystyle \alpha }de complejidad la función de estructura nos permite seleccionar el mejor modeloS{\displaystyle S}para la cadena individualincógnita{\displaystyle x}dentro de una franja deO(registronorte){\displaystyle O(\log n)}con certeza, no con gran probabilidad. [ 3 ]

La variante MDL y los modelos de probabilidad

La función MDL: La longitud del código mínimo de dos partes para x que consta del costo del modelo K(P) y la longitud deregistroPAG(incógnita){\displaystyle -\log P(x)}, en la clase modelo de funciones de masa de probabilidad computables de complejidad de Kolmogorov máxima dadaα{\displaystyle \alpha }, la complejidad de P está limitada superiormente porα{\displaystyle \alpha }, viene dado por la función MDL o el estimador MDL restringido:

λincógnita(α)=minPAG{Λ(PAG):PAG(incógnita)>0,K(PAG)α},{\displaystyle \lambda '_{x}(\alpha )=\min _{P}\{\Lambda (P):P(x)>0,\;K(P)\leq \alpha \},} dóndeΛ(PAG)=registroPAG(incógnita)+K(PAG)K(incógnita)O(1){\displaystyle \Lambda (P)=-\log P(x)+K(P)\geq K(x)-O(1)}es la longitud total del código de dos partes de x con la ayuda del modelo P.

Propiedad principal

Está demostrado que en cada nivelα{\displaystyle \alpha }de complejidad la función MDL nos permite seleccionar el mejor modelo P para la cadena individual x dentro de una tira deO(registronorte){\displaystyle O(\log n)}con certeza, no con gran probabilidad. [ 3 ]

Extensión a la distorsión de velocidad y eliminación de ruido

Resulta que el enfoque puede extenderse a una teoría de distorsión de tasa de secuencias finitas individuales y eliminación de ruido de secuencias finitas individuales [ 7 ] utilizando la complejidad de Kolmogorov. Se han realizado experimentos con éxito utilizando programas compresores reales [ 8 ] . Aquí se supone que, para datos naturales, la complejidad de Kolmogorov no está lejos de la longitud de una versión comprimida utilizando un buen compresor.

Referencias

  1. 1 2 3 Cover, Thomas M.; Thomas, Joy A. (1991). Elementos de la teoría de la información . Nueva York: Wiley. págs. 175–178 . ISBN  978-0471062592.
  2. Resumen de una charla para la Sociedad Matemática de Moscú en Uspekhi Mat. Nauk Volumen 29, Número 4(178) en las Comunicaciones de la Sociedad Matemática de Moscú página 155 (en la edición rusa, no traducida al inglés)
  3. 1 2 3 4 5 6 7 Vereshchagin, NK; Vitanyi, PMB (1 de diciembre de 2004). "Funciones estructurales de Kolmogorov y selección de modelos". IEEE Transactions on Information Theory . 50 (12): 3265– 3290. arXiv : cs/0204037 . doi : 10.1109/TIT.2004.838346 .
  4. Gacs, P.; Tromp, JT; Vitanyi, PMB (2001). "Algorithmic statistics". IEEE Transactions on Information Theory . 47 (6): 2443– 2463. arXiv : math/0006233 . doi : 10.1109/18.945257 .
  5. ^ Rissanen, Jorma (2007). Información y complejidad en el modelado estadístico (Online-Ausg. ed.). Nueva York: Springer. ISBN  978-0-387-36610-4.
  6. A.Kh. Shen, El concepto de estocasticidad ( α , β ) en el sentido de Kolmogorov y sus propiedades, Soviet Math. Dokl., 28:1(1983), 295-299
  7. Vereshchagin, Nikolai K.; Vitanyi, Paul MB (1 de julio de 2010). "Distorsión de tasa y eliminación de ruido de datos individuales mediante la complejidad de Kolmogorov". IEEE Transactions on Information Theory . 56 (7): 3438– 3454. arXiv : cs/0411014 . doi : 10.1109/TIT.2010.2048491 .
  8. de Rooij, Steven; Vitanyi, Paul (1 de marzo de 2012). "Aproximación de gráficos de tasa-distorsión de datos individuales: experimentos en compresión con pérdidas y eliminación de ruido". IEEE Transactions on Computers . 61 (3): 395– 407. arXiv : cs/0609121 . doi : 10.1109/TC.2011.25 .

Literatura

  • Cover, TM; P. Gacs; RM Gray (1989). "Contribuciones de Kolmogorov a la teoría de la información y la complejidad algorítmica" . Annals of Probability . 17 (3): 840– 865. doi : 10.1214/aop/1176991250 . JSTOR 2244387 . 
  • Kolmogorov, AN; Uspenskii, VA (1 de enero de 1987). "Algoritmos y aleatoriedad" . Teoría de la probabilidad y sus aplicaciones . 32 (3): 389– 412. doi : 10.1137/1132060 .
  • Li, M., Vitányi, PMB (2008). Introducción a la complejidad de Kolmogorov y sus aplicaciones (3.ª  ed.). Nueva York: Springer. ISBN 978-0387339986.{{cite book}}: CS1 maint: nombres múltiples: lista de autores ( enlace ) , especialmente las páginas  401-431 sobre la función de estructura de Kolmogorov y las páginas  613-629 sobre la distorsión de la tasa y la eliminación de ruido de secuencias individuales.
  • Shen, A. (1 de abril de 1999). "Discusión sobre la complejidad de Kolmogorov y el análisis estadístico". The Computer Journal . 42 (4): 340– 342. doi : 10.1093/comjnl/42.4.340 .
  • V'yugin, VV (1987). "Sobre el defecto de aleatoriedad de un objeto finito en relación con medidas con límites de complejidad dados" . Teoría de la probabilidad y sus aplicaciones . 32 (3): 508– 512. doi : 10.1137/1132071 .
  • V'yugin, VV (1 de abril de 1999). "Complejidad algorítmica y propiedades estocásticas de secuencias binarias finitas". The Computer Journal . 42 (4): 294– 317. doi : 10.1093/comjnl/42.4.294 .