Package BysSampler

El paquete BysDecision ofrece utilidades relacionadas con problemas de decisión bayesiana en ambiente de incertidumbre.

Definiciones

Sea el espacio de todas las decisiones o acciones posibles que puede tomar un decisor.

Sea una variable aleatoria vectorial con dominio cuya distribución conjunta es conocida y que recoge todas las situaciones de las cuales no es posible tener un conocimiento determinista antes de tomar la decisión.

Para cada posible decisión y cada posible situación , es decir cada posible realización de , el coste de la acción para el decisor es una función real conocida

Se pretende encontrar la decisión que minimiza la esperanza del coste

Nótese que en lugar de función de coste podría definirse una función de beneficio a maximizar sin más que multiplicarla por menos uno.

Decisión binaria

Cuando sólo hay dos posibles acciones, hablaremos de decisión binaria o booleana. Por ejemplo, un agricultor debe decidir si ir a regar un campo a manta o no hacerlo, pero no puede decidir qué cantidad de agua utilizar sino que ésta está fijada por el método de riego. El coste de la acción depende de que llueva o no y en qué medida lo haga:

• Si decide regar y
  • no llueve de forma significativa el coste de la acción es nulo,
  • llueve de forma moderada pierde el tiempo y el combustible necesarios para desplazarse,
  • llueve de forma excesiva pierde además parte de la cosecha de forma proporcional al exceso de agua.
• Si decide no regar y
  • no llueve de forma significativa pierde toda la cosecha
  • llueve de forma moderada el coste de la acción es nulo,
  • llueve de forma excesiva pierde parte de la cosecha de forma proporcional al exceso de agua.

La importancia de este tipo de problemas reside en que cualquier problema de decisión en el que haya un número finito de opciones cualitativas puede reformularse como un árbol de decisiones binarias.

Decisión continua

Si el conjunto de decisiones tiene medida no nula en diremos que se trata de una decisión continua. Esto no implica necesariamente que la variable aleatoria sea continua.

Por conveniencia, también se aplicará esta etiqueta cuando el espacio de opciones sea discreto y ordenado y el número de opciones sea suficientemente grande para despreciar el efecto de la discretización. En tales casos la solución continua se discretizará por redondeo directo o probando las soluciones discretas más cercanas a la continua.

Por ejemplo, en las apuestas deportivas se puede decidir apostar cualquier cantidad de dinero entre 0 y el máximo admitido por la casa de apuestas, mientras que el conjunto de resultados posibles suele ser finito y muy pequeño: {gana local,gana visitante}, {1,X,2}, etc.

Decisión biunívoca

Cuando cada posible decisión se puede asociar de forma biunívoca a cada posible situación, es decir, , entonces diremos que se trata de decisión biunívoca. Un caso típico es la distribución de material perecedero para su venta entre un conjunto de puntos de venta, en el que la demanda en cada punto de venta sería la variable aleatoria. Si se pone en servicio demasiado material en un punto dado, el exceso se perderá total o casi totalmente, si es posible reciclarlo en alguna medida, mientras que si se sirve poco se pierde el beneficio de la venta no realizada y se aumenta el riesgo de que la clientela se pase a la competencia. Tanto las decisiones posibles como las situaciones inciertas son en este caso

Cuando el dominio es continuo, este tipo de problemas pueden resolverse mediante técnicas de optimización muy eficaces y robustas , por lo que a veces es interesante convertir el problema al caso continuo.

Decisión con restricciones

Supongamos que en el caso anterior existen restricciones, como por ejemplo que la suma de las cantidades distribuidas no pueda superar las existencias disponibles.

Ya no se puede hablar propiamente de decisión biunívoca pues el conjunto de decisiones es un subconjunto del de situaciones, pero también están disponibles las mismas o gran parte de las técnicas de optimización aplicables a esta clase de problemas añadiendo las restricciones al problema de optimización.

Combinación de previsiones

Un caso particular de decisión con restricciones, que resulta muy interesante para los analistas de datos estadísticos, es la mal llamada combinación de previsiones, refiriéndose al problema de obtener estadísticos de un conjunto de variables aletarorias que sean congruentes con ciertas restricciones, las cuales deben cumplir en realidad las propias variables aleatorias, y por ende cada una de sus realizaciones conjuntas, pero no necesariamente sus estadísticos centrales: media, mediana y moda.

Si las distribuciones son todas normales y las restricciones son todas lineales, la media cumplirá dichas restricciones, pero en cualquier otro caso esto no ocurrirá así, y no siempre se tiene presente este hecho.

Un ejemplo típico de esta situación es la presencia de una variable aleatoria (a la que llamaremos global) que es la suma de las demás y para la que existe un modelo log-normal(o cualquier otra no normal) con menos varianza que en el de las partes también log-normales. Si el modelo se estima por simulación y conjuntamente, obligando a que cada muestra cumpla las restricciones, podremos ver cómo las medias de las partes no suman exactamente la media del global. Si cada modelo se estima por su lado ni las realizaciones ni los estadísticos cumplirán las restricciones. En este caso sería posible combinar las variables originales para dar otras que cumplieran las restricciones, pero como ya se ha dicho, los estadísticos centrales no tendrían porqué cumplirlas.

Como en muchas ocasiones se confunde el concepto de previsión, que es una variable aleatoria, con alguno de sus estadísticos centrales, que son simples números, los clientes suelen exigir que esos números sean coherentes con ciertas restricciones que son incompatibles con seguir siendo dichos estadísticos.

El analista debe tomar una decisión que consiste en decirle al cliente que la previsión es cierto número aún a sabiendas de que la previsión es una variable aleatoria y que lo que le entrega no es una previsión ni un estadístico de la msima, sino una decisión que es posible optimizar.

Lo ideal en estos casos sería convencer al cliente para investigar cuál es su función de coste real y calcularla o aproximarla en la medida de lo posible, o bien inferirla directamente si el negocio es suficientemente conocido y transparente.

Si el cliente es refractario a estos conceptos y no hay forma de averiguarla sin su colaboración, el analista sería muy libre de buscar su propia función de coste, entendido como el grado de rechazo por parte del cliente de las "previsiones" entregadas, o mejor aún, estableciendo en el contrato de servicio el precio según cada tipo de error. Si se observa que al cliente le gusta que el sesgo sea simétrico pues se pone una función simétrica, etc.; si prefiere equivocarse un poco más al alza que a la baja pues se va probando con distintos coeficientes relativos; etc. En cierta manera esto sería una forma de sonsacarle al cliente la función de coste de una forma inconsciente, pues muchas veces él tiene la intuición de lo que es bueno o malo para su negocio, aunque no sepa o no quiera trasladarlo a una formulación matemática.

En ausencia de información, el último recurso sería considerar como función de beneficio el logaritmo de la función de verosimilitud, o su opuesta como función de coste. En el caso de residuos normales estandarizadas e independientes la función de coste sería simplemente la suma de cuadrados de los los sesgos. En el ejemplo explicado anteriormente de las log-normales cuya suma debe dar otra log-normal, habría que transformar las variables mediante el logaritmo y premultiplicarlas por una descomposición simétrica de la matriz de covarianzas, para obtener residuos normales estandarizados e independientes, y sustituir en la ecuación de restricción la variable original por su expresión en términos de ellos. Esta situación es generalizable a cualquier tipo de restricciones de igualdad y/o desigualdad

sobre variables aleatorias transformadas de una multinormal mediante una función inversible

Los residuos normales estandarizados e independientes serían

La decisión óptima en la métrica transformada sería por tanto la solución del programa no lineal

Sujeto a:

Para obtener la solución en términos originales bastaría con calcular