Articulo de referencia

cadena de Markov

Diagrama que representa un proceso de Markov de dos estados. Los números indican la probabilidad de pasar de un estado a otro. En teoría de la probabilidad y estadística , una c...

Diagrama que representa un proceso de Markov de dos estados. Los números indican la probabilidad de pasar de un estado a otro.

En teoría de la probabilidad y estadística , una cadena de Markov o proceso de Markov es un proceso estocástico que describe una secuencia de eventos posibles en la que la probabilidad de cada evento depende únicamente del estado alcanzado en el evento anterior. De manera informal, esto puede entenderse como: "Lo que sucede a continuación depende únicamente del estado actual de las cosas ". Una secuencia infinita numerable , en la que la cadena cambia de estado en pasos de tiempo discretos, da lugar a una cadena de Markov de tiempo discreto (CMTD). Un proceso de tiempo continuo se denomina cadena de Markov de tiempo continuo (CMTC). Los procesos de Markov reciben su nombre en honor al matemático ruso Andrey Markov .

Las cadenas de Markov tienen muchas aplicaciones como modelos estadísticos de procesos del mundo real. [ 1 ] Proporcionan la base para métodos generales de simulación estocástica conocidos como Monte Carlo de cadena de Markov , que se utilizan para simular el muestreo de distribuciones de probabilidad complejas y han encontrado aplicación en áreas que incluyen estadística bayesiana , biología , química , economía , finanzas , teoría de la información , física , procesamiento de señales y procesamiento del habla . [ 1 ] [ 2 ] [ 3 ]

Los adjetivos markoviano y markoviano se utilizan para describir algo que está relacionado con un proceso de Markov. [ 4 ]

Principios

El matemático ruso Andrey Markov

Definición

Un proceso de Markov es un proceso estocástico que satisface la propiedad de Markov (a veces caracterizada como " falta de memoria "). En términos más sencillos, es un proceso para el cual se pueden hacer predicciones sobre resultados futuros basándose únicamente en su estado presente y, lo que es más importante, dichas predicciones son tan buenas como las que se podrían hacer conociendo el historial completo del proceso. [ 5 ] En otras palabras, condicionado al estado presente del sistema, sus estados pasados ​​y futuros son independientes .

Una cadena de Markov es un tipo de proceso de Markov que tiene un espacio de estados discreto o un conjunto de índices discretos (que a menudo representan el tiempo), pero la definición precisa de una cadena de Markov varía. [ 6 ] Por ejemplo, es común definir una cadena de Markov como un proceso de Markov en tiempo discreto o continuo con un espacio de estados numerable (por lo tanto, independientemente de la naturaleza del tiempo), [ 7 ] [ 8 ] [ 9 ] [ 10 ] pero también es común definir una cadena de Markov como que tiene tiempo discreto en un espacio de estados numerable o continuo (por lo tanto, independientemente del espacio de estados). [ 6 ]

Tipos de cadenas de Markov

Es necesario especificar el espacio de estados del sistema y el índice del parámetro de tiempo. La siguiente tabla ofrece una visión general de las diferentes instancias de procesos de Markov para distintos niveles de generalidad del espacio de estados, tanto para tiempo discreto como continuo:

Cabe señalar que no existe un acuerdo definitivo en la literatura sobre el uso de algunos términos que designan casos especiales de procesos de Markov. Generalmente, el término "cadena de Markov" se reserva para un proceso con un conjunto discreto de tiempos, es decir, una cadena de Markov de tiempo discreto (CMTD) [ 11 ] , pero algunos autores utilizan el término "proceso de Markov" para referirse a una cadena de Markov de tiempo continuo (CMTC) sin mención explícita [ 12 ] [ 13 ] [ 14 ] . Además, existen otras extensiones de procesos de Markov que se denominan así, pero que no necesariamente se incluyen en ninguna de estas cuatro categorías (véase modelo de Markov ). Asimismo, el índice de tiempo no tiene por qué ser necesariamente de valor real; al igual que con el espacio de estados, existen procesos concebibles que se mueven a través de conjuntos de índices con otras construcciones matemáticas. Nótese que la cadena de Markov de tiempo continuo con espacio de estados general es tan general que no tiene un término designado.

Si bien el parámetro de tiempo suele ser discreto, el espacio de estados de una cadena de Markov no tiene restricciones generalmente aceptadas: el término puede referirse a un proceso en un espacio de estados arbitrario. [ 15 ] Sin embargo, muchas aplicaciones de las cadenas de Markov emplean espacios de estados finitos o infinitos numerables , que presentan un análisis estadístico más directo. Además de los parámetros de índice de tiempo y espacio de estados, existen muchas otras variaciones, extensiones y generalizaciones (véase Variaciones ). Para simplificar, la mayor parte de este artículo se centra en el caso de tiempo discreto y espacio de estados discreto, salvo que se indique lo contrario.

Transiciones

Los cambios de estado del sistema se denominan transiciones. Las probabilidades asociadas a los distintos cambios de estado se denominan probabilidades de transición. El proceso se caracteriza por un espacio de estados, una matriz de transición que describe las probabilidades de transiciones específicas y un estado inicial (o distribución inicial) en dicho espacio. Por convención, se asume que todos los estados y transiciones posibles se han incluido en la definición del proceso, por lo que siempre existe un estado siguiente y el proceso no finaliza.

Un proceso aleatorio de tiempo discreto implica un sistema que se encuentra en un estado determinado en cada paso, y dicho estado cambia aleatoriamente entre pasos. Los pasos suelen considerarse como instantes en el tiempo, pero también pueden referirse a distancias físicas o cualquier otra medida discreta. Formalmente, los pasos son los números enteros o naturales , y el proceso aleatorio es una correspondencia entre estos y los estados. La propiedad de Markov establece que la distribución de probabilidad condicional del sistema en el siguiente paso (y, de hecho, en todos los pasos futuros) depende únicamente del estado actual del sistema, y ​​no del estado del sistema en pasos anteriores.

Dado que el sistema cambia aleatoriamente, generalmente es imposible predecir con certeza el estado de una cadena de Markov en un momento dado del futuro. Sin embargo, sí se pueden predecir las propiedades estadísticas del futuro del sistema. En muchas aplicaciones, son precisamente estas propiedades estadísticas las que resultan importantes.

Historia

Andrey Markov estudió los procesos de Markov a principios del siglo XX, publicando su primer artículo sobre el tema en 1906. [ 16 ] [ 17 ] [ 18 ] Los procesos de Markov en tiempo continuo fueron descubiertos mucho antes de su trabajo a principios del siglo XX en forma del proceso de Poisson . [ 19 ] [ 20 ] [ 21 ] Markov estaba interesado en estudiar una extensión de secuencias aleatorias independientes, motivado por un desacuerdo con Pavel Nekrasov, quien afirmaba que la independencia era necesaria para que se cumpliera la ley débil de los grandes números . [ 22 ] En su primer artículo sobre cadenas de Markov, publicado en 1906, Markov demostró que bajo ciertas condiciones los resultados promedio de la cadena de Markov convergerían a un vector fijo de valores, probando así una ley débil de los grandes números sin la suposición de independencia, [ 16 ] [ 17 ] [ 18 ] que se había considerado comúnmente como un requisito para que se cumplieran tales leyes matemáticas. [ 18 ] Más tarde, Markov utilizó cadenas de Markov para estudiar la distribución de vocales en Eugenio Oneguin , escrito por Alexander Pushkin , y demostró un teorema del límite central para dichas cadenas. [ 16 ]

En 1912, Henri Poincaré estudió cadenas de Markov en grupos finitos con el objetivo de estudiar el barajado de cartas. Otros usos tempranos de las cadenas de Markov incluyen un modelo de difusión, introducido por Paul y Tatyana Ehrenfest en 1907, y un proceso de ramificación, introducido por Francis Galton y Henry William Watson en 1873, anterior al trabajo de Markov. [ 16 ] [ 17 ] Después del trabajo de Galton y Watson, se reveló más tarde que su proceso de ramificación había sido descubierto y estudiado independientemente unas tres décadas antes por Irénée-Jules Bienaymé . [ 23 ] A partir de 1928, Maurice Fréchet se interesó en las cadenas de Markov, lo que finalmente lo llevó a publicar en 1938 un estudio detallado sobre las mismas. [ 16 ] [ 24 ]

Andrey Kolmogorov desarrolló en un artículo de 1931 gran parte de la teoría temprana de los procesos de Markov de tiempo continuo. [ 25 ] [ 26 ] Kolmogorov se inspiró en parte en el trabajo de Louis Bachelier de 1900 sobre las fluctuaciones en el mercado de valores, así como en el trabajo de Norbert Wiener sobre el modelo de Einstein del movimiento browniano. [ 25 ] [ 27 ] Introdujo y estudió un conjunto particular de procesos de Markov conocidos como procesos de difusión, donde derivó un conjunto de ecuaciones diferenciales que describen los procesos. [ 25 ] [ 28 ] Independientemente del trabajo de Kolmogorov, Sydney Chapman derivó en un artículo de 1928 una ecuación, ahora llamada ecuación de Chapman-Kolmogorov , de una manera menos rigurosa matemáticamente que Kolmogorov, mientras estudiaba el movimiento browniano. [ 29 ] Las ecuaciones diferenciales ahora se llaman ecuaciones de Kolmogorov [ 30 ] o ecuaciones de Kolmogorov-Chapman. [ 31 ] Otros matemáticos que contribuyeron significativamente a los fundamentos de los procesos de Markov incluyen a William Feller , a partir de la década de 1930, y luego Eugene Dynkin , a partir de la década de 1950. [ 26 ]

Ejemplos

  • Mark V. Shaney es un programa de cadena de Markov de tercer orden y un generador de texto de Markov . Ingiere el texto de muestra (el Tao Te Ching o las publicaciones de un grupo de Usenet ) y crea una lista masiva de todas las secuencias de tres palabras consecutivas (tripletes) que aparecen en el texto. Luego, elige dos palabras al azar y busca una palabra que siga a esas dos en uno de los tripletes de su lista masiva. Si hay más de una, elige una al azar (los tripletes idénticos se cuentan por separado, por lo que una secuencia que aparece dos veces tiene el doble de probabilidades de ser elegida que una que aparece solo una vez). Luego, agrega esa palabra al texto generado. A continuación, de la misma manera, elige un triplete que comienza con la segunda y la tercera palabra del texto generado, lo que da como resultado una cuarta palabra. Agrega la cuarta palabra, luego repite con la tercera y la cuarta palabra, y así sucesivamente. [ 32 ]
  • Los paseos aleatorios basados ​​en enteros y el problema de la ruina del jugador son ejemplos de procesos de Markov. [ 33 ] [ 34 ] Algunas variaciones de estos procesos se estudiaron cientos de años antes en el contexto de variables independientes. [ 35 ] [ 36 ] Dos ejemplos importantes de procesos de Markov son el proceso de Wiener , también conocido como proceso de movimiento browniano , y el proceso de Poisson , [ 19 ] que se consideran los procesos estocásticos más importantes y centrales en la teoría de procesos estocásticos. [ 37 ] [ 38 ] [ 39 ] Estos dos procesos son procesos de Markov en tiempo continuo, mientras que los paseos aleatorios sobre los enteros y el problema de la ruina del jugador son ejemplos de procesos de Markov en tiempo discreto. [ 33 ] [ 34 ]
  • Una cadena de Markov muy conocida es el llamado "paseo del borracho", un paseo aleatorio en la recta numérica donde, en cada paso, la posición puede cambiar en +1 o −1 con igual probabilidad. Desde cualquier posición, existen dos transiciones posibles: al siguiente o al anterior entero. Las probabilidades de transición dependen únicamente de la posición actual, no de cómo se alcanzó dicha posición. Por ejemplo, las probabilidades de transición de 5 a 4 y de 5 a 6 son ambas de 0,5, y todas las demás probabilidades de transición desde 5 son 0. Estas probabilidades son independientes de si el sistema se encontraba previamente en 4 o 6.
  • Una serie de estados independientes (por ejemplo, una serie de lanzamientos de moneda) satisface la definición formal de una cadena de Markov. Sin embargo, la teoría se suele aplicar solo cuando la distribución de probabilidad del siguiente estado depende del estado actual.

Un ejemplo no markoviano

