Guía Completa para Resolver Problemas de Estados Absorbentes en Cadenas de Markov

Guía Completa para Resolver Problemas de Estados Absorbentes en Cadenas de Markov

Nota: para usar esta calculadora, la matriz de transicion debe estar debidamente organizada, en este orden:

Calculadora de Cadenas de Markov

Calculadora de Cadenas de Markov con Estados Absorbentes

Introduce la matriz de transición P (dimensiones de 2×2 hasta 6×6):

Las cadenas de Markov con estados absorbentes son una herramienta fundamental en la modelación de procesos estocásticos, en los que ciertos estados actúan como "trampas" en las que, una vez ingresados, el proceso no puede salir. En esta guía te presentamos 4 pasos esenciales para resolver este tipo de problemas, ilustrándolos con ejemplos prácticos.

Tabla de Contenidos:

Estados Absorbentes

Una cadena de Markov describe la evolución de un sistema en pasos discretos, en el que la probabilidad de transitar a otro estado depende únicamente del estado actual. En algunos casos, el sistema contiene estados absorbentes: estados que, al alcanzarlos, el proceso se detiene.

Para resolver problemas de este tipo se deben calcular dos matrices clave:

  • La matriz fundamental N = (I − Q)−1, que indica el número esperado de visitas a los estados transitorios antes de ser absorbido.
  • La matriz de absorción P = N · R, que muestra la probabilidad de terminar en cada estado absorbente.

A continuación, se detalla el proceso en 4 pasos, acompañado de ejemplos.

1

Identificar los estados absorbentes

Definición: Un estado i es absorbente si la probabilidad de permanecer en él una vez que se entra es 1, es decir, Pii = 1 y Pij = 0 para todo j ≠ i.

Procedimiento:

  • Inspecciona cada fila de la matriz de transición P.
  • Si encuentras una fila cuyo valor en la diagonal es 1 y el resto de valores es 0, ese estado es absorbente.
2

Reordenar la matriz y extraer Q, R, O, I

Para aplicar la teoría, se reordena la matriz de transición de modo que primero aparezcan todos los estados transitorios y después todos los absorbentes. En esa forma canónica, la matriz queda dividida en cuatro bloques:

P reordenadaEstados transitoriosEstados absorbentes
Estados transitoriosQR
Estados absorbentesOI
  • Q (dimensión t×t, t = nº de estados transitorios): probabilidades de ir de un estado transitorio a otro estado transitorio.
  • R (dimensión t×m, m = nº de estados absorbentes): probabilidades de ir de un estado transitorio a un estado absorbente.
  • O (dimensión m×t): matriz de ceros — probabilidades de ir de un estado absorbente a uno transitorio. Siempre es cero porque de un estado absorbente no se sale.
  • I (dimensión m×m): matriz identidad — probabilidades de ir de un estado absorbente a sí mismo (una vez absorbido, el proceso no cambia).
3

Calcular la matriz fundamental N = (I − Q)−1

Definición: N = (I − Q)−1 se conoce como la matriz fundamental. Cada elemento Nij indica el número esperado de veces que el proceso estará en el estado transitorio j si empieza en el estado transitorio i, antes de ser absorbido.

Te puede interesar:  Criterio Maximin o Wald: Guía Definitiva para Decisiones Seguras sin Riesgos

Cálculo (método de Gauss-Jordan):

  • Se construye la matriz (I − Q).
  • Se aplica Gauss-Jordan para hallar su inversa.
  • El resultado es N.

Si Q es de dimensión t×t, entonces (I − Q) también es t×t, y su inversa N mantiene esa misma dimensión.

4

Calcular la matriz de absorción P = N · R

Una vez obtenida N, calculamos la probabilidad de que, partiendo de un estado transitorio i, el proceso sea absorbido en cada estado absorbente:

P = N · R

  • Cada fila de P corresponde a un estado transitorio.
  • Cada columna de P corresponde a un estado absorbente.
  • Por lo tanto, Pij es la probabilidad de que, empezando en el estado transitorio i, el proceso termine en el estado absorbente j.

