1. Introducción

Las empresas están expuestas a riesgos que se derivan de las fluctuaciones de los mercados económicos y financieros, ya que influyen directamente en los precios de activos financieros. El cambio generado en estos, puede conducir a ganancias o pérdidas increíblemente altas (Oro, 2004). Una propuesta para la administración de tales riesgos ha sido la creación y el desarrollo de instrumentos derivados. Un derivado es un instrumento financiero cuyo valor depende de otro activo llamado subyacente, su usopermite, entre otros, asegurar precios futuros y cubrirse contra potenciales pérdidas.

Las opciones plain vanilla son derivados que han alcanzado un gran desarrollo en los mercados estandarizados. La variedad de subyacentes sobre los que se suscribe este tipo de contratos es extensa, a saber, divisas, acciones, índices bursátiles, tasa deinterés, commodities, seguros y energía. Las opciones suelen clasificarse de distintas formas, una de ellas es por el derecho que dan al poseedor del contrato: si el contrato confiere al poseedor de la opción el derecho a comprar el activo subyacente,será una opción de compra también llamada call, pero si el contrato le confiere el derecho a vender el activo subyacente, será una opción de venta también llamada put. Otra clasificación es de acuerdo a la fecha de ejercicio del contrato, es europea si el contrato se ejerce justo al vencimiento de este y americana cuando se puede ejercer en cualquier momento durante la vida del contrato.

Sin embargo, las necesidades específicas de los inversionistas, dieron origen a la creación de las opciones de segunda generación, también conocidas como opciones exóticas, las cuales se negocian principalmente en mercados Over The Counter (OTC). Lasopciones asiáticas pertenecen al grupo de las llamadas opciones exóticas, su nombre se debe a que se comercializaron por primera vez en Tokio en 1987 (Venezia, 2010). La principal característica de este tipo de opciones es que para calcular su prima se toma en consideración el promedio de la evolución de precios del activo subyacente durante la vida del contrato, haciendo que su valuación sea menos vulnerable a la manipulación de precios; lo que las hace bastante competitivas como instrumento de especulación y arbitraje, y motiva su uso en portafolios de inversión. Por ejemplo Arregui y Vallejo (2001) muestran la relevancia de la valuación de opciones asiáticas que forman parte de un portafolio de inversión para calcular la rentabilidad que se puede garantizar además de valorar la situación patrimonial de un fondo de inversión garantizado de renta variable en cualquier momento.

Las opciones asiáticas además de ser de compra o venta, europea o americana, tienen su propia clasificación y esta se da por el tipo de promedio, aritmético o geométrico. Para un análisis de la valuación de opciones asiáticas con promedio geométrico pueden leerse los trabajos de Wang, Wang, Wang, y Jin, (2014), Zhang y Liu (2014), Kim y Wee (2014), Prakasa (2016), Jeon, Yoon, y Kang (2016), Zhang, Xiao, Kong y Zhang (2015) y Hubalek, Keller-Ressel y Sgarra (2017). Mientras que para opciones asiáticas con promedio aritmético refiérase a Sengupta (2014), Chung y Wong (2014), Mehrdoust y Saber (2015), Zeng y Kwok (2016), Hozman, Tichý y Cvejnová (2016) y Huang, O’Hara y Mataramvura (2017), por mencionar solo algunos. Es necesario precisar, que si el precio de ejercicio de una opción asiática depende de una cantidad fija, la opción se conoce como una opción asiática con precio de ejercicio constante. Pero si el precio de ejercicio es proporcional al precio del activo subyacente, entonces es una opción asiática con precio de ejercicio flotante o promedio.

