25 de mayo de 2011

CADENAS DE MARKOV

Una cadena de Markov es una serie de eventos estadísticamente determinados, en la cual la probabilidad de que ocurra un evento depende del evento inmediato anterior. Estas se dice que tienen memoria. Es decir, recuerdan el último evento y al mismo tiempo condiciona las posibilidades de los eventos futuros subsecuentes. Ésta dependencia del evento anterior distingue a las cadenas de Markov de las series de eventos independientes.

En 1907, A. A. Markov comenzó el estudio de un importante nuevo tipo de cambio de proceso. En este proceso, el resultado de un experimento dado puede afectar el resultado del próximo experimento o evento. A este tipo de proceso es llamado Cadena de Markov.
Uno de los métodos usuales para exhibir las probabilidades de transición de un evento o Estado es usar una Matriz de Transición. La cual debe cumplir con las siguientes condiciones:
1. La Matriz de Transición debe sr Cuadrara, es decir debe tener el mismo número de columnas como de filas.
2. En ella deben estar contenidos tanto en las filas como en las columnas los mismos Estados o Eventos trasitorios.
3. La Suma de los elementos de cada fila debe ser siempre igual a 1, cumpliendo con la teoría de Probabilidades.

4. Cada elemento de la matriz debe ser un número entre cero y 1.
A continuación se muestra la matriz de transición:




En esta están contenidos M estados transitorios con una P de probabilidad.
Las aplicaciones de las Cadenas de Markov son muy variadas, ya que describen satisfacteriamente cualquier tipo de evento estocástico, por lo cual estas son llamadas también Estudios de Comportamiento de un Sistema. Estos durante un período dado suele ser llevar al análisis de un proceso estocástico con la siguiente estructura. En puntos específicos del tiempo t , el sistema se encuentra exactamente en una de un número finito de estados mutuamente excluyentes y exhaustivos, etiquetados 0, 1, . . , S. Los periodos en el tiempo pueden encontrarse a intervalos iguales o su esparcimiento puede depender del comportamiento general del sistema en el que se encuentra sumergido el proceso estocástico. Aunque los estados pueden constituir una caracterización tanto cualitativa como cuantitativa del sistema, no hay pérdida de generalidad con las etiquetas numéricas 0, 1, . . , M , que se usarán en adelante para denotar los estados posibles del sistema. Así la representación matemática del sistema físico es la de un proceso estocástico {Xi}, en donde las variables aleatorias se observan en t = 0, 1, 2,. . ., y en donde cada variable aleatoria puede tomar el valor de cualquiera de los M + 1 enteros 0, 1, .. , M . Estos enteros son una caracterización de los M + 1 estados del proceso.(1)

Ejemplo.

Acorde a Kennedy, Snell y Thompson, la Tierra de Oz es bendecida por muchas cosas pero no por un buen clima. Ellos nunca tienes dos buenos días seguidos. Si ellos tienen un buen día, probablemente tendrán nieve o lluvia al siguiente. Si hay un cambio de lluvia o nieve, solamente la mitad de las veces es para cambiar a un buen día. Con esta información formamos una Cadena de Markov como la que sigue. Tomando como estados los diferentes tipos de clima R (lluvia), S (Snow) y N (normal). Con lo anterior se determina la matriz de transición de probabilidades.

Consideramos la mayoría de las veces la pregunta de determinar la probabilidad que, dada la cadena en un estado i hoy, esta llegue a un estado j dos días después. Normalmente esta probabilidad es denotada por



En el ejemplo planteado, vemos que si llueve hoy el evento que nevara dos días después es la unión separada de los siguientes tres eventos: 1) llueve mañana y neva dos días desde hoy; 2) es un buen día mañana y neva dos días después de hoy, y 3) neva hoy y también dos días después de hoy. La probabilidad del primero de estos eventos es el producto de la probabilidad condicional que llueva mañana, dado que llueve hoy, y la probabilidad condicional de que neve pasado mañana, dado que llueve hoy. Usando la matriz de transición P podemos escribir este producto como:

Los otros dos eventos también tienen probabilidades que pueden ser escritas como productos de las entradas de P. Así, tendremos que:

Para se plantea el siguiente teorema para probar lo observado anteriormente.

Teorema 1. Sea P la matriz de transición de una Cadena de Markov. La ij-ésimo entrada Pij^(n) de la matriz P^(n) da la probabilidad que la Cadena de Markov, comenzando en un estado inicial i esté en otro j después de n pasos.

Considerando de nuevo el ejemplo anterior de la Tierra de Oz, sabemos que el poder de la matriz de transición nos brinda información de interés acerca del proceso. Estaremos particularmente interesados en el estado de la cadena después de un largo número de pasos. Para ello es de gran utilidad utilizar cualquier programa que permita estos cálculos iterativos, entre ellos la hoja de cálculo EXCEL.

Tendremos por ejemplo después de 6 pasos lo siguiente:


Teorema 2. Sea P la matriz de transición de una Cadena de Markov y sea u el vector probabilidad el cual representa la distribución inicial de cada uno de los estados. Entonces la probabilidad de que la cadena esté en el estado i después de n pasos es la i-ésima en el vector: 

Además de ello, notamos que si queremos examinar el comportamiento de la cadena bajo el supuesto que comienza en un estado inicial i, simplemente se escoge siendo este el vector de probabilidad con la i-ésima  entrada igual a 1 y todas las otras igual a cero.
Ejemplo.

Siguiendo con el ejemplo de la Tierra de Oz, hacer el vector probabilidad inicial u igual (1/3,1/3,1/3). Entonces podremos calcular la distribución de los estados después de tres días usando el teorema 2 y nuestro cálculo previo de P^3, con lo cual obtenemos:



Visualicemos otro claro ejemplo de las aplicaciones de las cadenas de Markov:

Teniendo en cuenta el vector de Estados iniciales que describen el comportamiento de los abonados a las marcas de telefonía celular Movistar, Tigo y Comcel, y la matriz de transición dada, determine el estado estable para este sistema, modelando el problema como una cadena de Markov.

 
Donde:

M: Movistar
T: Tigo
C: Comcel
P0: Período inicial

Se establecen los estados en los diferentes cambios de período, para ello se utilizará la hoja de cálculo de Excel que facilitará el desarrollo del ejercicio por iteraciones.

Sabemos que el cambio de período se determina multiplicando el vector de estado con la matriz de transición así:

Tendremos entonces que:

Como podemos apreciar en Estado estable no hay cambias significativos en los resultados entre períodos, por lo cual siempre permanecerán en dicho estado.
Otra manera de conseguir dicho estable es apartir de resolución de ecuaciones por cualqueir método algebráico, entre ellos el de Gauss-Jordan. Esto consiste en establecer un sistema de ecuaciones que describa la situación problema. Tendremos a partir del vector de períodos como incógnitas y la matriz de transición dada que:
Multiplicando ambas matrices nos quedará el siguiente sistema:
Sumado a estas se consigna la sumatoria de las 3 incógnitas a partir de la teoría de la probabilidad, esto es:


Tendremos entonces un sistema de ecuaciones 4x3 que puede ser resuelto por cualquier método. Utilizando Excel tendremos que:


Como podemos apreciar se ha conseguido el estado estable directamente sin iterar. Este es otro método para hallar dicho estado.
Entre los conceptos más importantes tenemos lo que se conoce como estados Recurrentes que son aquellos que si después de haber entrado, el proceso definitivamente regresará a este. Por consiguiente un Estado es Recurrente si y sólo si no es transitorio. Además, ya que un estado Recurrente será visitado de nuevo de cada visita, podría ser visitado un número infinito de veces si el proceso continúa por siempre.