ℹ️ Si hay más de un estado absorbente, P da la distribución de probabilidades entre ellos. Si solo hay un estado absorbente, P es un vector columna de unos (100% de probabilidad).

Ejemplo 1 — Matriz 3×3 (1 estado absorbente)

Sistema con 3 estados donde el Estado 3 es la única "trampa"

Considera un sistema con 3 estados (E1, E2, E3). La siguiente matriz de transición describe el comportamiento del proceso:

EstadoE1E2E3
E10.50.30.2
E20.40.20.4
E3 (absorbente)001
1

Identificar los estados absorbentes

  • Fila E1: P11 = 0.5 ≠ 1 → Transitorio
  • Fila E2: P22 = 0.2 ≠ 1 → Transitorio
  • Fila E3: P33 = 1 y P31 = P32 = 0 → Absorbente ✓

✅ El Estado 3 es el único estado absorbente. E1 y E2 son estados transitorios.

2

Extraer las submatrices Q, R, O, I

La matriz ya está ordenada (transitorios primero). Identificamos los bloques:

Submatriz Q — Transitorios → Transitorios (2×2)

QE1E2
E10.50.3
E20.40.2

Submatriz R — Transitorios → Absorbentes (2×1)

RE3
E10.2
E20.4

O = [0   0] (absorbente → transitorios, siempre cero)  ·  I = [1] (absorbente → sí mismo)

3

Calcular la Matriz Fundamental N = (I − Q)−1

a) Construimos (I − Q):

I − Q =
10
01
0.50.3
0.40.2
=
0.5−0.3
−0.40.8

b) Invertimos por el método de Gauss-Jordan. Construimos la matriz aumentada [I − Q | I] y aplicamos operaciones de fila hasta que la parte izquierda sea la identidad. La parte derecha resultante es N.

Matriz aumentada inicial
0.5−0.310
−0.40.801

Izquierda = (I − Q)  |  Derecha = Identidad

Paso 1Fila 1 × 2 → para obtener 1 en posición [1,1]
1−0.620
−0.40.801

¿Por qué ×2? El valor en la posición [1,1] es 0.5. Para convertirlo en 1, necesitamos multiplicarlo por su inverso: 1 ÷ 0.5 = 2. Se multiplica toda la fila por ese mismo factor.
Resultado: 0.5×2 = 1  ·  −0.3×2 = −0.6  ·  1×2 = 2  ·  0×2 = 0

Paso 2Fila 2 → Fila 2 + 0.4 × Fila 1 → para obtener 0 en posición [2,1]
1−0.620
00.560.81

¿Por qué +0.4 × Fila 1? El valor en la posición [2,1] es −0.4. Para anularlo necesitamos sumarle su opuesto, es decir, +0.4. Como en la posición [1,1] ya hay un 1, multiplicar Fila 1 por 0.4 y sumárselo a Fila 2 produce exactamente: −0.4 + 0.4×1 = 0.
Resto de la fila: 0.8 + 0.4×(−0.6) = 0.56  ·  0 + 0.4×2 = 0.8  ·  1 + 0.4×0 = 1

Te puede interesar:  Predicción a Corto Plazo con Cadenas de Markov: Una Guía Paso a Paso
Paso 3Fila 2 ÷ 0.56 → para obtener 1 en posición [2,2]
1−0.620
011.4291.786

¿Por qué ÷0.56? El valor en la posición [2,2] es 0.56. Para convertirlo en 1, dividimos toda la fila entre ese mismo valor (es decir, multiplicamos por 1÷0.56).
Resultado: 0.56÷0.56 = 1  ·  0.8÷0.56 = 1.429  ·  1÷0.56 = 1.786

Paso 4Fila 1 → Fila 1 + 0.6 × Fila 2 → para obtener 0 en posición [1,2]
102.8571.071
011.4291.786

¿Por qué +0.6 × Fila 2? El valor en la posición [1,2] es −0.6. Para anularlo necesitamos sumarle +0.6. Como en la posición [2,2] ya hay un 1, multiplicar Fila 2 por 0.6 y sumárselo a Fila 1 produce: −0.6 + 0.6×1 = 0.
Resto de la fila: 2 + 0.6×1.429 = 2.857  ·  0 + 0.6×1.786 = 1.071

