Cadenas de Markov
Cadenas de Markov
Cadenas de Markov
MARKOV
INTRODUCCIN
t
j1
ij 1, con 0 tij 1
En los modelos ms simples de cadenas de Markov, los valores de las
probabilidades de transicin tij no dependen ni de cmo el sistema lleg al estado i,
ni del periodo n. Las probabilidades de ocupar un estado i dependen del nmero de
periodos o de transiciones efectuadas.
Por lo tanto una secuencia de intentos de un experimento es una cadena de
Markov.
Algunas veces nos interesa saber cmo cambia una variable aleatoria a travs
del tiempo. Por ejemplo, desearamos conocer cmo evoluciona el precio de las
acciones de una empresa en el mercado a travs del tiempo. El estudio de cmo
evoluciona una variable aleatoria incluye el concepto de procesos
estocsticos. En este captulo explicaremos esos procesos, en especial uno que
se conoce como cadena de Markov. Las cadenas de Markov se han aplicado en
reas tales como educacin, mercadotecnia, servicios de salud, contabilidad y
produccin.
Indica que la ley de probabilidad que relaciona el estado del siguiente perodo con
el estado actual no cambia, o que permanece estacionaria en el tiempo.
Toda cadena de Markov que cumple con esta condicin se llama cadena
estacionaria de Markov.
q= 1
j=1,n j
IMPORTANCIA:
APLICACIONES
Qu es un Proceso Estocstico?
Ejemplo
Xt = Condiciones de tiempo en una ciudad
Ejemplo
En los casos de la condicin del clima y del precio de venta de la carne de aves, su
estado en un tiempo determinado no es aleatorio por completo, sino que es afectado
en cierto grado por el tiempo de das previos. No obstante la sucesin es tan
compleja que dicho comportamiento cambia da a da de una manera que en
apariencia es algo aleatorio.
Ejemplo
Xt = Temperatura registrada en cada da del ao en la ciudad de Arequipa.
Ejemplo
Xt = Cantidad de alumnos que egresan semestralmente en la especialidad de
ingeniera industrial.
Ejemplo
Grfico:
Matriz de transicin:
Ejemplo
Supongamos que el estado de nimo de Enzo puede ser, alegre, triste o molesto.
Si est alegre, habr una probabilidad de 80% de que siga alegre al da siguiente
y 15% de que est triste. Si est triste, 30% de que siga triste y 50% de que est
molesto. Y, si est molesto, 80% de que no siga molesto y 50% de que est alegre
al da siguiente. Hallar la Matriz de Transicin.
Grfico:
Matriz de transicin:
Ejemplo:
Grfico:
Matriz de transicin:
Ejemplo:
Lnea telefnica
Grfico:
Matriz de transicin:
Se lanza un dado repetidas veces. Cada vez que sale menor que 5 se pierde 1 ,
y cada vez que sale 5 o 6 se gana 2 . El juego acaba cuando se tienen 0 o 100
.
Sea Xt=estado de cuentas en el instante t. Tenemos que { Xt } es una CM
S={0, 1, 2, , 100}
Grfico:
Ejemplo: Lanzamiento de una moneda
Lanzamos una moneda al aire 6 veces. El jugador gana 1 cada vez que sale
cara (C), y pierde 1 cada vez que sale cruz (F).
Xi = estado de cuentas del jugador despus de la i-sima jugada
La familia de variables aleatorias {X1, X2,, X6} constituye un proceso
estocstico
={CCCCCC,CCCCCF,}
Combinaciones () = 26 = 64
P()=1/64
T={1, 2, 3, 4, 5, 6}
S={6, 5, , 1, 0, 1, 2, , 5, 6}
Si fijo t, por ejemplo t0=3, obtengo una de las variables aleatorias del proceso:
X3 :
X 3
Los posibles valores que puede tomar el proceso en :X3()={3, 1, 1, 3}
APLICACIONES
Aplicacin 1:
Participaciones de mercado
Se desea saber:
a) El porcentaje de los clientes que consumen cada marca de cerveza despus
de un periodo.
b) El porcentaje de los clientes que consumen cada marca de cerveza despus
de 2 periodos.
c) A la larga cmo se reparte el mercado de bebedores de cerveza entre las
tres marcas?
Grficamente se tiene:
N
pj(1)= pi(0). tij
i 1
X3 = X2 T = X0T2 T = X0T3
Despus de n transiciones se tiene:
adems de:
A + B + C =1
El sistema de ecuaciones sera:
A + B + C = 1
0,7A + 0,1B + 0,24C = A
0,2A + 0,5B + 0,16C = B
0,1A + 0,4B + 0,6C = C
Este sistema es redundante y, para resolverlo, eliminamos una de las tres ltimas
ecuaciones (por ejemplo la ltima).
A + B + C = 1
-0,3A + 0,1B + 0,24C = 0
0,2A - 0,5B + 0,16C = 0
y la solucin es: = [ 0,376 0,265 0,359]
Observamos el aumento en la participacin de A, que pasa de 20% a 37.6%;
principalmente a costa de C, que cae de 50% a 35,9%. Entonces si C quiere
promover una campaa publicitaria para quebrar el proceso, debera,
principalmente, dirigirla hacia los actuales consumidores de A, ya que t AC=0,1 (muy
pequeo). Se observa que tBC es bastante grande.
La siguiente tabla muestra las distribuciones de estado para diferentes periodos de
transicin:
Se observa que a partir del periodo 7, las variaciones en los tres estados son casi
despreciables.
Ejemplos:
1 / 2 1 / 2
a) T= es regular ya que tij> 0
1 / 3 2 / 3
1 / 2 1 / 2 2 1 / 4 3 / 4
b) T= ,T = ,
0 1 0 1
1 / 8 7 / 8
T3= ,
0 1
entonces cualquier Tm no cumplir la condicin tij> 0 ya que siempre existir el
elemento t21= 0, por lo tanto T no es regular.
0 1 / 2 1 / 2 1 / 2 1 / 4 1 / 4
T= 1 / 2 0 1 / 2 , T2= 1 / 4 1 / 2 1 / 4 ,
1 / 2 1 / 2 0 1 / 4 1 / 4 1 / 2
Las probabilidades del estado del tiempo para la ciudad de Arequipa el mes de
enero del presente ao fueron extradas de la base de datos de SENAMHI Arequipa.
Dichos valores han sido tomados de un da anterior al da que iran a ser puestos a
conocimiento de la poblacin arequipea. Estos datos se pueden expresar mediante
la siguiente matriz de transicin:
Grficamente:
Pronosticando el clima.-
Xn = X(n-1)P o
Xn = XoPn
Xn = X(n-1)
Asi que;
0.9SO +0.5LL = SO - 0.1SO + 0.5LL = 0
0.1SO + 0.5LL = LL 0.1SO -0.5LL = 0
Adems SO + LL = 1
Resolviendo:
SO = 0.8333 y LL = 0.1667
Grficamente:
Usando diagramas de rbol:
Matricialmente
Grficamente:
pij (0) = 1, si i = j
pij (0) = 0, si i j
En general,
Ejemplo
Solucin
Grfico:
Matriz:
Definicin:
Definicin:
Definicin:
Definicin:
Definicin:
Un estado i es peridico con perodo k > 1, si k es el menor nmero tal que todas
las trayectorias que parten del estado i y regresan al estado i tienen una longitud
mltiplo de k.
Si k>1, entonces diremos que i es peridico de periodo j. El estado j ser peridico
de periodo k>1 si existen caminos que llevan desde j hasta i pero todos tienen
longitud mk, con m>0
Periodicidad
Definicin:
Cadenas ergdicas
ij = 1/i
ij = 1 + k j pik kj
Ejemplo
Solucin:
La matriz ser:
1 = 0.871 + 0.572
2 = 0.131 + 0.432
1 + 2 = 1
De donde:
1 = 57/70
2 = 13/70
Por tanto:
ij = 1 + k j pik kj
12 = 1 + p11 12
= 1 + 0.8712
21 = 1 + p22 21
= 1 + 0.4321
Dnde:
12 = 7.69 horas
21 =1.75 horas
El costo promedio a largo plazo, por unidad de tiempo est dado por:
g = j=1,s (j Cj)
Done:
Cj = Inventarios, Nivel de ventas, Mquinas malogradas, etc.
As tenemos:
En conclusin:
Aplicacin 1:
Planificacin de Personal
Solucin:
Los dos ltimos estados (a1, a2) son estados absorbentes y los dems son
transitorios (t1, t2, t3). Por ejemplo, Experimentado es estado transitorio, porque hay
una trayectoria de Experimentado a Sale sin ser socio, pero no hay trayectoria que
regrese de Sale sin ser socio ha Experimentado. Suponemos que una vez que un
abogado sale de la empresa nunca regresa.
Para dar respuesta a las preguntas formuladas anteriormente, es necesario obtener
las matrices: (I-Q)-1 y (I-Q)-1R, la informacin contenida en estas matrices
debidamente interpretadas, permite tomar decisiones.
Entonces s=5, m=2, y
Luego:
Interpretacin:
Muchas empresas, como por ejemplo Los Justicieros del ejemplo de planificacin
de personal, emplean varias categoras de personal. Con fines de planeacin a largo
plazo, a menudo es til poder predecir el nmero de empleados de cada categora
que, si las tendencias actuales continan, estarn disponibles en el estado estable.
Si existe censo de estado estable podemos encontrarlo al resolver un sistema de S
ecuaciones que se plantea como sigue: tan slo ntese que para que exista ese
estado, debe ser vlido, para i = 1, 2 .. S.
Nmero de personas que entran al grupo i durante cada periodo = Nmero de
personas que salen del grupo i durante cada periodo
Ejemplo: Regresemos al bufete de abogados Los Justicieros (Ejemplo anterior)
Supongamos que la meta a largo plazo de ese bufete es tener 50 abogados
principiantes, 30 con experiencia y 10 socios. Para alcanzar este censo de estado
estable, cuntos abogados de cada tipo deben contratar cada ao?
Solucin: Sean
H1= nmero de abogados principiantes a contratar
H2 = nmero de abogados con experiencia a contratar
H3 = nmero de abogados asociados a contratar
Entonces:
Generando el grfico de Markov
Aplicacin 3
Solucin:
Con los datos del problema, diferenciamos cada una de las matrices
Graficamos
Definimos las matrices absorbente (R) y transitoria (Q)
La concesin ser:
SOLUCIN:
Siendo los estados: F (funciona) y NF (no funciona), elaboramos las matrices de
transicin de estados respectivas.
Matriz de transicin de estados (T) para la Mquina 1
Luego hallamos los vectores de estado estable para ambas mquinas aplicando la
relacin:
= .T
Siendo = [ F NF ]
Dnde :
F: probabilidad de estado estable de que la mquina Funcione
NF: probabilidad de estado estable de que la mquina No Funcione
Adems F + NF =1
Para la mquina 1 tenemos:
[ F NF ]= [ F NF ]*T
F + NF =1
Reemplazando los datos de matriz de Transicin de estados de la mquina 1 y
resolviendo el sistema de ecuaciones tenemos:
F= 0.9375 y NF=0.0625
entonces 1= [ F NF ]= [ 0.9375 0.0625 ]
Para la mquina 2 tenemos 2 =[ F NF ]= [ 0.8889 0.1111 ]
Por lo tanto se observa que la mquina 1 tiene mayor probabilidad de
funcionamiento (93.75%) frente a 88.89% de la mquina 2, en consecuencia la
empresa debe rentar la mquina 1.
APLICACIN 5.- Una pequea tienda de videos lleva un control del nmero de
veces por semana que es rentado un video y estima las siguientes probabilidades
de transicin (ver matriz siguiente).
Grfico:
Dnde: los estados en orden son: 5 veces, 4 veces, 3 veces, 2 veces, 1 vez y 0
veces
Por ejemplo, si un video se rent 5 veces esta semana, entonces hay una
probabilidad de 80% de que sea rentado 5 veces la siguiente semana, 10% de
probabilidades de que sea rentado 4 veces y 10% de probabilidades de que sea
rentado 3 veces. Cuando un video es rentado 0 veces, este se desecha.
Se determina la matriz fundamental y expresan el nmero de semanas que sern
rentados ( 5, 4, 3, 2,1) veces:
APLICACIN 6.- Juan es propietario de un terreno con 5000 pinos. Cada ao Juan
permite a los detallistas de rboles de navidad seleccionar y cortar rboles para la
venta a clientes individuales. Juan protege los rboles pequeos (por lo general de
menos de 120 cm. de alto) de manera que estn disponibles para la venta en aos
futuros. Actualmente estn clasificados 1500 rboles como protegidos, en tanto que
los 3500 restantes estn disponibles para corte. Sin embargo, aunque en un ao
dado un rbol est disponible para corte, quizs no sea seleccionado sino hasta en
aos futuros. Aunque la mayora de los rboles que no se cortan en un ao dado
viven hasta el siguiente, todos los aos se pierden algunos pinos enfermos.
Al estudiar la operacin de los rboles de navidad de Juan como un proceso de
Markov con periodos anuales, definimos los cuatro estados siguientes:
Estado 1. Cortado y vendido.
Estado2. Perdido por enfermedad.
Estado3. Pequeo para cortarse
Estado4. Disponible para cortar, pero no cortado ni vendido
La siguiente matriz de transicin es apropiada
Grfico:
80% de probabilidad de que un rbol disponible para cortar sea cortado y vendido y
20% de probabilidad de que un rbol disponible para cortar se pierda por
enfermedad.
Estado Descripcin
1 Artculo en mquina 1
2 Artculo en inspeccin
1
3 Artculo en mquina 2
4 Artculo en inspeccin
2
5 Artculo en mquina 3
6 Artculo en inspeccin
3
7 Artculo en almacn
8 Artculo desechado
Grfico:
Matriz:
d. Con los datos de costos de mano de obra y de materiales, cules son los
costos directos esperados?
Resolviendo:
PE = PH = 0.2
PC = PS = 0.3
As que para un periodo de 5 horas (300 minutos)
tE = tH = 60 min
tC = tS = 90 min
Este tipo de distribucin es lmite ya que existe una sola clase final y es aperidica.
Debido a que ahora tenemos una nueva puerta debemos tener en cuenta la
aparicin de un nuevo estado (SL), genera el grfico siguiente:
En esta nueva situacin no tenemos una distribucin lmite ya que existen dos
posibles estados finales que pueden acabar con el ratn. Adems el estado inicial
va a repercutir en el resultado del problema ya que si el ratoncito se encuentra en
un principio en la cocina es ms probable que muera envenenado en la cocina,
mientras que si se encuentra en la entrada tiene ms posibilidades de salir de la
casa.
Para ello nos planteamos que los estados E, S, D son estados transitorios, mientras
que C y SL son estados absorbentes o finales.
La matriz de probabilidad:
APLICACIONES PROPUESTAS
APLICACIN 6.- Freezco, Inc., vende refrigeradores. La fbrica otorga una garanta
en todos los refrigeradores que especifica cambio gratis de cualquier unidad que se
descomponga antes de tres aos. Se nos da la siguiente informacin: (1) el 3% de
todos los refrigeradores nuevos falla durante su primer ao de funcionamiento; (2)
el 5% de todos los refrigeradores con 1 ao de funcionamiento falla durante el
segundo ao de trabajo, y (3) el 7% de todos los refrigeradores con dos aos de
funcionamiento falla durante su tercer ao. La garanta no vale para el refrigerador
de repuesto.
a) Use la teora de cadenas de Markov para predecir la fraccin de todos los
refrigeradores que deber cambiar Freezco.
b) Suponga que a Freezco le cuesta 500 dlares cambiar un refrigerador y que
vende 10000 refrigeradores al ao. Si la fbrica redujera el plazo de garanta a dos
aos, cunto dinero se ahorrara en costos de reemplazo?
APLICACIN 11.- El estado de las cuentas por cobrar en una empresa se modela
con frecuencia como una cadena absorbente de Markov. Suponga que una empresa
supone que una cuenta es incobrable si han pasado ms de tres meses de su fecha
de vencimiento. Entonces, al principio de cada mes, se puede clasificar cada cuenta
en uno de los siguientes estados especficos:
Estado 1 Cuenta nueva.
Estado 2 Los pagos de la cuenta estn retrasados un mes.
Estado 3 Los pagos de la cuenta estn retrasados dos meses.
Estado 4 Los pagos de la cuenta estn retrasados tres meses.
Estado 5 Se ha saldado una cuenta.
Estado 6 Se ha cancelado la cuenta por ser mal pagador.
Supongamos que los ltimos datos indican que la siguiente cadena de Markov
describe cmo cambia el estado de una cuenta de un mes al siguiente:
Por ejemplo si al principio de un mes una cuenta lleva dos meses de vencida, hay
40% de probabilidades de que no se pague al principio del mes siguiente y, por lo
tanto, que tenga tres meses de retraso y una probabilidad de 60% de que se pague.
Suponga ademn que despus de tres meses, la cuenta o se cobra o se considera
incobrable.
Una vez que una deuda se paga o se considera incobrable, se cierra y no se tiene
ms transiciones.
a) Cul es la probabilidad que una cuenta nueva sea cobrada alguna vez?
b) Cul es la probabilidad que una cuenta atrasada un mes se vuelva
finalmente incobrable?
c) Si las ventas de la empresa son 100 000 dlares en promedio mensual,
cunto dinero ser incobrable cada ao?
APLICACIN 13.- El 1 de enero (de este ao), las panaderas Klosman controlaban
el 40% de su mercado local, mientras que las otras dos panaderas, A y B, tenan
40 y 20 por ciento, respectivamente, del mercado. Basndose en un estudio de una
empresa de investigacin de mercado, se compilaron los siguientes datos: la
panadera Klosman retiene el 90% de sus clientes, y gana el 5% de los clientes de
A y el 10% de los de B. La panadera A retiene el 85% de sus clientes y gana 5%
de los clientes de Klosman y 7% de los de B. La panadera B retiene 83% de sus
clientes y gana 5% de los clientes de Klosman y 10% de los de A.
a) Cul ser la participacin de cada empresa en 1 de enero del ao
siguiente?
b) Klosman decide hacer una campaa publicitaria a efectos de ganar clientes,
dicha campaa altera las probabilidades de transicin de estados de la siguiente
manera: la panadera Klosman retiene el 90% de sus clientes, y gana el 15% de los
clientes de A y el 20% de los de B. La panadera A retiene el 75% de sus clientes y
gana 5% de los clientes de Klosman y 7% de los de B. La panadera B retiene 73%
de sus clientes y gana 5% de los clientes de Klosman y 10% de los de A. Si a
Klosman le cuesta 350 dlares por mes una campaa publicitaria y por cada cliente
ganado obtiene un ingreso igual a 10 dlares mensuales, por cuntos periodos
debe mantener su campaa publicitaria, sabiendo que se compite en un mercado
de 1000 clientes?