Supongamos que hay un monedero que contiene cinco monedas de 25¢ (cuartos de dólar), cinco monedas de 10¢ (diez centavos) y cinco monedas de 5¢ (cinco centavos). Una por una, se extraen monedas al azar del monedero y se colocan sobre una mesa. Siincógnitanorte{\displaystyle X_{n}}representa el valor total de las monedas colocadas sobre la mesa después de n extracciones, conincógnita0=0{\displaystyle X_{0}=0}, luego la secuencia{incógnitanorte:nortenorte}{\displaystyle \{X_{n}:n\in \mathbb {N} \}}no es un proceso de Markov.

Para ver por qué sucede esto, supongamos que en los primeros seis sorteos se extraen las cinco monedas de cinco centavos y una de veinticinco centavos. Por lo tanto,incógnita6=$0,50{\displaystyle X_{6}=\$0.50}Si no sabemos nada másincógnita6{\displaystyle X_{6}}, pero también los valores anteriores, entonces podemos determinar qué monedas se han sacado, y sabemos que la siguiente moneda no será un níquel; por lo tanto, podemos determinar queincógnita7$0,60{\displaystyle X_{7}\geq \$0.60}con probabilidad 1. Pero si no conocemos los valores anteriores, entonces basándonos únicamente en el valorincógnita6{\displaystyle X_{6}}podríamos suponer que sacamos cuatro monedas de diez centavos y dos de cinco centavos, en cuyo caso ciertamente sería posible sacar otra moneda de cinco centavos a continuación. Por lo tanto, nuestras suposiciones sobreincógnita7{\displaystyle X_{7}}se ven afectados por nuestro conocimiento de los valores anteriores aincógnita6{\displaystyle X_{6}}.

Sin embargo, es posible modelar este escenario como un proceso de Markov. En lugar de definirincógnitanorte{\displaystyle X_{n}}Para representar el valor total de las monedas sobre la mesa, podríamos definirincógnitanorte{\displaystyle X_{n}}para representar la cantidad de los distintos tipos de monedas en la mesa. Por ejemplo,incógnita6=1,0,5{\displaystyle X_{6}=1,0,5}podría definirse para representar el estado en el que hay una moneda de veinticinco centavos, cero monedas de diez centavos y cinco monedas de cinco centavos sobre la mesa después de 6 extracciones individuales. Este nuevo modelo podría representarse mediante6×6×6=216{\displaystyle 6\times 6\times 6=216}Estados posibles, donde cada estado representa la cantidad de monedas de cada tipo (de 0 a 5) que hay sobre la mesa. (No todos estos estados se pueden alcanzar en 6 extracciones).

Supongamos que el primer sorteo resulta en el estadoincógnita1=0,1,0{\displaystyle X_{1}=0,1,0}. La probabilidad de lograrincógnita2{\displaystyle X_{2}}ahora depende deincógnita1{\displaystyle X_{1}}; por ejemplo, el estadoincógnita2=1,0,1{\displaystyle X_{2}=1,0,1}no es posible. Después del segundo sorteo, el tercer sorteo depende de las monedas que se hayan extraído hasta el momento, pero ya no solo de las monedas que se extrajeron para el primer estado (ya que se ha añadido información probabilísticamente importante al escenario). De esta manera, la probabilidad de laincógnitanorte=i,j,k{\displaystyle X_{n}=i,j,k}El estado depende exclusivamente del resultado de laincógnitanorte1=,metro,pag{\displaystyle X_{n-1}=\ell ,m,p}estado.

Definición formal

Cadena de Markov de tiempo discreto

Una cadena de Markov de tiempo discreto es una secuencia de variables aleatorias X 1 , X 2 , X 3 , ... con la propiedad de Markov , es decir, que la probabilidad de pasar al siguiente estado depende solo del estado actual y no de los estados anteriores:

Pr(incógnitanorte+1=incógnitaincógnita1=incógnita1,incógnita2=incógnita2,,incógnitanorte=incógnitanorte)=Pr(incógnitanorte+1=incógnitaincógnitanorte=incógnitanorte),{\displaystyle \Pr(X_{n+1}=x\mid X_{1}=x_{1},X_{2}=x_{2},\ldots ,X_{n}=x_{n})=\Pr(X_{n+1}=x\mid X_{n}=x_{n}),}si ambas probabilidades condicionales están bien definidas, es decir, siPr(incógnita1=incógnita1,,incógnitanorte=incógnitanorte)>0.{\displaystyle \Pr(X_{1}=x_{1},\ldots ,X_{n}=x_{n})>0.}

Los posibles valores de X i forman un conjunto numerable S llamado espacio de estados de la cadena.

Variaciones

  • Las cadenas de Markov homogéneas en el tiempo son procesos dondePr(incógnitanorte+1=incógnitaincógnitanorte=y)=Pr(incógnitanorte=incógnitaincógnitanorte1=y){\displaystyle \Pr(X_{n+1}=x\mid X_{n}=y)=\Pr(X_{n}=x\mid X_{n-1}=y)}para todo n . La probabilidad de la transición es independiente de n .
  • Las cadenas de Markov estacionarias son procesos dondePr(incógnita0=incógnita0,incógnita1=incógnita1,,incógnitak=incógnitak)=Pr(incógnitanorte=incógnita0,incógnitanorte+1=incógnita1,,incógnitanorte+k=incógnitak){\displaystyle \Pr(X_{0}=x_{0},X_{1}=x_{1},\ldots ,X_{k}=x_{k})=\Pr(X_{n}=x_{0},X_{n+1}=x_{1},\ldots ,X_{n+k}=x_{k})}para todo n y k . Se puede demostrar que toda cadena estacionaria es homogénea en el tiempo mediante la regla de Bayes.
    Una condición necesaria y suficiente para que una cadena de Markov homogénea en el tiempo sea estacionaria es que la distribución deincógnita0{\displaystyle X_{0}}es una distribución estacionaria de la cadena de Markov.
  • Una cadena de Markov con memoria (o una cadena de Markov de orden m ), donde m es finito, es un proceso que satisfacePr(incógnitanorte=incógnitanorteincógnitanorte1=incógnitanorte1,incógnitanorte2=incógnitanorte2,,incógnita1=incógnita1)=Pr(incógnitanorte=incógnitanorteincógnitanorte1=incógnitanorte1,incógnitanorte2=incógnitanorte2,,incógnitanortemetro=incógnitanortemetro) para norte>metro{\displaystyle {\begin{aligned}{}&\Pr(X_{n}=x_{n}\mid X_{n-1}=x_{n-1},X_{n-2}=x_{n-2},\dots ,X_{1}=x_{1})\\=&\Pr(X_{n}=x_{n}\mid X_{n-1}=x_{n-1},X_{n-2}=x_{n-2},\dots ,X_{nm}=x_{nm}){\text{ para }}n>m\end{aligned}}}En otras palabras, el estado futuro depende de los m estados pasados. Es posible construir una cadena(Ynorte){\displaystyle (Y_{n})}de(incógnitanorte){\displaystyle (X_{n})}que tiene la propiedad de Markov 'clásica' al tomar como espacio de estados las m -tuplas ordenadas de valores X , es decir,Ynorte=(incógnitanorte,incógnitanorte1,,incógnitanortemetro+1){\displaystyle Y_{n}=\left(X_{n},X_{n-1},\ldots ,X_{n-m+1}\right)}.

Espacio de estados finito

Si el espacio de estados es finito , la distribución de probabilidad de transición se puede representar mediante una matriz , llamada matriz de transición, con el elemento ( i , j ) de P igual a

pagij=Pr(incógnitanorte+1=jincógnitanorte=i).{\displaystyle p_{ij}=\Pr(X_{n+1}=j\mid X_{n}=i).}

Dado que la suma de cada fila de P es igual a uno y todos sus elementos son no negativos, P es una matriz estocástica derecha .

Relación de la distribución estacionaria con los autovectores y símplices

Una distribución estacionaria π es un vector (fila) cuyas entradas son no negativas y suman 1, no se ve afectada por la operación de la matriz de transición P sobre él y, por lo tanto, se define por

πPAG=π.{\displaystyle \pi \mathbf {P} =\pi .}

Al comparar esta definición con la de un vector propio, vemos que los dos conceptos están relacionados y que

π=miimii{\displaystyle \pi ={\frac {e}{\sum _{i}{e_{i}}}}}

es un normalizado (iπi=1{\textstyle \sum _ {i} \pi _ {i} = 1}) múltiplo de un vector propio izquierdo e de la matriz de transición P con un valor propio de 1. Si hay más de un vector propio unitario, entonces una suma ponderada de los estados estacionarios correspondientes también es un estado estacionario. Pero para una cadena de Markov, generalmente se está más interesado en un estado estacionario que sea el límite de la secuencia de distribuciones para alguna distribución inicial.

Los valores de una distribución estacionariaπi{\displaystyle \textstyle \pi _ {i}}están asociados con el espacio de estados de P y sus autovectores conservan sus proporciones relativas. Dado que los componentes de π son positivos y la restricción de que su suma sea la unidad se puede reescribir comoi1πi=1{\textstyle \sum _{i}1\cdot \pi _{i}=1}vemos que el producto escalar de π con un vector cuyas componentes son todas 1 es la unidad y que π se encuentra en un simplex .

Cadena de Markov homogénea en el tiempo con un espacio de estados finito

Si la cadena de Markov es homogénea en el tiempo, entonces la matriz de transición P es la misma después de cada paso, por lo que la probabilidad de transición de k pasos se puede calcular como la k -ésima potencia de la matriz de transición, P k .

Si la cadena de Markov es irreducible y aperiódica, entonces existe una única distribución estacionaria π . [ 40 ] Además, en este caso P k converge a una matriz de rango uno en la que cada fila es la distribución estacionaria π :

límitekPAGk=1π{\displaystyle \lim _{k\to \infty }\mathbf {P} ^{k}=\mathbf {1} \pi }

donde 1 es el vector columna con todas las entradas iguales a 1. Esto se establece en el teorema de Perron-Frobenius . Si, por cualquier medio,límitekPAGk{\textstyle \lim _{k\to \infty }\mathbf {P} ^{k}}Una vez encontrada, la distribución estacionaria de la cadena de Markov en cuestión se puede determinar fácilmente para cualquier distribución inicial, como se explicará a continuación.

Para algunas matrices estocásticas P , el límitelímitekPAGk{\textstyle \lim _{k\to \infty }\mathbf {P} ^{k}}no existe, mientras que la distribución estacionaria sí, como lo demuestra este ejemplo:

PAG=(0110)PAG2k=IPAG2k+1=PAG{\displaystyle \mathbf {P} ={\begin{pmatrix}0&1\\1&0\end{pmatrix}}\qquad \mathbf {P} ^{2k}=I\qquad \mathbf {P} ^{2k+1}=\mathbf {P} }
(1212)(0110)=(1212){\displaystyle {\begin{pmatrix}{\frac {1}{2}}&{\frac {1}{2}}\end{pmatrix}}{\begin{pmatrix}0&1\\1&0\end{pmatrix}}={\begin{pmatrix}{\frac {1}{2}}&{\frac {1}{2}}\end{pmatrix}}}

(Este ejemplo ilustra una cadena de Markov periódica).

Debido a que hay varios casos especiales diferentes a considerar, el proceso de encontrar este límite, si existe, puede ser una tarea larga. Sin embargo, hay muchas técnicas que pueden ayudar a encontrar este límite. Sea P una matriz n × n y definamosQ=límitekPAGk.{\textstyle \mathbf {Q} =\lim _{k\to \infty }\mathbf {P} ^{k}.}

Siempre es cierto que

QPAG=Q.{\displaystyle \mathbf {QP} =\mathbf {Q}.}

Restando Q de ambos lados y factorizando se obtiene

Q(PAGInorte)=0norte,norte,{\displaystyle \mathbf {Q} (\mathbf {P} -\mathbf {I} _{n})=\mathbf {0} _{n,n},}

donde I n es la matriz identidad de tamaño n , y 0 n , n es la matriz cero de tamaño n × n . Multiplicar matrices estocásticas siempre produce otra matriz estocástica, por lo que Q debe ser una matriz estocástica (ver la definición anterior). A veces es suficiente usar la ecuación matricial anterior y el hecho de que Q es una matriz estocástica para resolver Q. Incluyendo el hecho de que la suma de cada una de las filas en P es 1, hay n+1 ecuaciones para determinar n incógnitas, por lo que es computacionalmente más fácil si por un lado se selecciona una fila en Q y se sustituye cada uno de sus elementos por uno, y por otro se sustituye el elemento correspondiente (el de la misma columna) en el vector 0 , y luego se multiplica por la izquierda este último vector por el inverso de la matriz anterior transformada para encontrar Q.