Resultado: la parte derecha es N
102.8571.071
011.4291.786

Parte izquierda = identidad ✓  ·  Parte derecha (verde) = N = (I − Q)⁻¹

NE1E2
E12.8571.071
E21.4291.786

📌 Interpretación: Partiendo de E1, se esperan ~2.86 visitas a E1 y ~1.07 visitas a E2 antes de ser absorbido. Partiendo de E2, se esperan ~1.43 visitas a E1 y ~1.79 visitas a E2.

4

Calcular la Matriz de Absorción P = N · R

Multiplicamos N (2×2) por R (2×1):

P = N · R =
2.8571.071
1.4291.786
×
0.2
0.4

Fila E1:   2.857 × 0.2  +  1.071 × 0.4  =  0.571 + 0.428  =  1.000


Fila E2:   1.429 × 0.2  +  1.786 × 0.4  =  0.286 + 0.714  =  1.000

PP(absorción en E3)
Desde E11.000 = 100%
Desde E21.000 = 100%

Conclusión: Con un solo estado absorbente, la probabilidad de terminar en E3 es del 100% sin importar el estado de inicio. Esto confirma que todo estado transitorio termina inevitablemente en el absorbente.

Ejemplo 2 — Matriz 4×4 (2 estados absorbentes)

Caso de negocio: clientes que pueden terminar en quiebra o fidelización

Una empresa clasifica a sus clientes en dos categorías de riesgo (E1, E2) que pueden terminar en quiebra (E3) o en fidelización permanente (E4):

EstadoE1E2E3 (quiebra)E4 (fideliz.)
E10.40.20.30.1
E20.10.50.20.2
E3 (absorbente)0010
E4 (absorbente)0001
1

Identificar los estados absorbentes

  • E1: P11 = 0.4 ≠ 1 → Transitorio
  • E2: P22 = 0.5 ≠ 1 → Transitorio
  • E3: P33 = 1, resto = 0 → Absorbente ✓
  • E4: P44 = 1, resto = 0 → Absorbente ✓

✅ Hay 2 estados absorbentes (E3 y E4) y 2 estados transitorios (E1 y E2).

2

Extraer las submatrices Q, R, O, I

La matriz ya está ordenada. Extraemos los cuatro bloques:

Submatriz Q (2×2) — Transitorios → Transitorios

QE1E2
E10.40.2
E20.10.5

Submatriz R (2×2) — Transitorios → Absorbentes

RE3E4
E10.30.1
E20.20.2

O = matriz 2×2 de ceros (absorbentes → transitorios)  ·  I = identidad 2×2 (absorbentes → sí mismos)

3

Calcular la Matriz Fundamental N = (I − Q)−1

a) Construimos (I − Q):

I − Q =
10
01
0.40.2
0.10.5
=
0.6−0.2
−0.10.5

b) Invertimos por el método de Gauss-Jordan. Construimos la matriz aumentada [I − Q | I]:

Matriz aumentada inicial
0.6−0.210
−0.10.501
Paso 1Fila 1 ÷ 0.6 → para obtener 1 en posición [1,1]
1−0.3331.6670
−0.10.501
Te puede interesar:  Cadenas de Markov: Guía Completa y Práctica para Entenderlas Fácilmente

¿Por qué ÷0.6? El valor en la posición [1,1] es 0.6. Para convertirlo en 1, dividimos toda la fila entre ese mismo valor (es decir, multiplicamos por 1÷0.6 ≈ 1.667).
Resultado: 0.6÷0.6 = 1  ·  −0.2÷0.6 = −0.333  ·  1÷0.6 = 1.667

Paso 2Fila 2 → Fila 2 + 0.1 × Fila 1 → para obtener 0 en posición [2,1]
1−0.3331.6670
00.4670.1671