Si el proceso entra a cierto estado y permanece en este estado al siguiente paso, se considera un regreso a este estado. Entonces, el siguiente tipo de estado se considera un tipo especial de estado recuerrente, el Estado Absorbente.
Un estado se llama Absorbente, si después de haber estrado en él, el proceso nunca saldrá de este. Por consiguiente, el estado i es un estado absorbente si y sólo si Pij=1.
CADENAS DE MARKOV ABSORBENTE (ABSORBING MARKOV CHAIN)

Un estado inicial i de una cadena de Markov es llamado Absorbente si es imposible dejarlo (por ejemplo Pii=1). Una Cadena de Markov es absorbente si ésta tiene al menos un estado absorbente y si desde todos los estados es posible llegar a dicho estado (no necesariamente en un paso).

Ejemplo:

Almacenes Juanchi´s vende partes de automóviles y camiones a empresas que cuentan con flotas. Cuan las empresas compran a Juanchi´s se le dan 3 meses para pagar si las cuentas no se saldan en ese período, Juanchi´s cancela la cuenta, la remite a una agencia de cobranzas y da por terminadas las transacciones. Por lo tanto Juanchi´s clasifica sus cuentas en nuevas, de un mes de retraso, de dos meses de retraso, de tres meses de retraso, pagadas e incobrables. Juanchi´s investigó sus antiguos recursos y descubrió que:

a)      70% de las cuentas nuevas se pagan en un mes.

b)      60% de las cuentas de 1 mes de retraso se liquidan al final de mes.

c)       50% de las cuentas con 2 meses de retraso se pagan al final de ese último mes.

d)      60% con 3 meses de retraso se remiten a una agencia de cobranza.

1.  Forme la matriz de transición con estos datos.
2. ¿Cuál es la probabilidad de una que una cuenta nueva se liquide?
3. ¿Cuál es la probabilidad de que una cuenta con 1 mes de retraso se convierta en incobrable?

4. ¿En cuántos meses debe esperar Juanchi´s para que un cliente nuevo promedio liquidara su cuenta?
5. Si las ventas de Juanchi´s son en promedio US$125000/mes, ¿cuánto dinero se aceptaría como deuda incobrable al mes y cada año?

En primera instancia se plantea a continuación la matriz:


Se debe ahora considerar las matrices absorbente y no absorbe, así como la identidad


Además de ello, se debe tener en cuenta,


A partir de la anterior matriz podemos hallar la matriz fundamental para conseguir las probabilidades, para esto hallamos la matriz inversa de (N - I) por la matriz absorbente, obteniendo las probabilidades de ocurrencia de cada estado contenido en la matriz de no absorbente.


Ahora multiplicamos por la matriz absorbente así:


Con lo cual obtenemos las probabilidades para cada uno de los estados para llegar a cada uno de los estados absorbentes:


2. La probabilidad de que una cuenta nueva se liquide es 96%.
3. La probabilidad de que una cuenta de 1 mes de retraso se convierta en inconbrable es 12%.

4. El tiempo que se espera para que un cliente nuevo promedio liquidara su cuenta está dado por la suma de las probabilidades de la matriz (N - I) inversa, correspondiente a la primera fila. Esto es:


Sumando tendremos que este tiempo en meses es 1+ 0.3 + 0.12 + 0.06 = 1.48 meses.

5.  El dinero que se aceptará como deuda incobrable es de:



REFERENCIAS.
(1) Tutorial de Investigación de Operaciones II. Markov Chain. [En línea] [PDF] [Disponible en:


(2) GONZÁLEZ, Medardo. Clase Presencial. Investigación de Operaciones II . Universidad del Atlántico. Mayo de 2011.



REGLAS DE ORO EN LA GESTIÓN DE INVENTARIO

Cuando se hace parte del área de Inventarios de una empresa es importante tener en cuenta los siguientes aspectos:

1. Las tablas de Entrada y Salida de inventario deben estar dibidamente documentadas.

2. Todo Item (elemento que está guardado en inventario), debe estar debidamente codificado. Para esto se tiene en cuenta Ubicación (General) y Localización (Específica).

3. En cuanto sea posible todos los items deben estar guardados en el mismo lugar.

4. Nunca jamás recibir comisiones de algún proveedor.

5. En cuanto sea posible el lugar físico donde se efectúe la recepción de materiales debe ser diferente de donde se hace la salida.

6. Los items deben ser almacenados en orden con respecto a su peso, es decir, los más livianos arriba y los más pesados abajo.

7. Ningún miembro del equipo del almacén se puede ir hasta que no haya un conteo físico de los items que se movieron ese día.

8. Los items deben ser contados por 3 auditores diferentes y se consigna aquella lectura que sea coicidente entre al menos dos de ellos.

9. El punto más lejos donde debe haber un extintor en el almacen es máximo cada 10 metros.

10. Los reportes de inventario deben estar listos como máximo 3 días después del cierre.

PUNTO DE REORDEN CON DEMANDA PROBABILÌSTICA (MODELO EOQ)


El punto de Reorden r es un factor que sirve para garantizar que el de inventario manejado por la empresa en el almacen satisfaga el nivel de servicio al cliente. Para esto se debe tener en cuenta un Stock o inventario de seguridad  dado por un factor estocàstico definido como:
Donde sigma t representa la desviación estándar del producto en el tiempo de entrega o:


Cuando no se conoce la desviaión del tiempo de entrega, o no son iguales los períodos dados de la demanda y del tiempo de entrega, se resuelve diviendo ambos períodos como se mostró en la anterior ecuación.

Este ùltimo se usa para cubrir las variaciones en la demanda, ademàs de prever situaciones tales como: faltas del proveedor, mal manejo de materiales, etc. Además de ello representan niveles de inventario que cubran la variabilidad de la demanda, así como ajustar las demás variables de manera que se vuelvan irrelevantes. A pesar de ello, la mayor incertidumbre es causada por los clientes debido a que sus requerimientos son difíciles de predecir.

La fórmula de Punto de reorden se muestra a continuación:

Ejemplo:

A partir de los siguientes parámetros estocástico de la demanda de focos, una empresa decide definir un Punto de Reorden adecuado para cumplir con los requerimientos de sus clientes.

La demanda en 1 semana presenta una distribución gaussiana de media 154 focos/semana y desviación estándar de 25 Focos/semana.

En primera instancia calculamos el EOQ con la demanda semanal promedio 154 focos. Entonces tendremos que:


No obstante, hay que tener en cuenta que estas unidades sólo suplen la media de la demanda, es decir el percentil de 50%. Pero ello quiere decir que no tendremos existencias para suplir aquellas que superen el valor medio.

El director de Operaciones de la empresa sugiere la política de permitir que sólo fallen en 1 pedido. Entonces debemos determinar el número de pedidos ideal. Esto es:


Ahora, considerando la política optada tendremos que el nivel de significancia será:

A partir de lo anterior tendremos que el Punto de

La empresa debe encargar a su proveedor cada vez que haya sólo 195 unidades en almacen.



REFERENCIAS

(1) HILLIER, Frederick S., LIEBERMAN, Gerald J. Investigación de Operaciones. Edición Séptima: McGraw-Hill, 2002. Teoría de Inventarios: Páginas 935-987.

(2) TAHA, Hamdy A. Investigación de Operaciones. Edición Séptima: Prentice Hall, 2004. Modelos Determinísticos de Inventario. Páginas 429-440.

(3) GONZÁLEZ, Medardo. Clases Presenciales: Modelos de Inventario. Investigación de Operaciones II. Universidad del Atlántico. Febrero de 2011.










20 de marzo de 2011

Modelo de EOQ con Descuento por Cantidades

Hasta este momento se han analizado los costos en que se incurre por pedir o por producir con una demanda constante, ya sea permitiendo o no faltantes en nuestros parámetros de inventario; asumiendo además que para el caso de producir la rata de producción debe ser mayor que la demanda del período. Esto dejando a un lado el efecto del costo unitario que se mantiene fijo en el tiempo y no varía si hay alguna modificación en la cantidad a pedir o producir. Dicho efecto no debe dejarse a un lado, ya que muchas empresas trabajan bajo este criterio, siempre, cuando se varía el nivel de pedido habrá un descuento adicional sobre el costo del mismo. Es por esto que a continuación se hablará acerca del Modelo EOQ con Descuento por Cantidades. Para ello se palnteará la siguiente situación.

Suponga que una empresa informática se dedica a la venta de computadores y actualmente no cuenta con un modelo de inventario óptimo. En estos momentos trata de determinar cómo minimizar los costos anuales relacionados con la compra de tarjetas de video para los computadores. Cada vez que se hace un pedido se incurre en un costo de $20. El precio de las taretas depende del número pedido de ellas. Si el costo anual de almacenamiento es de 20% del valor de inventario y la empresa demanda 80 tarjetas al mes, ¿cuál es la cantidad óptima de pedido? ¿cuál es el Costo Total Anual mínimo? ¿Cuál es la cantidad óptima anual? y ¿cuál es el período por pedido según las cantidades anuales?

Este tipo de situación es común cuando se compra al por mayor, para esto se debe analizar según los parámetros del modelo EOQ. 

La siguiente tabla muestra los descuentos según el volumen pedido:

A partir de lo anterior se analizan las cantidades óptimas por pedido (se asume que no se aceptan faltantes) para cada descuento, siempre y cuando se tenga en cuenta lo siguiente:

- La cantidad óptima de pedido será aquella que, al momento de calcular las cantidades a partir del modelo adecuado y según los costos unitarios de cada descuento, este valor se halle dentro del rango de cantidades especificado (en la tabla, la primera columna). Caso contrario, si dicho valor no se encuentra en su rango especificado, se dice que No Cumple para el período, y la cantidad óptima será el valor más cercano a su límite de especificación.

Los datos conocidos son:

Cp = $20/pedido
Cmi = 20% Cu
d = 80 tarj/mes * 12meses/1año = 960 tajertas/año

Así, se procede en primera instancia a calcular las cantidades óptimas para cada rango especificado en la tabla, esto es:



Los comentarios en paréntesis están relacionados al criterio antes señalado, y nos permite concluir que la cantidad óptima de pedido en un período es de 139 tarjetas de video.

No obstante, no siempre la cantidad óptima de pedido para un período representa la cantidad que minimiza los Costos Totales Anuales. Recordemos que la cantidad óptima de pedido sólo permite que los costos de pedir y de mantener en inventario sean iguales, y si el costo unitario es constante, esta llega a ser la cantidad adecuada para optimizar. Pero, en esta situación, la variación de los costos unitarios generan un segundo enfoque, ya que dependerá principalmente de la relación cantidad-precio unitario. Con base en ello debe resolver la función CTA para las tres cantidades óptimas optenidas anteriormente. Esto es:


Así tendremos que:

A partir de los resultados obtenidos, podrmos apreciar que el Costo Total Anual mínimo es el incurrido cuando se toman pedidos de 300 tarjetas y es de $9766,00.

Por otro lado, analizando el comportamiento de los costos nos podemos dar cuenta que:

- A medida que se pide una cantidad mayor la diferencia entre los costos de pedir y de mantener en inventario es cada vez más grande, siendo los costos de mantener por obvias razones la que aumente y de pedir el que disminuye.

Finalmente, conociendo la Cantidad óptima anual, podemos determinar el período para cuando son más de un pedido. Así tendremos que:



Como se pudo apreciar, el modelo de EOQ con Descuentos por Cantidades representa una herramienta útil cuando existe un porcentaje de reducción en el costo unitario de los artículos que se ordenan. Análogamente al EOQ con y sin faltantes, se manejan los mismos parámetros de evaluación, la diferencia radica en que se debe hacer un análisis para cada uno de los precios unitarios especificados por el proveedor tanto para establecer la cantidad óptima en un pedido, como para establecer la cantidad óptima anual y los Costos Totales Anuales mínimos, siendo los últimos los más representativos para cualquier organización.



REFERENCIAS

(1) HILLIER, Frederick S., LIEBERMAN, Gerald J. Investigación de Operaciones. Edición Séptima: McGraw-Hill, 2002. Teoría de Inventarios: Páginas 935-987.

(2) TAHA, Hamdy A. Investigación de Operaciones. Edición Séptima: Prentice Hall, 2004. Modelos Determinísticos de Inventario. Páginas 429-440.

(3) GONZÁLEZ, Medardo. Clases Presenciales: Modelos de Inventario. Investigación de Operaciones II. Universidad del Atlántico. Febrero de 2011.




Modelo de Lote Económico de Producción (LEP) Con Faltantes

Para este caso se manejan los mismo parámetros establecidos en el modelo de LEP sin faltantes, entre ellos se debe cumplir que la rata de producción (R) sea mayor que la demanda (D), además de incluir un costo por los artículos o unidades faltantes en los costos totales. Además de ello, se debe considerar dos nuevos períodos que son el perído desde que se agota el inventario y empieza haber faltantes y otro que se considera como el tiempo en que se nivelan dichos faltantes. A continuación se muestra el gráfico que ilustra lo dicho anteriormente:


Donde:

R: rata o radio de producción
d: Demanda del período
t1: Período de producción hasta el inventario máximo
t2: Período en que la demanda consume el inventario máximo
t3: Período en se comienza a incurrir en pedidos en espera
t4: Período en que se sustentan todos los pedidos atrasados del período anterior

Teniendo en cuenta esto, nuestra función de costos para un período es:



Siendo Cf  el costo por faltantes.

A partir de la gráfica antes señalada se obtienen por relación de triángulos las siguientes ecuaciones



Teniendo en cuenta esto, se deben expresar los Costos totales en función de las variables conocidas, es decir, Q y S. Utilizando este sistema de ecuaciones tendremos que:

Despejando t3 y t4 de (3) y (4), y posteriormente sumándolas tendremos que:



Ahora, si se considera que:


Entonces despejando la suma de t1 más t2, nos quedará que:


Siendo t1 igual a:


Con este lograremos deducir  Imáx:


Con base en esto podemos deducir la siguiente Función de Costos por período:

 De la misma forma como en los modelos anteriores proseguimos con la determinación del Costo Total anual multiplicando la expresión por N, esto es:


Lo importante sería optimizar la función para determinar cuales serían nuestas cantidades óptimas a producir y cual sería el nivel tanto de inventarios como de faltantes que nos proporcionarían los Costos mínimos en un período dado. Con base en este análisis, y teniendo en cuenta que la función de Costos depende de Q y S. Para ello primero consideramos la derivada parcial con respecto a S así:


Asimismo, igualando a cero la expresión y despejando a Q tendremos que:

 Se procede con la derivada parcial con respecto a Q y esto es:

Igualando a cero y despejando S cuadrado, esto es:



Ahora, reemplazando en la expresión el Q hallado anteriormente y simplificando la expresión, nos quedará lo siguiente:



Con lo cual nos queda que nuestro S óptimo será:


Reemplazando esta expresión en el Q hallado previamente nos permitirá encontrar nuestra candidad óptima por corrida de producción, esta es:



A partir de estos también podremos deducir nuestro Período óptimo en que incurre este modelo, así como el Inventario máximo:




REFERENCIAS

(1) HILLIER, Frederick S., LIEBERMAN, Gerald J. Investigación de Operaciones. Edición Séptima: McGraw-Hill, 2002. Teoría de Inventarios: Páginas 935-987.

(2) TAHA, Hamdy A. Investigación de Operaciones. Edición Séptima: Prentice Hall, 2004. Modelos Determinísticos de Inventario. Páginas 429-440.

(3) GONZÁLEZ, Medardo. Clases Presenciales: Modelos de Inventario. Investigación de Operaciones II. Universidad del Atlántico. Febrero de 2011.