En lo que concierne a métodos de valuación, se encuentran disponibles en la literatura métodos tanto numéricos como analíticos. Respecto a métodos numéricos, pueden leerse los trabajos de: Yor (1992), en donde expresa el valor de la opción como una integral triple difícil de evaluar numéricamente, Dufresne (2000) quien obtiene una fórmula que expresa el valor de la opción asiática como una serie infinita de polinomios de Laguerre con cada coeficiente en la serie dada por una única integral que necesita ser calculada numéricamente, la limitación que encontró es que solo aplica a promedios continuos con ponderaciones iguales además de que no hay límite en el error que se produce al truncar la serie. Linetsky (2004) desarrolla dos representaciones, la primera es una serie infinita de términos relacionados con las funciones de Whittaker M y W. y la segunda es una integral real única de una expresión que involucra la función Whittaker W. plus y funciones gamma incompletas, así como polinomios de Laguerre encontrando una precisión de 6 decimales al comparar con Vecer (2002). Autores como Turnbull y Wakeman (1991) encuentran la aproximación lognormal con los primeros y segundos momentos combinados obteniendo el precio de una opción conmedia aritmética.

Particularmente, el método de simulación de Monte Carlo es utilizado frecuentemente en Finanzas, porque permite encontrar soluciones aproximadas de problemas matemáticos que involucran variables aleatorias dependientes del tiempo (Venegas-Martínez, 2008), al respecto se puede consultar a: Boyle (1977), quien usa el método de simulación de Monte Carlo para resolver problemas de valoración de opciones, simula el proceso que generan los rendimientos en el activo subyacente suponiendo neutralidad de riesgo para derivar el valor de la opción e introduce una variable de control apropiada con la cual se aumenta significativamente la precisión de los resultados, Kemna y Vorst (1990) ofrecen un método de precios mejorado incluyendo la simulación de Monte Carlo con elementos de reducción de varianza. Por su parte Levy (1992), proporciona aproximaciones analíticas de forma cerrada para valuar opciones europeas que involucran el promedio aritmético de las tasas de cambio futuras. Asimismo se han desarrollado investigaciones cuya característica es utilizar el método de Monte Carlo con variables de control, por ejemplo, Fu, Madan y Wang (1999) utilizan variables de control para corregir el sesgo de discretización inherente a la simulación, también pueden consultarse: Datey, Gauthier y Simonato(2003), Ben y Kebaier (2014), Mehrdoust (2015), Shiraya y Takahashi (2017), entre otros.

En lo que a métodos analíticos compete, el enfoque de ecuaciones diferenciales parciales ha sido recurrentemente utilizado para valuar instrumentos derivados, particularmente opciones. Véanse los trabajos de Zvan et al. (1997), Rogers y Shi (1995) quienes realizan la valuación de la opción asiática resolviendo una ecuación diferencial parcial parabólica en dos variables además de proporcionar un límite inferior que es tan preciso que es esencialmente el precio real. Por otra parte, Vecer (2001 y 2002), Marcozzi (2003), Foufas y Larson (2008) obtienen el precio de la opción asiática con el modelo de Black-Scholes, las ecuaciones diferenciales parciales resultantes son de tipo parabólico en una dimensión espacial y el cálculo numérico se realiza mediante el uso de un método de elementos finitos adaptativos.

De particular mención para esta investigación son los modelos que utilizan métodos mixtos, en la búsqueda de valuación de derivados en escenarios más apegados a la realidad. Un instrumento muy valioso en este rubro es la Teoría de Control Óptimo Estocástico en Optimización Dinámica en tiempo continuo, en este contexto existen modelos Microeconómicos o Macroeconómicos disponibles en la literatura, en los que realizan la valuación de activos financieros y/o instrumentos derivados. Véanse por ejemplo Sierra (2007), Cruz-Aké y Venegas-Martínez (2010), Martínez-Palacios, Sánchez-Daza y Venegas-Martínez (2012), Cai, Song y Kou (2015), Venegas-Martínez y Aguilar-Sánchez (2005) y Martínez-Palacios, Venegas-Martínez y Martínez-Sánchez (2015), quienes modelan en tiempo continuo la toma de decisiones óptima de consumo y una cartera de activos en un horizonte temporal finito, de longitud estocástica, y en el proceso de solución, obtienen el modelo equivalente a Black-Scholes-Merton para valuar una opción de venta de tipo americano; Sierra (2007) mediante un modelo de agente representativo y suponiendo la distribución de activos modelados mediante Brownianos fraccionales, obtiene el precio de una opción europea de compra. Por su parte Cai, Song y Kou (2015) propusieron fórmulas de aproximación de doble transformación de forma cerrada para los precios de las opciones asiáticas controladas de forma discreta y continúa estableciendo una cadena general de Markov de tiempo continuo.