Aquí hay un método para hacerlo: primero, defina la función f ( A ) para devolver la matriz A con su columna más a la derecha reemplazada con todos 1. Si [ f ( PI n )] −1 existe entonces [ 41 ] [ 40 ]

Q=F(0norte,norte)[F(PAGInorte)]1.{\displaystyle \mathbf {Q} =f(\mathbf {0} _{n,n})[f(\mathbf {P} -\mathbf {I} _{n})]^{-1}.}
Explicación: La ecuación matricial original es equivalente a un sistema de n × n ecuaciones lineales con n × n variables. Además, existen n ecuaciones lineales adicionales debido a que Q es una matriz estocástica derecha cuyas filas suman 1. Por lo tanto, se necesitan n × n ecuaciones lineales independientes de las ( n × n + n ) ecuaciones para resolver las n × n variables. En este ejemplo, las n ecuaciones de " Q multiplicada por la columna más a la derecha de ( P - I n )" se han reemplazado por las n ecuaciones estocásticas.

Un aspecto a tener en cuenta es que si P tiene un elemento P i , i en su diagonal principal que es igual a 1 y la i -ésima fila o columna está llena de ceros, entonces esa fila o columna permanecerá sin cambios en todas las potencias subsiguientes P k . Por lo tanto, la i -ésima fila o columna de Q tendrá los 1 y los 0 en las mismas posiciones que en P .

Velocidad de convergencia a la distribución estacionaria

Como se indicó anteriormente, a partir de la ecuaciónπ=πPAG,{\displaystyle {\boldsymbol {\pi }}={\boldsymbol {\pi }}\mathbf {P},}Si existe, la distribución estacionaria (o de estado estable) π es un vector propio izquierdo de la matriz estocástica por filas P. Entonces, suponiendo que P es diagonalizable o, equivalentemente, que P tiene n vectores propios linealmente independientes, la velocidad de convergencia se desarrolla de la siguiente manera. (Para matrices no diagonalizables, es decir, matrices defectuosas , se puede comenzar con la forma normal de Jordan de P y proceder con un conjunto de argumentos un poco más complejos de manera similar. [ 42 ] )

Sea U la matriz de autovectores (cada uno normalizado para tener una norma L2 igual a 1) donde cada columna es un autovector izquierdo de P y sea Σ la matriz diagonal de autovalores izquierdos de P , es decir, Σ = diag( λ 1 , λ 2 , λ 3 ,..., λ n ). Entonces, por descomposición en autovalores

PAG=UΣU1.{\displaystyle \mathbf {P} =\mathbf {U\Sigma U} ^{-1}.}

Sean los autovalores enumerados de tal manera que:

1=|λ1|>|λ2||λ3||λnorte|.{\displaystyle 1=|\lambda _ {1}|>|\lambda _ {2}|\geq |\lambda _ {3}|\geq \cdots \geq |\lambda _ {n}|.}

Dado que P es una matriz estocástica por filas, su mayor valor propio izquierdo es 1. Si existe una distribución estacionaria única, entonces el mayor valor propio y el vector propio correspondiente también son únicos (porque no hay otro π que resuelva la ecuación de distribución estacionaria anterior). Sea u i la i -ésima columna de la matriz U , es decir, u i es el vector propio izquierdo de P correspondiente a λ i . Además, sea x un vector fila de longitud n que representa una distribución de probabilidad válida; dado que los vectores propios u i abarcanRnorte,{\displaystyle \mathbb {R} ^{n},}podemos escribir

incógnitaT=i=1norteaii,aiR.{\displaystyle \mathbf {x} ^{\mathsf {T}}=\sum _{i=1}^{n}a_{i}\mathbf {u} _{i},\qquad a_{i}\in \mathbb {R} .}

Si multiplicamos x por P desde la derecha y continuamos esta operación con los resultados, al final obtenemos la distribución estacionaria π . En otras palabras, π = a 1 u 1xPP ... P = xP k cuando k → ∞. Eso significa

π(k)=incógnita(UΣU1)(UΣU1)(UΣU1)=incógnitaUΣkU1=(a11T+a22T++anortenorteT)UΣkU1=a1λ1k1T+a2λ2k2T++anorteλnorteknorteTij para ij=λ1k{a11T+a2(λ2λ1)k2T+a3(λ3λ1)k3T++anorte(λnorteλ1)knorteT}{\displaystyle {\begin{aligned}{\boldsymbol {\pi }}^{(k)}&=\mathbf {x} \left(\mathbf {U\Sigma U} ^{-1}\right)\left(\mathbf {U\Sigma U} ^{-1}\right)\cdots \left(\mathbf {U\Sigma U} ^{-1}\right)\\&=\mathbf {xU\Sigma } ^{k}\mathbf {U} ^{-1}\\&=\left(a_{1}\mathbf {u} _{1}^{\mathsf {T}}+a_{2}\mathbf {u} _{2}^{\mathsf {T}}+\cdots +a_{n}\mathbf {u} _{n}^{\mathsf {T}}\right)\mathbf {U\Sigma } ^{k}\mathbf {U} ^{-1}\\&=a_{1}\lambda _{1}^{k}\mathbf {u} _{1}^{\mathsf {T}}+a_{2}\lambda _{2}^{k}\mathbf {u} _{2}^{\mathsf {T}}+\cdots +a_{n}\lambda _{n}^{k}\mathbf {u} _{n}^{\mathsf {T}}&&u_{i}\bot u_{j}{\text{ for }}i\neq j\\&=\lambda _{1}^{k}\left\{a_{1}\mathbf {u} _{1}^{\mathsf {T}}+a_{2}\left({\frac {\lambda _{2}}{\lambda _{1}}}\right)^{k}\mathbf {u} _{2}^{\mathsf {T}}+a_{3}\left({\frac {\lambda _{3}}{\lambda _{1}}}\right)^{k}\mathbf {u} _{3}^{\mathsf {T}}+\cdots +a_{n}\left({\frac {\lambda _{n}}{\lambda _{1}}}\right)^{k}\mathbf {u} _{n}^{\mathsf {T}}\right\}\end{aligned}}}

Dado que π es paralelo a u 1 (normalizado por la norma L2) y π ( k ) es un vector de probabilidad, π ( k ) se aproxima a a 1 u 1 = π cuando k → ∞ con una velocidad del orden de λ 2 / λ 1 exponencialmente. Esto se deduce de que|λ2||λnorte|,{\displaystyle |\lambda _{2}|\geq \cdots \geq |\lambda _{n}|,}Por lo tanto, λ 2 / λ 1 es el término dominante. Cuanto menor sea la razón, más rápida será la convergencia. [ 43 ] El ruido aleatorio en la distribución de estado π también puede acelerar esta convergencia a la distribución estacionaria. [ 44 ]

Cadena de Markov en tiempo continuo

Una cadena de Markov de tiempo continuo(incógnitat)t0{\displaystyle (X_{t})_{t\geq 0}}Se define mediante un espacio de estados finito o numerable S , una matriz de tasas de transición Q con dimensiones iguales a las del espacio de estados y una distribución de probabilidad inicial definida en dicho espacio. Para i j , los elementos q ij son no negativos y describen la tasa de transición del proceso del estado i al estado j . Los elementos q ii se eligen de tal manera que la suma de cada fila de la matriz de tasas de transición sea cero, mientras que las sumas de las filas de una matriz de probabilidad de transición en una cadena de Markov (discreta) son todas iguales a uno. 

Hay tres definiciones equivalentes del proceso. [ 45 ]

Definición infinitesimal

La cadena de Markov de tiempo continuo se caracteriza por las tasas de transición, las derivadas con respecto al tiempo de las probabilidades de transición entre los estados i y j.

Dejarincógnitat{\displaystyle X_{t}}Sea la variable aleatoria que describe el estado del proceso en el instante t , y supongamos que el proceso se encuentra en un estado i en el instante t . Entonces, sabiendoincógnitat=i{\displaystyle X_{t}=i},incógnitat+h=j{\displaystyle X_{t+h}=j}es independiente de los valores anteriores(incógnitas:s<t){\displaystyle \left(X_{s}:s<t\right)}y cuando h → 0 para todo j y para todo t , Pr(incógnita(t+h)=jincógnita(t)=i)=δij+qijh+o(h),{\displaystyle \Pr(X(t+h)=j\mid X(t)=i)=\delta _{ij}+q_{ij}h+o(h),} dóndeδij{\displaystyle \delta _{ij}}es la delta de Kronecker , utilizando la notación de o minúscula .qij{\displaystyle q_{ij}}puede considerarse como una medida de la rapidez con que se produce la transición de i a j .

Definición de cadena de saltos/tiempo de espera

Definimos una cadena de Markov de tiempo discreto Y n para describir el n -ésimo salto del proceso y variables S 1 , S 2 , S 3 , ... para describir los tiempos de retención en cada uno de los estados donde S i sigue la distribución exponencial con parámetro de tasa − q Y i Y i .

Definición de probabilidad de transición

Para cualquier valor n = 0, 1, 2, 3, ... y tiempos indexados hasta este valor de n : t 0 , t 1 , t 2 , ... y todos los estados registrados en estos tiempos i 0 , i 1 , i 2 , i 3 , ... se cumple que

Pr(incógnitatnorte+1=inorte+1incógnitat0=i0,incógnitat1=i1,,incógnitatnorte=inorte)=paginorteinorte+1(tnorte+1tnorte){\displaystyle \Pr(X_{t_{n+1}}=i_{n+1}\mid X_{t_{0}}=i_{0},X_{t_{1}}=i_{1},\ldots ,X_{t_{n}}=i_{n})=p_{i_{n}i_{n+1}}(t_{n+1}-t_{n})}

donde p ij es la solución de la ecuación directa (una ecuación diferencial de primer orden ).