¿Por qué +0.1 × Fila 1? El valor en la posición [2,1] es −0.1. Para anularlo necesitamos sumarle +0.1. Como en [1,1] ya hay un 1, multiplicar Fila 1 por 0.1 y sumárselo a Fila 2 produce: −0.1 + 0.1×1 = 0.
Resto: 0.5 + 0.1×(−0.333) = 0.467  ·  0 + 0.1×1.667 = 0.167  ·  1 + 0.1×0 = 1

Paso 3Fila 2 ÷ 0.467 → para obtener 1 en posición [2,2]
1−0.3331.6670
010.3572.143

¿Por qué ÷0.467? El valor en la posición [2,2] es 0.467. Para convertirlo en 1, dividimos toda la fila entre ese mismo valor.
Resultado: 0.467÷0.467 = 1  ·  0.167÷0.467 = 0.357  ·  1÷0.467 = 2.143

Paso 4Fila 1 → Fila 1 + 0.333 × Fila 2 → para obtener 0 en posición [1,2]
101.7860.714
010.3572.143

¿Por qué +0.333 × Fila 2? El valor en la posición [1,2] es −0.333. Para anularlo necesitamos sumarle +0.333. Como en [2,2] ya hay un 1, multiplicar Fila 2 por 0.333 y sumárselo a Fila 1 produce: −0.333 + 0.333×1 = 0.
Resto: 1.667 + 0.333×0.357 = 1.786  ·  0 + 0.333×2.143 = 0.714

Resultado: la parte derecha es N
101.7860.714
010.3572.143

Parte izquierda = identidad ✓  ·  Parte derecha (verde) = N = (I − Q)⁻¹

NE1E2
E11.7860.714
E20.3572.143

📌 Interpretación: Un cliente en E1 visitará su estado ~1.79 veces y E2 ~0.71 veces antes de terminar. Un cliente en E2 permanecerá en E2 ~2.14 veces antes de ser absorbido.

4

Calcular la Matriz de Absorción P = N · R

Multiplicamos N (2×2) por R (2×2). El resultado P es también 2×2:

P = N · R =
1.7860.714
0.3572.143
×
0.30.1
0.20.2

Calculamos cada celda de P por separado:

Fila E1, columna E3:   1.786×0.3 + 0.714×0.2 = 0.536 + 0.143 = 0.679

Fila E1, columna E4:   1.786×0.1 + 0.714×0.2 = 0.179 + 0.143 = 0.321

Suma fila E1: 0.679 + 0.321 = 1.000 ✓


Fila E2, columna E3:   0.357×0.3 + 2.143×0.2 = 0.107 + 0.429 = 0.536

Fila E2, columna E4:   0.357×0.1 + 2.143×0.2 = 0.036 + 0.429 = 0.464

Suma fila E2: 0.536 + 0.464 = 1.000 ✓

PE3 — QuiebraE4 — Fidelización
Desde E10.679 = 67.9%0.321 = 32.1%
Desde E20.536 = 53.6%0.464 = 46.4%

Conclusión: Un cliente en E1 tiene 67.9% de probabilidad de quiebra y 32.1% de fidelización. Un cliente en E2 tiene 53.6% de quiebra y 46.4% de fidelización. Las filas suman 1, lo que confirma la corrección de todos los cálculos.

📋 Resumen del método — 4 pasos

  1. Identificar absorbentes: Pii = 1 con todos los demás valores en 0.
  2. Extraer Q, R, O e I: Ordena la matriz con transitorios primero, luego absorbentes.
  3. Calcular N = (I − Q)−1: Construye I − Q y aplica Gauss-Jordan para invertirla.
  4. Calcular P = N · R: Cada Pij es la probabilidad de absorción en el estado j partiendo de i.

🔎 Verificación: La suma de cada fila de la matriz P debe ser exactamente 1. Si no es así, hay un error en N o en R.

Arturo

Ingeniero Industrial con +20 años de experiencia en optimizar procesos y garantizar la calidad y seguridad en la industria. Fundador de aprendeindustrial.com, donde comparte conocimiento práctico para los ingenieros del futuro.

Te Puede Interesar:

Go up