En la presente investigación el marco de la Teoría de Control Óptimo Estocástico en Optimización Dinámica en tiempo continuo, y en un ambiente de riesgos de mercado, se desarrolla un modelo de un agente económico racional que dispone de una riquezainicial y enfrenta la decisión de cómo distribuir su riqueza entre consumo e inversión en portafolio de activos, uno de ellos una opción asiática de venta suscrita sobre una acción que tiene volatilidad estocástica, de forma tal que maximice su utilidad total esperada por el consumo e inversión, en un horizonte de planeación temporal finito. En el proceso de solución, mediante un sistema de ecuaciones diferenciales se deduce el modelo para valuar la opción asiática de venta de tipo europeo que tiene precio de ejercicio variable igual a la media aritmética. Para dar solución al modelo, y con el objetivo de un análisis profundo de solución, se establece una metodología mediante simulación Montecarlo en la que se supone que las distribuciones de los precios del subyacente y la volatilidad son conducidos por un movimiento geométrico browniano, la volatilidad del subyacente y sus parámetros se calibran mediante una adaptación del modelo de tasa corta de Cox, Ingersoll y Ross (1985) (CIR) a volatilidad realizada, la tasa de interés libre de riesgo se obtiene de la página de Banxico para los cuatro plazos de las opciones del boletín que se publica en la página web del Mercado Mexicano de Derivados (MexDer).

En el contexto de las investigaciones mencionadas, las características distintivas del presente trabajo son: 1) Se obtiene una ecuación diferencial parcial equivalente a la ecuación obtenida por Black-Scholes-Merton introduciendo el supuesto de volatilidad estocástica conducida por un movimiento geométrico browniano, 2) La valuación de la opción asiática de venta con promedio aritmético y europea se establece mediante micro fundamentos y 3) Se desarrolla una metodología mediante Simulación Monte Carlo para la valuación de opciones asiáticas de compra y de venta con el supuesto de que la volatilidad es estocástica y conducida por una ecuación diferencial estocástica con reversión a la media, con parámetros calibrados a partir de información de mercado.

Este trabajo está organizado de la siguiente forma, en la siguiente sección se definen los activos de inversión a los que tiene acceso el agente y se establece su representación analítica en el modelo, en la sección 3 se establece la ecuación de riqueza del agente representativo de esta economía, así como el modelo de Control Óptimo Estocástico por resolver. En el transcurso de la sección 4 se proporciona el proceso de solución al problema planteado utilizando programación dinámica estocástica en tiempo continuo y en la cual se obtiene los resultados centrales: i) ecuación diferencial parcial de Hamilton-Jacobi-Bellman y ii) las condiciones de primer orden, que conducen a las trayectorias óptimas. En la sección 5 se propone una funciónHARA y se da la solución del problema planteado, obteniendo las proporciones óptimas de consumo e inversión. En la sección 6, se deduce la ecuación diferencial parcial lineal de segundo orden equivalente a la de Black-Sholes-Merton, la cual es el modelo mediante el que se valuara una opción asiática de venta del tipo europeo con precio de ejercicio variable media aritmética suscrita sobre una acción cuya volatilidad. En la sección 7 se aproximan los precios de opciones europeas y asiáticas de compra y venta por simulación Monte Carlo con parámetros calibrados adaptando el modelo de Cox-Ingersoll-Ross con volatilidad realizada y su correspondiente análisis de resultados. En la sección 8 se presentan las conclusiones de la presente investigación indicando aquellas extensiones previstas en la agenda futura de trabajo. Por último, se presentan las referencias bibliográficas utilizadas y un Apéndice en el que se dispone en un cuadro los precios obtenidos de opciones europeas de compra y de venta sobre WALMEXV.