PAG(t)=PAG(t)Q{\displaystyle P'(t)=P(t)Q}

con condición inicial P(0) es la matriz identidad .

Cadenas de Markov con interacción local

Las «cadenas de Markov con interacción local» son cadenas de Markov cuya evolución tiene en cuenta el estado de otras cadenas de Markov. Esto corresponde a la situación en la que el espacio de estados tiene forma de producto (cartesiano). Véase sistema de partículas interactuantes y autómatas celulares estocásticos (autómatas celulares probabilísticos). Véase, por ejemplo, Interacción de procesos de Markov [ 46 ] o [ 47 ] .

Proceso de Markov de tiempo discreto con espacio de estados general

Cadenas Harris

Muchos resultados para cadenas de Markov de tiempo discreto con espacio de estados finito pueden generalizarse a cadenas con espacio de estados no numerable a través de cadenas de Harris .

El uso de cadenas de Markov en los métodos de Monte Carlo de cadenas de Markov abarca casos en los que el proceso sigue un espacio de estados continuo.

Proceso de Markov en tiempo continuo con espacio de estados general

La definición de procesos de Markov en tiempo continuo con espacio de estados general es más técnica que la anterior.

Un proceso de Markov en tiempo continuoincógnita=(incógnitat)t0{\displaystyle X=(X_{t})_{t\geq 0}}es un proceso estocástico adaptado a una filtraciónF=(Ft)t0{\displaystyle \mathbb {F} =({\mathcal {F}}_{t})_{t\geq 0}}con valores en un espacio polaco localmente compacto(S,B(S)){\displaystyle (S,{\mathcal {B}}(S))}(p.ej,(R,B(R)){\displaystyle (\mathbb {R} ,{\mathcal {B}}(\mathbb {R} ))}). Esto último esencialmente garantiza que las expectativas condicionales deincógnitat{\displaystyle X_{t}}son regulares , lo que, en términos sencillos, significa que se comportan "bien". Entoncesincógnita{\displaystyle X}Se denomina proceso de Markov si satisface la propiedad de Markov , es decir, para todots0{\displaystyle t\geq s\geq 0}yAB(S){\displaystyle A\in {\mathcal {B}}(S)}[ 5 ]

PAG(incógnitatAFs)=PAG(incógnitatAincógnitas){\displaystyle P(X_{t}\in A\mid {\mathcal {F}}_{s})=P(X_{t}\in A\mid X_{s})}.

Además,incógnita{\displaystyle X}Se denomina homogéneo en el tiempo si satisface la propiedad débil de Markov para todot,s0{\displaystyle t,s\geq 0}:

PAG(incógnitat+sAFs)=PAG(incógnitatAincógnita0=incógnita)|incógnita=incógnitas=:PAGt(incógnitas,A){\displaystyle P(X_{t+s}\in A\mid {\mathcal {F}}_{s})=P(X_{t}\in A\mid X_{0}=x)|_{x=X_{s}}=:P_{t}(X_{s},A)}.

La función(t,incógnita,A)PAGt(incógnita,A){\displaystyle (t,x,A)\mapsto P_{t}(x,A)}es la llamada función de transición deincógnita{\displaystyle X}y(PAGt)t0{\displaystyle (P_{t})_{t\geq 0}}el semigrupo de transición del proceso. Las funciones de transición son generalizaciones de las matrices de transición utilizadas en el contexto de un espacio de estados finito.

De una manera más abstracta, los procesos de Markov también pueden definirse o construirse al revés: Sea(PAGt)t0{\displaystyle (P_{t})_{t\geq 0}}ser un semigrupo de transición, es decir,

  1. PAGt{\displaystyle P_{t}}es el núcleo de Markov para todost0{\displaystyle t\geq 0},
  2. PAGt+s(incógnita,A)=SPAGt(y,A)PAGs(incógnita,dy)t,s0,incógnitaR,AB(S){\displaystyle P_{t+s}(x,A)=\int _{S}P_{t}(y,A)P_{s}(x,dy)\quad \forall t,s\geq 0,x\in \mathbb {R} ,A\in {\mathcal {B}}(S)}(ecuación de Chapman-Kolmogorov),
  3. PAG0(incógnita,)=δincógnita{\displaystyle P_{0}(x,\cdot )=\delta _{x}},

dóndeδincógnita{\displaystyle \delta _{x}}es la medida de Dirac enincógnita{\displaystyle x}, yincógnita:Ω×[0,)S{\displaystyle X:\Omega \times [0,\infty )\to S}. Entoncesincógnita{\displaystyle X}es un proceso de Markov homogéneo con respecto a la filtración naturalFincógnita=(σ(incógnitas:0st))t0{\displaystyle \mathbb {F} ^{X}=(\sigma (X_{s}:0\leq s\leq t))_{t\geq 0}}, si para todos0t1<...<tnorte{\displaystyle 0\leq t_{1}<...<t_{n}},A1,...,AnorteB(S){\displaystyle A_{1},...,A_{n}\in {\mathcal {B}}(S)}la medida de probabilidad subyacentePAG{\displaystyle P}Satisface

PAG(incógnitat1A1,...,incógnitatnorteAnorteincógnita0=incógnita)=A1...Anorte1PAGtnortetnorte1(incógnitanorte1,Anorte)PAGt1(incógnita,dincógnita1){\displaystyle P(X_{t_{1}}\in A_{1},...,X_{t_{n}}\in A_{n}\mid X_{0}=x)=\int _{A_{1}}...\int _{A_{n-1}}P_{t_{n}-t_{n-1}}(x_{n-1},A_{n})\cdots P_{t_{1}}(x,dx_{1})}.

O bien, si no hay medida de probabilidadPAG{\displaystyle P}Se ha especificado, la ecuación anterior define una medidaPAGincógnita:=PAG(incógnita0=incógnita){\displaystyle P^{x}:=P(\cdot \mid X_{0}=x)}enσ(incógnitas:s0){\displaystyle \sigma (X_{s}:s\geq 0)}bajo el cual el procesoincógnita{\displaystyle X}comenzó enincógnita{\displaystyle x}es un proceso de Markov por construcción.

En otras palabras, los procesos de Markov pueden definirse como procesos estocásticos.incógnita{\displaystyle X}en un espacio de probabilidad filtrado, o indirectamente en términos de un semigrupo de transición (es decir, las probabilidades de transición del proceso), que induce un espacio de probabilidad bajo el cualincógnita{\displaystyle X}tiene la propiedad de Markov.

Propiedades

Se dice que dos estados se comunican entre sí si ambos son alcanzables el uno del otro mediante una secuencia de transiciones con probabilidad positiva. Esta es una relación de equivalencia que da lugar a un conjunto de clases comunicantes. Una clase es cerrada si la probabilidad de abandonarla es cero. Una cadena de Markov es irreducible si existe una única clase comunicante: el espacio de estados.

Un estado i tiene un período k si k es el máximo común divisor del número de transiciones por las que se puede llegar a i , partiendo de i . Es decir:

k=mcd{norte>0:Pr(incógnitanorte=iincógnita0=i)>0}{\displaystyle k=\gcd\{n>0:\Pr(X_{n}=i\mid X_{0}=i)>0\}}

El estado es periódico sik>1{\displaystyle k>1}; de lo contrariok=1{\displaystyle k=1}y el estado es aperiódico .

Se dice que un estado i es transitorio si, partiendo de i , existe una probabilidad distinta de cero de que la cadena nunca regrese a i . En caso contrario, se denomina recurrente (o persistente ). [ 48 ] Para un estado recurrente i , el tiempo medio de llegada se define como:

METROi=mi[Ti]=norte=1norteFii(norte){\displaystyle M_{i}=E[T_{i}]=\sum _{n=1}^{\infty }n\cdot f_{ii}^{(n)}}dóndeFii(norte):=Pr(min{metro>0:incógnitametro=i}=norteincógnita0=i){\displaystyle f_{ii}^{(n)}:=\Pr(\min\{m>0:X_{m}=i\}=n\mid X_{0}=i)}.

El estado i es recurrente positivo siMETROi{\displaystyle M_{i}}es finito y recurrente nulo en caso contrario. Periodicidad, transitoriedad, recurrencia y recurrencia positiva y nula son propiedades de clase; es decir, si un estado tiene la propiedad, entonces todos los estados de su clase comunicante tienen la propiedad. [ 49 ]

Un estado i se denomina absorbente si no hay transiciones salientes desde ese estado.

Irreductibilidad

Dado que la periodicidad es una propiedad de clase, si una cadena de Markov es irreducible, entonces todos sus estados tienen el mismo período. En particular, si un estado es aperiódico, entonces toda la cadena de Markov es aperiódica. [ 50 ]

Si una cadena de Markov finita es irreducible, entonces todos los estados son recurrentes positivos y tiene una distribución estacionaria única dada porπi=1/mi[Ti]{\displaystyle \pi _{i}=1/E[T_{i}]}.

Ergodicidad

Se dice que un estado i es ergódico si es aperiódico y recurrente positivo. En otras palabras, un estado i es ergódico si es recurrente, tiene un período de 1 y un tiempo medio de recurrencia finito.

Si todos los estados en una cadena de Markov irreducible son ergódicos, entonces se dice que la cadena es ergódica. Equivalentemente, existe algún enterok{\displaystyle k}de tal manera que todas las entradas deMETROk{\displaystyle M^{k}}son positivos.

Se puede demostrar que una cadena de Markov irreducible de estado finito es ergódica si tiene un estado aperiódico.

Una cadena de Markov con más de un estado y solo una transición saliente por estado no es irreducible o no es aperiódica, por lo tanto no puede ser ergódica.

Terminología

Algunos autores denominan ergódicas, incluso periódicas, a cualquier cadena de Markov recurrente positiva e irreducible. [ 51 ] De hecho, las cadenas de Markov meramente irreducibles corresponden a procesos ergódicos , definidos según la teoría ergódica . [ 52 ]

Algunos autores llaman primitiva a una matriz si existe algún número enterok{\displaystyle k}de tal manera que todas las entradas deMETROk{\displaystyle M^{k}}son positivos. [ 53 ] Algunos autores lo llaman regular . [ 54 ]

Índice de primitividad

El índice de primitividad , o exponente , de una matriz regular, es el más pequeñok{\displaystyle k}de tal manera que todas las entradas deMETROk{\displaystyle M^{k}}son positivos. El exponente es una propiedad puramente teórica de grafos, ya que depende únicamente de si cada entrada deMETRO{\displaystyle M}es cero o positivo, y por lo tanto se puede encontrar en un grafo dirigido consigramonorte(METRO){\displaystyle \mathrm {sign} (M)}como su matriz de adyacencia.

Existen varios resultados combinatorios sobre el exponente cuando hay un número finito de estados.norte{\displaystyle n}sea ​​el número de estados, entonces [ 55 ]

  • El exponente es(norte1)2+1{\displaystyle \leq (n-1)^{2}+1}. El único caso en el que es una igualdad es cuando la gráfica deMETRO{\displaystyle M}va así12norte1 y 2{\displaystyle 1\to 2\to \dots \to n\to 1{\text{ and }}2}.
  • SiMETRO{\displaystyle M}tienek1{\displaystyle k\geq 1}entradas diagonales, entonces su exponente es2nortek1{\displaystyle \leq 2n-k-1}.
  • Sisigramonorte(METRO){\displaystyle \mathrm {sign} (M)}es simétrico, entoncesMETRO2{\displaystyle M^{2}}tiene entradas diagonales positivas, lo que por la proposición anterior significa que su exponente es2norte2{\displaystyle \leq 2n-2}.
  • (Teorema de Dulmage-Mendelsohn) El exponente esnorte+s(norte2){\displaystyle \leq n+s(n-2)}dóndes{\displaystyle s}es la circunferencia del gráfico . Se puede mejorar a(d+1)+s(d+12){\displaystyle \leq (d+1)+s(d+1-2)}, dónded{\displaystyle d}es el diámetro del gráfico . [ 56 ]

Sistema dinámico que conserva la medida

Si una cadena de Markov tiene una distribución estacionaria, entonces se puede convertir en un sistema dinámico que conserva la medida : Sea el espacio de probabilidadΩ=Σnorte{\displaystyle \Omega =\Sigma ^{\mathbb {N} }}, dóndeΣ{\displaystyle \Sigma }es el conjunto de todos los estados para la cadena de Markov. Sea la sigma-álgebra en el espacio de probabilidad generada por los conjuntos de cilindros. Sea la medida de probabilidad generada por la distribución estacionaria y la transición de la cadena de Markov. SeaT:ΩΩ{\displaystyle T:\Omega \to \Omega }ser el operador de turno:T(incógnita0,incógnita1,)=(incógnita1,){\displaystyle T(X_{0},X_{1},\dots )=(X_{1},\dots )}. De manera similar podemos construir un sistema dinámico de este tipo conΩ=ΣZ{\displaystyle \Omega =\Sigma ^{\mathbb {Z} }}en cambio. [ 57 ]

Dado que las cadenas de Markov irreducibles con espacios de estados finitos tienen una distribución estacionaria única, la construcción anterior es inequívoca para las cadenas de Markov irreducibles.

En la teoría ergódica , un sistema dinámico que preserva la medida se llama ergódico si cualquier subconjunto medibleS{\displaystyle S}de tal manera queT1(S)=S{\displaystyle T^{-1}(S)=S}implicaS={\displaystyle S=\emptyset }oΩ{\displaystyle \Omega }(hasta un conjunto nulo).

La terminología es inconsistente. Dada una cadena de Markov con una distribución estacionaria estrictamente positiva en todos los estados, la cadena de Markov es irreducible si su sistema dinámico correspondiente que conserva la medida es ergódico . [ 52 ]

Representaciones markovianas

En algunos casos, los procesos aparentemente no markovianos aún pueden tener representaciones markovianas, construidas mediante la expansión del concepto de estados "actuales" y "futuros". Por ejemplo, sea X un proceso no markoviano. Entonces definamos un proceso Y , tal que cada estado de Y represente un intervalo de tiempo de estados de X. Matemáticamente, esto toma la forma:

Y(t)={incógnita(s):s[a(t),b(t)]}.{\displaystyle Y(t)={\big \{}X(s):s\in [a(t),b(t)]\,{\big \}}.}

Si Y tiene la propiedad de Markov, entonces es una representación markoviana de X.

Un ejemplo de un proceso no markoviano con una representación markoviana es una serie temporal autorregresiva de orden mayor que uno. [ 58 ]

Tiempos de golpeo

El tiempo de llegada es el tiempo que transcurre, partiendo de un conjunto de estados dado, hasta que la cadena alcanza un estado o conjunto de estados determinado. La distribución de dicho período de tiempo sigue una distribución de tipo fase. La distribución más simple corresponde a una transición con distribución exponencial.

Horarios de llegada previstos

Para un subconjunto de estados A S , el vector k A de tiempos de llegada (donde el elemento kiA{\displaystyle k_{i}^{A}}representa el valor esperado , comenzando en el estado i que la cadena entra en uno de los estados en el conjunto A ) es la solución mínima no negativa a [ 59 ]

kiA=0 para iAjSqijkjA=1 para iA.{\displaystyle {\begin{aligned}k_{i}^{A}=0&{\text{ for }}i\in A\\-\sum _{j\in S}q_{ij}k_{j}^{A}=1&{\text{ for }}i\notin A.\end{aligned}}}

inversión del tiempo

Para un proceso de Markov generalincógnita{\displaystyle X}en tiempo continuo (un CTMC o un proceso con espacio de estados general), el proceso inversoincógnita=(incógnitaTt)t[0,T]{\displaystyle {\overleftarrow {X}}=(X_{T-t})_{t\in [0,T]}}desde un tiempo fijoT>0{\displaystyle T>0}es de nuevo un proceso de Markov. Esto se deduce directamente de la propiedad de Markov : En términos informales, el futuro y el pasado son independientes dado el presente. Bajo inversión temporal, sus roles simplemente se intercambian. Sin embargo, el proceso inverso no es homogéneo en el tiempo en general. Si para algún tiempo aleatorioτ{\displaystyle \tau }(no necesariamente un momento de parada ) el proceso detenidoincógnitaτ=(incógnitatτ)t0{\displaystyle X^{\tau }=(X_{t\land \tau })_{t\geq 0}}es un proceso de Markov homogéneo en el tiempo, entonces el proceso inversoincógnitaτ=(incógnitaτtτ1{τ<})t0{\displaystyle {\overleftarrow {X^{\tau }}}=(X_{\tau -t\land \tau }1_{\{\tau <\infty \}})_{t\geq 0}}es nuevamente homogéneo en el tiempo. [ 60 ]

Siincógnita{\displaystyle X}es un CTMC, entonces por el lema de Kellyincógnita{\displaystyle {\overleftarrow {X}}}tiene la misma distribución estacionaria que el proceso directo.

Se dice que una cadena es reversible si el proceso inverso es idéntico al proceso directo (en cuanto a su distribución). El criterio de Kolmogorov establece que la condición necesaria y suficiente para que una cadena de Markov sea reversible es que el producto de las tasas de transición alrededor de un bucle cerrado sea el mismo en ambas direcciones.

Cadena de Markov embebida

Un método para encontrar la distribución de probabilidad estacionaria , π , de una cadena de Markov ergódica de tiempo continuo, Q , consiste en encontrar primero su cadena de Markov embebida (CME) . Estrictamente hablando, la CME es una cadena de Markov regular de tiempo discreto, a veces denominada proceso de salto . Cada elemento de la matriz de probabilidad de transición de un paso de la CME, S , se denota por s ij , y representa la probabilidad condicional de transición del estado i al estado j . Estas probabilidades condicionales pueden encontrarse mediante

sij={qijkiqiksi ij0de lo contrario.{\displaystyle s_{ij}={\begin{cases}{\frac {q_{ij}}{\sum _{k\neq i}q_{ik}}}&{\text{if }}i\neq j\\0&{\text{otherwise}}.\end{cases}}}

A partir de esto, S puede escribirse como

S=I(diagnóstico(Q))1Q{\displaystyle S=I-\left(\operatorname {diag} (Q)\right)^{-1}Q}

donde I es la matriz identidad y diag( Q ) es la matriz diagonal formada al seleccionar la diagonal principal de la matriz Q y establecer todos los demás elementos a cero.

Para encontrar el vector de distribución de probabilidad estacionaria, debemos encontrar a continuaciónφ{\displaystyle \varphi }de tal manera que

φS=φ,{\displaystyle \varphi S=\varphi ,}

conφ{\displaystyle \varphi }siendo un vector fila, de tal manera que todos los elementos enφ{\displaystyle \varphi }son mayores que 0 yφ1{\displaystyle \|\varphi \|_{1}}= 1. A partir de esto, se puede hallar π como

π=φ(diagnóstico(Q))1φ(diagnóstico(Q))11.{\displaystyle \pi ={-\varphi (\operatorname {diag} (Q))^{-1} \over \left\|\varphi (\operatorname {diag} (Q))^{-1}\right\|_{1}}.}

( S puede ser periódico, incluso si Q no lo es. Una vez que se encuentra π , debe normalizarse a un vector unitario ).

Otro proceso de tiempo discreto que puede derivarse de una cadena de Markov de tiempo continuo es un esqueleto δ : la cadena de Markov (de tiempo discreto) formada al observar X ( t ) a intervalos de δ unidades de tiempo. Las variables aleatorias X (0), X (δ), X (2δ), ... dan la secuencia de estados visitados por el esqueleto δ.   

Tipos especiales de cadenas de Markov

modelo de Markov

Los modelos de Markov se utilizan para modelar sistemas cambiantes. Existen cuatro tipos principales de modelos que generalizan las cadenas de Markov dependiendo de si cada estado secuencial es observable o no, y de si el sistema debe ajustarse en función de las observaciones realizadas:

Esquema de Bernoulli

Un esquema de Bernoulli es un caso especial de cadena de Markov donde la matriz de probabilidad de transición tiene filas idénticas, lo que significa que el siguiente estado es independiente incluso del estado actual (además de ser independiente de los estados anteriores). Un esquema de Bernoulli con solo dos estados posibles se conoce como proceso de Bernoulli .

Sin embargo, según el teorema de isomorfismo de Ornstein , toda cadena de Markov aperiódica e irreducible es isomorfa a un esquema de Bernoulli; [ 61 ] por lo tanto, también se podría afirmar que las cadenas de Markov son un "caso especial" de los esquemas de Bernoulli. El isomorfismo generalmente requiere una recodificación compleja. El teorema de isomorfismo es incluso un poco más fuerte: afirma que cualquier proceso estocástico estacionario es isomorfo a un esquema de Bernoulli; la cadena de Markov es solo un ejemplo de ello.

Subdesplazamiento de tipo finito

Cuando la matriz de Markov se reemplaza por la matriz de adyacencia de un grafo finito , el desplazamiento resultante se denomina cadena de Markov topológica o subdesplazamiento de tipo finito . [ 61 ] Una matriz de Markov compatible con la matriz de adyacencia puede entonces proporcionar una medida sobre el subdesplazamiento. Muchos sistemas dinámicos caóticos son isomorfos a cadenas de Markov topológicas; ejemplos de ello son los difeomorfismos de variedades cerradas , el sistema de Prouhet-Thue-Morse , el sistema de Chacon , los sistemas sóficos , los sistemas libres de contexto y los sistemas de codificación por bloques . [ 61 ]

Aplicaciones

Las cadenas de Markov se han utilizado en una amplia gama de temas en las ciencias naturales y sociales, así como en aplicaciones tecnológicas.

Física

Los sistemas markovianos aparecen ampliamente en termodinámica y mecánica estadística , siempre que se utilicen probabilidades para representar detalles desconocidos o no modelados del sistema, si se puede asumir que la dinámica es invariante en el tiempo y que no es necesario considerar ningún historial relevante que no esté ya incluido en la descripción del estado. [ 62 ] [ 63 ] Por ejemplo, un estado termodinámico opera bajo una distribución de probabilidad que es difícil o costosa de obtener. Por lo tanto, el método de Monte Carlo de cadena de Markov se puede utilizar para extraer muestras aleatoriamente de una caja negra para aproximar la distribución de probabilidad de los atributos en un rango de objetos. [ 63 ]

Las cadenas de Markov se utilizan en simulaciones de QCD en la red . [ 64 ]

Química

mi+SmiSustratovinculanteSmiCatalíticopaso+PAG{\displaystyle {\ce {{E}+{\underset {Substrate \atop binding}{S<=>E}}{\overset {Catalytic \atop step}{S->E}}+P}}}
Cinética de Michaelis-Menten . La enzima (E) se une a un sustrato (S) y produce un producto (P). Cada reacción es una transición de estado en una cadena de Markov.

Una red de reacciones es un sistema químico que involucra múltiples reacciones y especies químicas . Los modelos estocásticos más simples de tales redes tratan el sistema como una cadena de Markov de tiempo continuo donde el estado es el número de moléculas de cada especie y con las reacciones modeladas como posibles transiciones de la cadena. [ 65 ] Las cadenas de Markov y los procesos de Markov de tiempo continuo son útiles en química cuando los sistemas físicos se aproximan mucho a la propiedad de Markov. Por ejemplo, imaginemos un gran número n de moléculas en solución en el estado A, cada una de las cuales puede experimentar una reacción química al estado B con una cierta velocidad promedio. Tal vez la molécula sea una enzima , y ​​los estados se refieren a cómo está plegada . El estado de cualquier enzima individual sigue una cadena de Markov, y dado que las moléculas son esencialmente independientes entre sí, el número de moléculas en el estado A o B en un momento dado es n veces la probabilidad de que una molécula dada esté en ese estado.

El modelo clásico de actividad enzimática, la cinética de Michaelis-Menten , puede verse como una cadena de Markov, donde en cada paso de tiempo la reacción procede en alguna dirección. Si bien la cinética de Michaelis-Menten es bastante sencilla, también se pueden modelar redes de reacción mucho más complejas con cadenas de Markov. [ 66 ]

También se utilizó un algoritmo basado en una cadena de Markov para enfocar el crecimiento de compuestos químicos mediante la fragmentación in silico hacia una clase deseada de compuestos, como fármacos o productos naturales. [ 67 ] A medida que crece una molécula, se selecciona un fragmento de la molécula naciente como el estado "actual". Este no conoce su pasado (es decir, no sabe qué está ya unido a él). Luego, pasa al siguiente estado cuando se le une un fragmento. Las probabilidades de transición se entrenan con bases de datos de clases auténticas de compuestos. [ 68 ]

Asimismo, el crecimiento (y la composición) de los copolímeros puede modelarse mediante cadenas de Markov. A partir de las relaciones de reactividad de los monómeros que componen la cadena polimérica en crecimiento, se puede calcular la composición de la cadena (por ejemplo, si los monómeros tienden a añadirse de forma alternada o en largas secuencias del mismo monómero). Debido a efectos estéricos , los efectos de Markov de segundo orden también pueden influir en el crecimiento de algunas cadenas poliméricas.

De manera similar, se ha sugerido que la cristalización y el crecimiento de algunos materiales de óxido de superred epitaxial pueden describirse con precisión mediante cadenas de Markov. [ 69 ]

Biología

Las cadenas de Markov se utilizan en diversas áreas de la biología. Algunos ejemplos notables son:

teoría de la información

Las cadenas de Markov se utilizan en todo el procesamiento de la información. El famoso artículo de Claude Shannon de 1948 , "Una teoría matemática de la comunicación" , que en un solo paso creó el campo de la teoría de la información , comienza introduciendo el concepto de entropía al modelar textos en un lenguaje natural (como el inglés) generados por un proceso de Markov ergódico, donde cada letra puede depender estadísticamente de las letras anteriores. [ 72 ] Estos modelos idealizados pueden capturar muchas de las regularidades estadísticas de los sistemas. Incluso sin describir la estructura completa del sistema a la perfección, estos modelos de señales pueden hacer posible una compresión de datos muy efectiva a través de técnicas de codificación de entropía como la codificación aritmética . También permiten una estimación de estado y un reconocimiento de patrones efectivos . Las cadenas de Markov también juegan un papel importante en el aprendizaje por refuerzo .

Las cadenas de Markov también son la base de los modelos ocultos de Markov, que son una herramienta importante en campos tan diversos como las redes telefónicas (que utilizan el algoritmo de Viterbi para la corrección de errores), el reconocimiento de voz y la bioinformática (como en la detección de reordenamientos [ 73 ] ).

El algoritmo de compresión de datos sin pérdidas LZMA combina cadenas de Markov con la compresión Lempel-Ziv para lograr índices de compresión muy altos.

teoría de colas

Las cadenas de Markov son la base del tratamiento analítico de las colas ( teoría de colas ). Agner Krarup Erlang inició el tema en 1917. [ 74 ] Esto las hace fundamentales para optimizar el rendimiento de las redes de telecomunicaciones, donde los mensajes a menudo deben competir por recursos limitados (como el ancho de banda). [ 75 ]

Numerosos modelos de colas utilizan cadenas de Markov de tiempo continuo. Por ejemplo, una cola M/M/1 es una CTMC sobre los enteros no negativos, donde las transiciones ascendentes de i a i  +  1 ocurren a una tasa λ según un proceso de Poisson y describen las llegadas de trabajos, mientras que las transiciones de i a i  1 (para i  >  1) ocurren a una tasa μ (los tiempos de servicio de los trabajos se distribuyen exponencialmente) y describen los servicios completados (salidas) de la cola.

aplicaciones de Internet

Un diagrama de estados que representa el algoritmo PageRank con una probabilidad de transición de M, oαki+1αnorte{\displaystyle {\frac {\alpha }{k_{i}}}+{\frac {1-\alpha }{N}}}

El PageRank de una página web, tal como lo utiliza Google , se define mediante una cadena de Markov. [ 76 ] [ 77 ] [ 78 ] Es la probabilidad de estar en la páginai{\displaystyle i}en la distribución estacionaria en la siguiente cadena de Markov en todas las páginas web (conocidas). Sinorte{\displaystyle N}es el número de páginas web conocidas, y una páginai{\displaystyle i}tieneki{\displaystyle k_{i}}enlaces salientes desde él entonces tiene probabilidad de transiciónαki+1αnorte{\displaystyle {\frac {\alpha }{k_{i}}}+{\frac {1-\alpha }{N}}}para todas las páginas que están vinculadas y1αnorte{\displaystyle {\frac {1-\alpha }{N}}}para todas las páginas que no están enlazadas. El parámetroα{\displaystyle \alpha }se considera que es aproximadamente 0,85. [ 79 ]

Los modelos de Markov también se han utilizado para analizar el comportamiento de navegación web de los usuarios. La transición de un usuario a través de los enlaces de un sitio web específico puede modelarse mediante modelos de Markov de primer o segundo orden, y pueden utilizarse para realizar predicciones sobre la navegación futura y para personalizar la página web para cada usuario.

Estadística

Los métodos de cadena de Markov también se han vuelto muy importantes para generar secuencias de números aleatorios que reflejen con precisión distribuciones de probabilidad deseadas muy complejas, mediante un proceso llamado Monte Carlo de cadena de Markov (MCMC). En los últimos años, esto ha revolucionado la viabilidad de los métodos de inferencia bayesiana , permitiendo simular una amplia gama de distribuciones posteriores y hallar numéricamente sus parámetros.

Economía y finanzas

Las cadenas de Markov se utilizan en finanzas y economía para modelar una variedad de fenómenos diferentes, incluyendo la distribución del ingreso, la distribución del tamaño de las empresas, los precios de los activos y las caídas del mercado. DG Champernowne construyó un modelo de cadena de Markov de la distribución del ingreso en 1953. [ 80 ] Herbert A. Simon y su coautor Charles Bonini utilizaron un modelo de cadena de Markov para derivar una distribución estacionaria de Yule de tamaños de empresas. [ 81 ] Louis Bachelier fue el primero en observar que los precios de las acciones seguían un paseo aleatorio. [ 82 ] El paseo aleatorio fue visto más tarde como evidencia a favor de la hipótesis del mercado eficiente y los modelos de paseo aleatorio fueron populares en la literatura de la década de 1960. [ 83 ] Los modelos de cambio de régimen de los ciclos económicos fueron popularizados por James D. Hamilton (1989), quien utilizó una cadena de Markov para modelar cambios entre períodos de alto y bajo crecimiento del PIB (o, alternativamente, expansiones y recesiones económicas). [ 84 ] Un ejemplo más reciente es el modelo multifractal de conmutación de Markov de Laurent E. Calvet y Adlai J. Fisher, que se basa en la conveniencia de modelos anteriores de cambio de régimen. [ 85 ] [ 86 ] Utiliza una cadena de Markov arbitrariamente grande para impulsar el nivel de volatilidad de los rendimientos de los activos.

La macroeconomía dinámica hace un uso extensivo de las cadenas de Markov. Un ejemplo es el uso de cadenas de Markov para modelar exógenamente los precios de las acciones en un entorno de equilibrio general . [ 87 ]

Las agencias de calificación crediticia elaboran tablas anuales de las probabilidades de transición para bonos de diferentes calificaciones crediticias. [ 88 ]

ciencias sociales

Las cadenas de Markov se utilizan generalmente para describir argumentos dependientes de la trayectoria , donde las configuraciones estructurales actuales condicionan los resultados futuros. Un ejemplo es la reformulación de la idea, originalmente debida a El Capital de Karl Marx , que vincula el desarrollo económico con el auge del capitalismo . En la investigación actual, es común utilizar una cadena de Markov para modelar cómo, una vez que un país alcanza un nivel específico de desarrollo económico, la configuración de factores estructurales, como el tamaño de la clase media , la proporción de residencia urbana frente a rural, la tasa de movilización política , etc., generará una mayor probabilidad de transición de un régimen autoritario a uno democrático . [ 89 ]

Música

Las cadenas de Markov se emplean en la composición musical algorítmica , particularmente en software como Csound , Max y SuperCollider . En una cadena de primer orden, los estados del sistema se convierten en valores de notas o tonos, y se construye un vector de probabilidad para cada nota, completando una matriz de probabilidad de transición (véase más abajo). Se construye un algoritmo para producir valores de notas de salida basados ​​en las ponderaciones de la matriz de transición, que pueden ser valores de notas MIDI , frecuencia ( Hz ) o cualquier otra métrica deseada. [ 90 ]

Una cadena de Markov de segundo orden puede introducirse considerando el estado actual y también el estado anterior, como se indica en la segunda tabla. Las cadenas de orden superior, de orden n, tienden a "agrupar" notas particulares, mientras que ocasionalmente se "desvían" hacia otros patrones y secuencias. Estas cadenas de orden superior tienden a generar resultados con una estructura fraseológica , en lugar del "vagabundeo sin rumbo" producido por un sistema de primer orden. [ 91 ]

Las cadenas de Markov pueden utilizarse estructuralmente, como en Analogique A y B de Xenakis. [ 92 ] Las cadenas de Markov también se utilizan en sistemas que emplean un modelo de Markov para reaccionar interactivamente a la entrada musical. [ 93 ]

Por lo general, los sistemas musicales necesitan imponer restricciones de control específicas a las secuencias de longitud finita que generan, pero estas restricciones no son compatibles con los modelos de Markov, ya que inducen dependencias de largo alcance que violan la hipótesis de memoria limitada de Markov. Para superar esta limitación, se ha propuesto un nuevo enfoque. [ 94 ]

Juegos y deportes

Las cadenas de Markov pueden utilizarse para modelar muchos juegos de azar. Los juegos infantiles Serpientes y Escaleras y " ¡Hola, Cherry-O! ", por ejemplo, se representan con precisión mediante cadenas de Markov. En cada turno, el jugador comienza en un estado determinado (en una casilla específica) y, a partir de ahí, tiene una probabilidad fija de pasar a otros estados (casillas).

Los modelos de cadena de Markov se han utilizado en el análisis avanzado del béisbol desde 1960, aunque su uso sigue siendo poco frecuente. Cada media entrada de un partido de béisbol se ajusta al estado de la cadena de Markov cuando se considera el número de corredores y outs. Durante cualquier turno al bate, hay 24 combinaciones posibles de número de outs y posición de los corredores. Mark Pankin demuestra que los modelos de cadena de Markov pueden utilizarse para evaluar las carreras generadas tanto por jugadores individuales como por un equipo. [ 95 ] También analiza diversos tipos de estrategias y condiciones de juego: cómo se han utilizado los modelos de cadena de Markov para analizar estadísticas en situaciones de juego como el toque de bola y el robo de bases , y las diferencias al jugar en césped natural frente a césped artificial . [ 96 ]

generadores de texto de Markov

Los procesos de Markov también pueden utilizarse para generar texto de apariencia superficialmente real a partir de un documento de muestra. Estos procesos se emplean en diversos programas de generación de parodias con fines recreativos (véase Dissociated Press , Jeff Harrison, [ 97 ] Mark V. Shaney , [ 98 ] [ 99 ] y Academias Neutronium). Existen varias bibliotecas de código abierto para la generación de texto mediante cadenas de Markov.

Véase también

Notas

  1. 1 2 Sean Meyn; Richard L. Tweedie (2 de abril de 2009). Cadenas de Markov y estabilidad estocástica . Cambridge University Press. pág.  3. ISBN 978-0-521-73182-9.
  2. Reuven Y. Rubinstein; Dirk P. Kroese (20 de septiembre de 2011). Simulación y Método Montecarlo . John Wiley e hijos. pag. 225.ISBN  978-1-118-21052-9.
  3. Dani Gamerman; Hedibert F. Lopes (10 de mayo de 2006). Markov Chain Monte Carlo: Stochastic Simulation for Bayesian Inference, Second Edition . CRC Press. ISBN 978-1-58488-587-0.
  4. "Markoviano" . Oxford English Dictionary ( edición en línea). Oxford University Press. (Se requiere suscripción o ser miembro de una institución participante ).
  5. 1 2 Øksendal, BK (Bernt Karsten) (2003). Ecuaciones diferenciales estocásticas: una introducción con aplicaciones (6.ª ed.). Berlín: Springer. ISBN  3-540-04758-1OCLC 52203046 
  6. ^ Søren Asmussen (15 de mayo de 2003) . Probabilidad Aplicada y Colas . Medios de ciencia y negocios de Springer. pag. 7.ISBN  978-0-387-00211-8.
  7. Emanuel Parzen (17 de junio de 2015). Procesos estocásticos . Courier Dover Publications. pág. 188. ISBN  978-0-486-79688-8.
  8. Samuel Karlin; Howard E. Taylor (2 de diciembre de 2012). Un primer curso sobre procesos estocásticos . Academic Press. págs. 29 y 30. ISBN  978-0-08-057041-9.
  9. John Lamperti (1977). Procesos estocásticos: una revisión de la teoría matemática . Springer-Verlag. págs. 106–121 . ISBN  978-3-540-90275-1.
  10. Sheldon M. Ross (1996). Procesos estocásticos . Wiley. págs. 174 y 231. ISBN  978-0-471-12062-9.
  11. Everitt, BS (2002) The Cambridge Dictionary of Statistics . CUP. ISBN 0-521-81099-X
  12. Parzen, E. (1962) Procesos estocásticos , Holden-Day. ISBN 0-8162-6664-6(Tabla 6.1)
  13. Dodge, Y. (2003) The Oxford Dictionary of Statistical Terms , OUP. ISBN 0-19-920613-9(entrada para "cadena de Markov")
  14. Dodge, Y. The Oxford Dictionary of Statistical Terms , OUP. ISBN 0-19-920613-9
  15. Meyn, S., Sean P. y Richard L. Tweedie. (2009) Cadenas de Markov y estabilidad estocástica . Cambridge University Press. (Prefacio, pág. iii)
  16. 1 2 3 4 5 Charles Miller Grinstead; James Laurie Snell (1997). Introducción a la probabilidad . American Mathematical Soc. págs. 464–466 . ISBN  978-0-8218-0749-1.
  17. 1 2 3 Pierre Bremaud (9 de marzo de 2013). Cadenas de Markov: campos de Gibbs, simulación de Monte Carlo y colas . Springer Science & Business Media. pág. ix. ISBN  978-1-4757-3124-8.
  18. 1 2 3 Hayes, Brian (2013). "Primeros eslabones en la cadena de Markov". American Scientist . 101 (2): 92– 96. doi : 10.1511/2013.101.92 .
  19. 1 2 Sheldon M. Ross (1996). Procesos estocásticos . Wiley. págs. 235 y 358. ISBN  978-0-471-12062-9.
  20. Jarrow, Robert; Protter, Philip (2004). «Una breve historia de la integración estocástica y las finanzas matemáticas: Los primeros años, 1880-1970». Un volumen conmemorativo para Herman Rubin . págs. 75-91 . CiteSeerX 10.1.1.114.632 . doi : 10.1214/lnms/1196285381 . ISBN   978-0-940600-61-4.
  21. Guttorp, Peter; Thorarinsdottir, Thordis L. (2012). "¿Qué pasó con el caos discreto, el proceso de Quenouille y la propiedad de Markov aguda? Algo de historia de los procesos puntuales estocásticos". International Statistical Review . 80 (2): 253– 268. doi : 10.1111/j.1751-5823.2012.00181.x .
  22. Seneta, E. (1996). "Markov y el nacimiento de la teoría de la dependencia en cadena". International Statistical Review . 64 (3): 255– 257. doi : 10.2307/1403785 . JSTOR 1403785 . 
  23. Seneta, E. (1998). "IJ Bienaymé [1796–1878]: Criticidad, desigualdad e internacionalización". International Statistical Review . 66 (3): 291– 292. doi : 10.2307/1403518 . JSTOR 1403518 . 
  24. ^ Bru B, Hertz S (2001). "Maurice Fréchet". En Heyde CC , Seneta E, Crépel P, Fienberg SE, Gani J (eds.). Estadísticos de los siglos . Nueva York, Nueva York: Springer. págs. 331– 334. doi : 10.1007/978-1-4613-0179-0_71 . ISBN  978-0-387-95283-3.
  25. 1 2 3 Kendall, DG; Batchelor, GK; Bingham, NH; Hayman, WK; Hyland, JME; Lorentz, GG; Moffatt, HK; Parry, W.; Razborov, AA; Robinson, CA; Whittle, P. (1990). "Andrei Nikolaevich Kolmogorov (1903–1987)". Boletín de la Sociedad Matemática de Londres . 22 (1): 33. doi : 10.1112/blms/22.1.31 .
  26. 1 2 Cramér, Harald (1976). "Medio siglo con la teoría de la probabilidad: algunos recuerdos personales" . The Annals of Probability . 4 (4): 509– 546. doi : 10.1214/aop/1176996025 .
  27. Marc Barbut; Bernard Locker; Laurent Mazliak (23 de agosto de 2016). Paul Lévy y Maurice Fréchet: 50 años de correspondencia en 107 cartas . Springer Londres. pag. 5.ISBN  978-1-4471-7262-8.
  28. Valeriy Skorokhod (5 de diciembre de 2005). Principios básicos y aplicaciones de la teoría de la probabilidad . Springer Science & Business Media. pág. 146. ISBN  978-3-540-26312-8.
  29. Bernstein, Jeremy (2005). "Bachelier". American Journal of Physics . 73 (5): 395– 398. Bibcode : 2005AmJPh..73..395B . doi : 10.1119/1.1848117 .
  30. William J. Anderson (6 de diciembre de 2012). Cadenas de Markov de tiempo continuo: un enfoque orientado a las aplicaciones . Springer Science & Business Media. pág. vii. ISBN  978-1-4612-3038-0.
  31. Kendall, DG; Batchelor, GK; Bingham, NH; Hayman, WK; Hyland, JME; Lorentz, GG; Moffatt, HK; Parry, W.; Razborov, AA; Robinson, CA; Whittle, P. (1990). "Andrei Nikolaevich Kolmogorov (1903–1987)". Boletín de la Sociedad Matemática de Londres . 22 (1): 57. doi : 10.1112/blms/22.1.31 .
  32. Subramanian, Devika (otoño de 2008). "El curioso caso de Mark V. Shaney" (PDF) . Ciencias de la Computación. Apuntes del curso Comp 140, otoño de 2008. Universidad William Marsh Rice . Consultado el 30 de noviembre de 2024 .
  33. 1 2 Ionut Florescu (7 de noviembre de 2014). Probabilidad y procesos estocásticos . John Wiley & Sons. págs. 373 y 374. ISBN  978-1-118-59320-2.
  34. 1 2 Samuel Karlin; Howard E. Taylor (2 de diciembre de 2012). Un primer curso sobre procesos estocásticos . Academic Press. pág. 49. ISBN  978-0-08-057041-9.
  35. Weiss, George H. (2006). "Random Walks". Encyclopedia of Statistical Sciences . p. 1. doi : 10.1002/0471667196.ess2180.pub2 . ISBN  978-0-471-66719-3.
  36. Michael F. Shlesinger (1985). El maravilloso mundo de la estocástica: un homenaje a Elliott W. Montroll . North-Holland. págs. 8–10 . ISBN  978-0-444-86937-1.
  37. Emanuel Parzen (17 de junio de 2015). Procesos estocásticos . Courier Dover Publications. págs. 7, 8. ISBN  978-0-486-79688-8.
  38. Joseph L. Doob (1990). Procesos estocásticos . Wiley. págs. 46, 47. 
  39. Donald L. Snyder; Michael I. Miller (6 de diciembre de 2012). Procesos puntuales aleatorios en el tiempo y el espacio . Springer Science & Business Media. pág. 32. ISBN  978-1-4612-3166-0.
  40. 1 2 Serfozo, Richard (2009). Fundamentos de los procesos estocásticos aplicados . Probabilidad y sus aplicaciones. Berlín: Springer. doi : 10.1007/978-3-540-89332-5 . ISBN 978-3-540-89331-8.
  41. "Capítulo 11 "Cadenas de Markov"" (PDF) . Archivado del original (PDF) el 15-02-2017 . Recuperado el 02-06-2017 .
  42. Schmitt, Florian; Rothlauf, Franz (2001). "Sobre la importancia del segundo autovalor más grande en la tasa de convergencia de los algoritmos genéticos". Actas del 14.º Simposio sobre Sistemas Distribuidos Confiables . CiteSeerX 10.1.1.28.6191 . 
  43. Rosenthal, Jeffrey S. (1995). "Tasas de convergencia para cadenas de Markov". SIAM Review . 37 (3): 387– 405. doi : 10.1137/1037083 . JSTOR 2132659 . 
  44. Franzke, Brandon; Kosko, Bart (1 de octubre de 2011). "El ruido puede acelerar la convergencia en cadenas de Markov". Physical Review E. 84 ( 4) 041112. Bibcode : 2011PhRvE..84d1112F . doi : 10.1103/PhysRevE.84.041112 . PMID 22181092 . 
  45. Norris, JR (1997). "Cadenas de Markov de tiempo continuo I". Cadenas de Markov . págs. 60–107 . doi : 10.1017/CBO9780511810633.004 . ISBN  978-0-511-81063-3.
  46. Spitzer, Frank (1970). "Interacción de procesos de Markov" . Advances in Mathematics . 5 (2): 246– 290. Bibcode : 1970AdMat...5..246S . doi : 10.1016/0001-8708(70)90034-4 .
  47. Dobrushin, RL ; Kryukov, VI; Toom, AL (1978). Sistemas celulares estocásticos: ergodicidad, memoria, morfogénesis . Manchester University Press. ISBN 978-0-7190-2206-7. Consultado el 4 de marzo de 2016 .
  48. Heyman, Daniel P.; Sobel, Mathew J. (1982). Modelos estocásticos en investigación operativa, volumen 1. Nueva York: McGraw-Hill. pág. 230. ISBN  0-07-028631-0.
  49. Peres, Yuval . "Demuestra que la recurrencia positiva es una propiedad de clase" . Mathematics Stack Exchange . Consultado el 1 de febrero de 2024 .
  50. Lalley, Steve (2016). "Cadenas de Markov: Teoría básica" (PDF) . Consultado el 22 de junio de 2024 .
  51. Parzen, Emanuel (1962). Procesos estocásticos . San Francisco: Holden-Day. pág. 145. ISBN  0-8162-6664-6.{{cite book}}: Incompatibilidad de ISBN/Fecha ( ayuda )
  52. 1 2 Shalizi, Cosma (1 de diciembre de 2023). "Teoría ergódica" . bactra.org . Recuperado el 1 de febrero de 2024 .
  53. Seneta, E. (Eugene) (1973). Matrices no negativas: una introducción a la teoría y las aplicaciones . Archivo de Internet. Nueva York, Wiley. ISBN 978-0-470-77605-6.
  54. "10.3: Cadenas de Markov regulares" . Matemáticas LibreTexts . 22 de marzo de 2020. Consultado el 1 de febrero de 2024 .
  55. Seneta, E. (Eugene) (1973). "2.4. Propiedades combinatorias". Matrices no negativas: una introducción a la teoría y las aplicaciones . Archivo de Internet. Nueva York, Wiley. ISBN 978-0-470-77605-6.
  56. Shen, Jian (1996-10-15). "Una mejora del teorema de Dulmage-Mendelsohn" . Matemáticas Discretas . 158 (1): 295– 297. doi : 10.1016/0012-365X(95)00060-A .
  57. Kallenberg, Olav (2002). Fundamentos de la probabilidad moderna . Probabilidad y sus aplicaciones (2.ª ed., [Nachdr.] ed.). Nueva York, NY Berlín Heidelberg: Springer. Proposición 8.6 (página 145). ISBN  978-0-387-95313-7.
  58. Doblinger, G. (septiembre de 1998). "Suavizado de señales AR ruidosas mediante un filtro de Kalman adaptativo" (PDF) . 9.ª Conferencia Europea de Procesamiento de Señales (EUSIPCO 1998) : 781–784 .
  59. Norris, JR (1997). "Cadenas de Markov de tiempo continuo II". Cadenas de Markov . págs. 108–127 . doi : 10.1017/CBO9780511810633.005 . ISBN  978-0-511-81063-3.
  60. Chung, Kai Lai; Walsh, John B. (2006). Procesos de Markov, movimiento browniano y simetría temporal (2.ª ed.). Springer Nueva York. pág. 304. ISBN   978-0-387-28696-9.
  61. 1 2 3 Matthew Nicol y Karl Petersen, (2009) " Teoría ergódica: ejemplos básicos y construcciones ", Enciclopedia de la complejidad y la ciencia de sistemas , Springer https://doi.org/10.1007/978-0-387-30440-3_177
  62. Fitzpatrick, Richard. "Termodinámica y mecánica estadística" (PDF) . Archivado del original (PDF) el 30-11-2016 . Consultado el 02-06-2017 .
  63. 1 2 van Ravenzwaaij, Don; Cassey, Pete; Brown, Scott D. (2016-03-11). "Una introducción sencilla al muestreo de Monte Carlo de cadenas de Markov" . Psychonomic Bulletin & Review . 25 (1): 143– 154. doi : 10.3758/s13423-016-1015-8 . PMC 5862921. PMID 26968853 .  
  64. Gattringer, Christof; Lang, Christian B (2010). Cromodinámica cuántica en la red . Lecture Notes in Physics. Vol. 788. Springer-Verlag Berlin Heidelberg. doi : 10.1007/978-3-642-01850-3 . ISBN  978-3-642-01849-7.
  65. Anderson, David F.; Kurtz, Thomas G. (2011), "Continuous Time Markov Chain Models for Chemical Reaction Networks", Design and Analysis of Biomolecular Circuits , Springer New York, pp. 3–42 , doi : 10.1007/978-1-4419-6766-4_1 , ISBN  978-1-4419-6765-7{{citation}}: CS1 mantenimiento: parámetro de trabajo con ISBN ( enlace )
  66. Du, Chao; Kou, SC (septiembre de 2012). "Análisis de correlación de la reacción enzimática de una sola molécula de proteína" . The Annals of Applied Statistics . 6 (3): 950– 976. arXiv : 1209.6210 . Bibcode : 2012arXiv1209.6210D . doi : 10.1214/12-aoas541 . PMC 3568780. PMID 23408514 .  
  67. Kutchukian, Peter; Lou, David; Shakhnovich, Eugene (2009). "FOG: Algoritmo de crecimiento optimizado de fragmentos para la generación de novo de moléculas que ocupan sitios químicos similares a fármacos". Journal of Chemical Information and Modeling . 49 (7): 1630– 1642. doi : 10.1021/ci9000458 . PMID 19527020 . 
  68. Kutchukian, PS; Lou, D.; Shakhnovich, Eugene I. (2009-06-15). "FOG: Algoritmo de crecimiento optimizado de fragmentos para la generación de novo de moléculas que ocupan un espacio químico similar al de los fármacos". Journal of Chemical Information and Modeling . 49 (7): 1630– 1642. doi : 10.1021/ci9000458 . PMID 19527020 . 
  69. Kopp, VS; Kaganer, VM; Schwarzkopf, J.; Waidick, F.; Remmele, T.; Kwasniewski, A.; Schmidbauer, M. (2011). "Difracción de rayos X de estructuras en capas no periódicas con correlaciones: cálculo analítico y experimento en películas mixtas de Aurivillius". Acta Crystallographica Sección A . 68 (Pt 1): 148– 155. Bibcode : 2012AcCrA..68..148K . doi : 10.1107/S0108767311044874 . PMID 22186291 . 
  70. George, Dileep; Hawkins, Jeff (2009). Friston, Karl J. (ed.). "Hacia una teoría matemática de los microcircuitos corticales" . PLOS Comput Biol . 5 (10) e1000532. Bibcode : 2009PLSCB...5E0532G . doi : 10.1371/journal.pcbi.1000532 . PMC 2749218. PMID 19816557 .  
  71. Gupta, Ankur; Rawlings, James B. (abril de 2014). "Comparación de métodos de estimación de parámetros en modelos cinéticos químicos estocásticos: ejemplos en biología de sistemas" . AIChE Journal . 60 (4): 1253– 1268. Bibcode : 2014AIChE..60.1253G . doi : 10.1002/aic.14409 . PMC 4946376. PMID 27429455 .  
  72. Thomsen, Samuel W. (2009), "Algunas evidencias sobre la génesis de la teoría de la información de Shannon", Studies in History and Philosophy of Science , 40 (1): 81– 91, Bibcode : 2009SHPSA..40...81T , doi : 10.1016/j.shpsa.2008.12.011
  73. Pratas, D; Silva, R; Pinho, A; Ferreira, P (18 de mayo de 2015). "Un método sin alineación para encontrar y visualizar reordenamientos entre pares de secuencias de ADN" . Informes científicos . 5 (10203) 10203. Código Bib : 2015NatSR...510203P . doi : 10.1038/srep10203 . PMC 4434998 . PMID 25984837 .  
  74. O'Connor, John J.; Robertson, Edmund F. , "Cadena de Markov" , Archivo MacTutor de Historia de las Matemáticas , Universidad de St Andrews
  75. SP Meyn, 2007. Control Techniques for Complex Networks Archived 2015-05-13 at the Wayback Machine , Cambridge University Press, 2007.
  76. Patente estadounidense 6,285,999
  77. Gupta, Brij; Agrawal, Dharma P.; Yamaguchi, Shingo (16 de mayo de 2016). Manual de investigación sobre soluciones criptográficas modernas para la seguridad informática y cibernética . IGI Global. págs. 448–. ISBN  978-1-5225-0106-0.
  78. Langville, Amy N.; Meyer, Carl D. (2006). "Una reordenación para el problema PageRank" (PDF) . SIAM Journal on Scientific Computing . 27 (6): 2112– 2113. Bibcode : 2006SJSC...27.2112L . CiteSeerX 10.1.1.58.8652 . doi : 10.1137/040607551 . Archivado del original (PDF) el 21-09-2017 . Recuperado el 07-11-2017 . 
  79. Page, Lawrence; Brin, Sergey; Motwani, Rajeev; Winograd, Terry (1999). El PageRank Citation Ranking: Trayendo orden a la web (Informe técnico). CiteSeerX 10.1.1.31.1768 . 
  80. Champernowne, D (1953). "Un modelo de distribución del ingreso". The Economic Journal . 63 (250): 318– 51. doi : 10.2307/2227127 . JSTOR 2227127 . 
  81. Simon, Herbert; C Bonini (1958). "La distribución del tamaño de las empresas comerciales". Am. Econ. Rev. 42 : 425–40 .
  82. ^ Bachelier, Luis (1900). "Teoría de la especulación". Annales Scientifiques de l'École Normale Supérieure . 3 : 21– 86. doi : 10.24033/asens.476 . hdl : 2027/coo.31924001082803 .
  83. p. ej. Fama, E (1965). "El comportamiento de los precios del mercado de valores". Journal of Business . 38 .
  84. Hamilton, James (1989). "Un nuevo enfoque para el análisis económico de series temporales no estacionarias y el ciclo económico". Econometrica . 57 (2): 357–84 . CiteSeerX 10.1.1.397.3582 . doi : 10.2307/1912559 . JSTOR 1912559 .  
  85. Calvet, Laurent E.; Fisher, Adlai J. (2001). "Pronóstico de la volatilidad multifractal" . Journal of Econometrics . 105 (1): 27– 58. Bibcode : 2001JEcon.105...27C . doi : 10.1016/S0304-4076(01)00069-0 .
  86. Calvet, Laurent; Adlai Fisher (2004). "Cómo pronosticar la volatilidad a largo plazo: cambio de régimen y estimación de procesos multifractales". Journal of Financial Econometrics . 2 : 49–83 . CiteSeerX 10.1.1.536.8334 . doi : 10.1093/jjfinec/nbh003 . 
  87. Brennan, Michael; Xiab, Yihong. "Volatilidad del precio de las acciones y prima de riesgo de las acciones" (PDF) . Departamento de Finanzas, Escuela de Administración Anderson, UCLA . Archivado del original (PDF) el 28 de diciembre de 2008.
  88. "Un ejemplo de cadena de Markov en la modelización del riesgo crediticio" (PDF) . Universidad de Columbia . Archivado del original (PDF) el 24 de marzo de 2016.
  89. Acemoglu, Daron; Georgy Egorov; Konstantin Sonin (2011). "Modelo político de evolución social" . Actas de la Academia Nacional de Ciencias . 108 ( Supl. 4): 21292–21296 . Bibcode : 2011PNAS..10821292A . CiteSeerX 10.1.1.225.6090 . doi : 10.1073/pnas.1019454108 . PMC 3271566. PMID 22198760 .   
  90. K McAlpine; E Miranda; S Hoggar (1999). "Creación musical con algoritmos: un sistema de estudio de caso". Computer Music Journal . 23 (2): 19– 30. doi : 10.1162/014892699559733 .
  91. Curtis Roads, ed. (1996). The Computer Music Tutorial . MIT Press. ISBN 978-0-262-18158-7.
  92. Xenakis, Iannis; Kanach, Sharon (1992) Música formalizada: Matemáticas y pensamiento en la composición , Pendragon Press. ISBN 1576470792
  93. "Continuador" . Archivado del original el 13 de julio de 2012.
  94. Pachet, F.; Roy, P.; Barbieri, G. (2011) "Procesos de Markov de longitud finita con restricciones" Archivado el 14 de abril de 2012 en Wayback Machine , Actas de la 22.ª Conferencia Internacional Conjunta sobre Inteligencia Artificial , IJCAI, páginas 635–642, Barcelona, ​​España, julio de 2011
  95. Pankin, Mark D. "MODELOS DE CADENAS DE MARKOV: FUNDAMENTOS TEÓRICOS" . Archivado del original el 9 de diciembre de 2007. Consultado el 26 de noviembre de 2007 .
  96. Pankin, Mark D. "EL BÉISBOL COMO UNA CADENA DE MARKOV" . Archivado del original el 13 de mayo de 2001. Consultado el 24 de abril de 2009 .
  97. "El rincón del poeta – Fieralingue" . Archivado del original el 6 de diciembre de 2010.
  98. Kenner, Hugh; O'Rourke, Joseph (noviembre de 1984). "Un generador de parodias para micros". BYTE . 9 (12): 129– 131, 449– 469.
  99. Hartman, Charles (1996). Virtual Muse: Experiments in Computer Poetry . Hanover, NH: Wesleyan University Press. ISBN 978-0-8195-2239-9.

Referencias

  • AA Markov (1906) "Rasprostranenie zakona bol'shih chisel na velichiny, zavisyaschie drug ot druga". Izvestiya Fiziko-matematicheskogo obschestva pri Kazanskom universitete , 2-ya seriya, tom 15, págs  .
  • AA Markov (1971). «Extensión de los teoremas límite de la teoría de la probabilidad a una suma de variables conectadas en una cadena». Reimpreso en el Apéndice B de: R. Howard. Sistemas probabilísticos dinámicos, volumen 1: Cadenas de Markov . John Wiley and Sons.
  • Texto clásico en traducción: Markov, AA (2006). «Un ejemplo de investigación estadística del texto Eugenio Oneguin sobre la conexión de muestras en cadenas». Science in Context . 19 (4). Traducido por Link, David: 591– 600. doi : 10.1017/s0269889706001074 .
  • Leo Breiman (1992) [1968] Probabilidad . Edición original publicada por Addison-Wesley; reimpresa por la Society for Industrial and Applied Mathematics ISBN 0-89871-296-3(Véase el capítulo 7)
  • JL Doob (1953) Procesos estocásticos . Nueva York: John Wiley and Sons ISBN 0-471-52369-0.
  • SP Meyn y RL Tweedie (1993) Cadenas de Markov y estabilidad estocástica . Londres: Springer-Verlag ISBN 0-387-19832-6. en línea: MCSS . Segunda edición de próxima publicación, Cambridge University Press, 2009.
  • Dynkin, Evgeny Borisovich (1965). Procesos de Markov . Grundlehren der mathematischen Wissenschaften. vol.  Yo (121). Traducido por Fabio, Jaap; Greenberg, Vida Lázaro; Maitra, Ashok Prasad; Majone, Giandomenico . Berlín: Springer-Verlag . doi : 10.1007/978-3-662-00031-1 . ISBN 978-3-662-00033-5Título n.° 5104.; Procesos de Markov . Grundlehren der mathematischen Wissenschaften. vol. II (122). 1965. doi : 10.1007/978-3-662-25360-1 . ISBN  978-3-662-23320-7. Título-Núm. 5105.(Nota: Este texto fue publicado originalmente en ruso como Марковские процессы ( Markovskiye protsessy ) por Fizmatgiz en 1963 y traducido al inglés con la ayuda del autor).
  • SP Meyn. Técnicas de control para redes complejas . Cambridge University Press, 2007. ISBN 978-0-521-88441-9El apéndice contiene una versión abreviada de Meyn y Tweedie. Disponible en línea en CTCN.
  • Booth, Taylor L. (1967). Máquinas secuenciales y teoría de autómatas (1.ª  ed.). Nueva York, NY: John Wiley and Sons, Inc. Número de catálogo de la Biblioteca del Congreso: 67-25924.Libro extenso y de amplio alcance dirigido a especialistas, tanto informáticos teóricos como ingenieros eléctricos. Incluye explicaciones detalladas de técnicas de minimización de estados, máquinas de estados finitos (FSM), máquinas de Turing, procesos de Markov e indecidibilidad. Excelente análisis de los procesos de Markov (págs.  449 y siguientes). Se abordan las transformadas Z y D en su contexto.
  • Kemeny, John G.; Hazleton Mirkil; J. Laurie Snell; Gerald L. Thompson (1959). Estructuras matemáticas finitas (1.ª  ed.). Englewood Cliffs, NJ: Prentice-Hall, Inc. Número de catálogo de la Biblioteca del Congreso: 59-12841.Texto clásico. Véase el capítulo 6, Cadenas de Markov finitas, págs.  384 y siguientes.
  • John G. Kemeny y J. Laurie Snell (1960) Cadenas finitas de Markov , D. van Nostrand Company ISBN 0-442-04328-7
  • E. Nummelin. «Cadenas de Markov irreducibles generales y operadores no negativos». Cambridge University Press, 1984, 2004. ISBN 0-521-60494-X
  • Seneta, E. Matrices no negativas y cadenas de Markov . 2.ª ed. revisada, 1981, XVI, 288 p., Tapa blanda. Serie Springer en Estadística. (Publicado originalmente por Allen & Unwin Ltd., Londres, 1973) . ISBN 978-0-387-29765-1
  • Kishor S. Trivedi , Probabilidad y estadística con aplicaciones en fiabilidad, teoría de colas e informática , John Wiley & Sons, Inc., Nueva York, 2002. ISBN 0-471-33341-7.
  • KS Trivedi y RASahner, SHARPE a los veintidós años , vol. 36, n.º 4, págs.  52-57, ACM SIGMETRICS Performance Evaluation Review, 2009.
  • RA Sahner, KS Trivedi y A. Puliafito, Análisis del rendimiento y la fiabilidad de los sistemas informáticos: un enfoque basado en ejemplos utilizando el paquete de software SHARPE , Kluwer Academic Publishers, 1996. ISBN 0-7923-9650-2.
  • G. Bolch, S. Greiner, H. de Meer y KS Trivedi, Redes de colas y cadenas de Markov , John Wiley, 2.ª edición, 2006. ISBN 978-0-7923-9650-5.
  • "Cadena de Markov" , Enciclopedia de Matemáticas , EMS Press , 2001 [1994]
  • Capítulo sobre cadenas de Markov en el libro de introducción a la probabilidad de la Sociedad Matemática Americana. Archivado el 22 de mayo de 2008 en la Wayback Machine.
  • Introducción a las cadenas de Markov en YouTube
  • Una explicación visual de las cadenas de Markov
  • Artículo original de AA Markov (1913): Un ejemplo de investigación estadística del texto Eugenio Oneguin sobre la conexión de muestras en cadenas (traducido del ruso).