2. Descripción del marco teórico

Se supone una economía pequeña y cerrada, en la que se comercializa un solo bien de consumo de carácter perecedero. En esta economía existe un sistema bancario o mercado de crédito en el que son posibles las negociaciones de forma continua, es decir,se opera los 365 días del año, las 24 horas del día. En este sistema, las ventas en corto son permitidas e ilimitadas.

3. Ecuación de riqueza y planteamiento del Modelo

4. Proceso de solución del problema: ecuación diferencial parcial de Hamilton-Jacobi-Bellman y condiciones de primer orden

5. Función de utilidad y proporciones óptimas de inversión

6. Ecuación Diferencial Parcial para valuar una opción asiática de compra

Al suponer la solución de esquina con α 1 t = 1 y α 2 t = 0 y sustituir estos valores en los premios al riesgo de la acción y la opción, se obtiene la fórmula de valuación para una opción asiática de tipo europeo análoga a la EDP de Black-Scholes-Merton, a la que se debe de imponer la condición de frontera, que corresponde al valor intrínseco de la opción:

7. Calibración del modelo con volatilidad estocástica realizada y Simulación Montecarlo

En esta sección se describe un sistema de ecuaciones diferenciales estocásticas que conducen la dinámica de los rendimientos y la volatilidad del activo subyacente, además de una manera de calibrar sus parámetros con datos históricos y de mercado. Se supone que el precio del subyacente y su volatilidad son conducidos movimientos geométricos brownianos con o sin correlación, con la característica de que la ecuación que conduce la volatilidad es con un proceso de tipo CIR, posteriormente se simulan trayectorias del subyacente con la función de pago de la opción y se descuenta con la tasa libre de riesgo el valor final del promedio de las trayectorias. De esta manera se determina el precio de una opción asiática de tipo europea. La simulación de Monte Carlo es muy útil, en especial cuando no se tiene una fórmula cerrada para la valuación.

Una manera de considerar la naturaleza estocástica en los rendimientos de los activos financieros es modelarla por medio de procesos de reversión a la media, como por ejemplo con modelos de volatilidad estocástica de Hull y White (1987) y Heston (1993), entre otros. La calibración de los parámetros de estos modelos puede hacerse por series de tiempo, modelos con funciones de pérdida, diversos modelos de la familia GARCH, y con volatilidad realizada, entre otros.

A continuación, se describe la metodología para determinar el precio de una opción asiática con subyacente promedio sin pago de dividendos mediante simulación Monte Carlo.

Suponga que el precio del subyacente es conducido por un movimiento geométrico browniano con volatilidad estocástica σ t conducida por un proceso de reversión a la media y ambos procesos están correlacionados como sigue:

donde

W σ , t es un proceso de Wiener correlacionado con W S , t , es decir, Cov ( d W σ , t , d W S , t ) = ρ d t .

Si λ = 0 , entonces la dinámica de la volatilidad es conducida por un modelo de tipo Vasicek y se considera λ = 0.5 para un proceso de tipo CIR.

Los parámetros θ, κ y σ σ , se interpretan como la tasa de largo plazo, la velocidad de reversión a la tasa de largo plazo y la volatilidad de la varianza de la tasa de interés (denominada también como la volatilidad de la volatilidad), respectivamente.

Para ejecutar la simulación de ambos procesos se generan trayectorias con una estructura de relación dada por:

con:

Para hacer lo anterior, se calcula dado que ∑ = L L simular d Z ~ N ( 0 , I 2 ) para obtener d W ˜ = L d Z , se elige un número para la partición del plazo al vencimiento de la opción, por ejemplo N = 100 y mediante el método de remuestreo cuadrático (Barraquand, 1995) generar d Z , que involucra a d W ˜ por construcción. Sean μ Z y ∑ Z , la media teórica y la matriz de covarianza de d Z respectivamente:

Con las ecuaciones anteriores se generan las trayectorias del sistema de ecuaciones diferenciales estocásticas planteado en (24).

El precio de una opción de compra con subyacente promedio está dado por la función de pago:

y el precio de una opción de venta con subyacente promedio está dado por:

Si el subyacente es conducido por el sistema dado en (24), entonces el algoritmo para determinar los precios de las opciones de compra y de venta es:

1. 1. Generar d W S , t y d W σ , t como:

1. 2. Discretizar el sistema de ecuaciones diferenciales estocásticas como:

1. donde hace referencia a la k-ésima simulación.

1. 3. Definir la media aritmética de las trayectorias generadas:

1. 4. Calcular el precio de la opción como:

1. donde denota el número de trayectorias simuladas y el número de precios generados.

7.1 Análisis y discusión de resultados

En esta sección se calculan por medio de simulación Monte Carlo precios de opciones asiáticas con subyacente promedio, el subyacente es el precio de la acción de WALMEXV y se comparan con los precios de opciones publicadas en el boletín de MexDer deldía 16/05/2018. Se supone que los parámetros de volatilidad estocástica se calculan mediante una adaptación del modelo de CIR al método de volatilidad realizada. La tasa de interés se obtiene de la página de Banxico para los cuatro plazos de las opciones del boletín.

La muestra de la volatilidad realizada para calibrar los modelos comprende del 30 de diciembre de 2016 al 16 de mayo de 2018. Los parámetros estimados con el modelo de CIR mediante máxima verosimilitud se observan en el cuadro siguiente:

Cuadro 1.

Cuadro resumen de estimación de parámetros de volatilidad realizada con el modelo CIR por el método de máxima verosimilitud.

Parámetros	CIR
κ	2.280414
θ	0.013979
σ	0.002504
No. Obs. n	344
Likelihood Ratio	9.433697

^TFN1Fuente: Elaboración propia.

De los resultados del cuadro anterior, se puede verificar que se cumple la condición de Feller, es decir, las trayectorias del proceso que conduce la volatilidad de la volatilidad son siempre positivas. Los resultados de la simulación se muestran en la gráfica 1, en la cual se muestran la serie original de la volatilidad realizada y cien trayectorias simuladas con los parámetros dados por el modelo CIR. Se observa la tendencia de la volatilidad hacia la baja, explicado por la misma dinámica del precio de la acción de WALMEXV.

[Figure ID: f1] Gráfica 1.

Trayectorias simuladas.

  —Fuente: Elaboración propia..

[Figure ID: f2] Gráfica 2.

Comparación entre precios de opciones de MexDer, opciones europeas y asiáticas con parámetros calibrados con el modelo de CIR. Fuente: Elaboración propia.

La Gráfica 2 muestra los precios de opciones de compra y de venta europeas, y precios de opciones con subyacente promedio ambas con volatilidad estocástica con parámetros calibrados con volatilidad realizada conducida por el modelo de CIR, en ambos casos los plazos son desde T=30, 128, 219 y 303 días y los precios de ejercicio varían desde $43 hasta $57 con incrementos de $1, tales precios se comparan con los publicados por MexDer el día 16/05/2018. El precio de WALMEXV=$49.39 y número de trayectorias que se simularon para determinar los precios de las opciones fue de 10,000.

Se tienen quince precios de ejercicio de los cuales siete están dentro del dinero para opciones de compra y ocho están dentro del dinero para opciones de venta. Se observa que para opciones al plazo de T=30 días en los primeros tres precios de ejercicio, los precios de opciones de compra con subyacente promedio resultaron mayores que los precios publicados por MexDer. Para los plazos de T=128 días y T=219 días solo para el primer precio de ejercicio los precios de opciones de compra son menores en comparación con los precios de las publicadas por MexDer, mientras que para el plazo de T=303 días solamente los dos primeros precios de ejercicio. En cuanto a las opciones de venta el análisis del cuadro A.1 y la gráfica 2 resulta que los precios de opciones de venta con subyacente promedio son menores para los cuatro plazos considerados.

8. Conclusiones

En esta investigación se planteó y resolvió un problema de Control Óptimo Estocástico en Optimización Dinámica que modela la toma de decisiones de un consumidor racional quien dispone de una riqueza inicial y enfrenta la decisión de cómo distribuir su riqueza entre consumo y un portafolio de activos en horizonte de planeación finito, de tal forma que maximice su utilidad total esperada por el consumo e inversión. Entre los supuestos del modelo están que el agente puede invertir en una acción con volatilidad estocástica y una opción asiática de venta con precio de ejercicio media aritmética y de tipo europea, cuyas distribuciones respectivas de precios y volatilidad son modeladas por movimientos geométricos Brownianos. Una particularidad de la programación dinámica como herramienta de solución al modelo, es que este se resuelve analíticamente de manera determinista al resolver la EDP de HJB, pero con soluciones estocásticas. Por otra parte en el proceso de solución del modelo, a través de un sistema de ecuaciones diferenciales resultante se deduce una ecuación diferencial parcial de segundo orden equivalente a la de Black-Scholes-Merton con la que se valúa una opción asiática de venta de tipo europeo, con precio de ejercicio variable igual a la media aritmética mediate microfundamentos.

Se hace importante observar que la fórmula alternativa de valuación deducida en este trabajo no estaba disponible en la literatura, mediante fundamentos de racionalidad económica.

Para verificar las ventajas del modelo propuesto se desarrolla una metodología mediante simulación Montecarlo con los siguientes supuestos: la volatilidad del subyacente es estocástica, los parámetros se pueden obtener a partir de la volatilidad del subyacente, el precio del subyacente y la volatilidad σ t son conducidos por un movimiento geométrico browniano. El subyacente es el precio de la acción de WALMEXV y se comparan con los precios de opciones publicadas en el boletín del mercado mexicano de derivados (MEXDER) del día dieciséis de mayo de dos mil dieciocho. Se supone que los parámetros de volatilidad estocástica se calculan mediante una adaptación del modelo de CIR al método de volatilidad realizada. La tasa de interés se obtiene de la página de Banxico para los cuatro plazos de las opciones del boletín.

A partir de la volatilidad anualizada se determinaron los parámetros necesarios para la simulación del precio de la opción. Un resultado relevante es que, al comparar los precios de opciones europeas, tanto simulados como con los publicados por el Mercado Mexicano de Derivados con sus contrapartes asiáticas, los precios de opciones asiáticas son muy cercanos en el caso de opciones de compra y de venta dentro del dinero, en el dinero y fuera del dinero para el corto plazo, pero conforme el plazo al vencimiento de la opción se incrementa, la diferencia entre las europeas y las asiáticas aumenta ya que la prima es menor. Lo anterior se debe a que por el supuesto de volatilidad estocástica, la volatilidad del promedio aritmético tiene variaciones menos significativas, lo que genera diferencia en los precios de opciones europeas versus opciones asiáticas y esta diferencia se acentúa más en el largo plazo. Esto muestra que las opciones asiáticas presentan ventaja respecto a las opciones europeas. Para futuras líneas de investigación se plantea extender la valoración de opciones de tipo asiático con promedio geométrico, tasa de interés y/o volatilidad estocástica y de tipo americana, y en un contexto microeconómico y/o macroeconómico.