Estimación de la distribución multivariada de los rendimientos de los tipos de cambio contra el dólar de las criptomonedas Bitcoin, Ripple y Ether

Beatriz Mota Aragón1*; José Antonio Núñez2

1. Universidad Autónoma Metropolitana, México, Universidad Autónoma Metropolitana, Universidad Autónoma Metropolitana, Mexico , 2. Tecnológico de Monterrey, México , Tecnológico de Monterrey, Tecnológico de Monterrey, Mexico

Correspondence: *. E-mail:


Resumen

En este artículo se estima la distribución multivariada para analizar la dependencia del Bitcoin (BTC), Ripple (XRP) y Ether (ETH). Se utiliza la familia Hiperbólica Generalizada de distribuciones (GH) y en particular la distribución Varianza Gamma. El procedimiento para la estimación de los parámetros de la GH es a través del algoritmo EM (Expectation-Maximization). Los resultados muestran que existe una dependencia positiva entre los tres tipos de cambio respecto del dólar americano y seestima una distribución Varianza-Gamma de dimensión tres. Esta distribución es muy flexible para el ajuste de series de los rendimientos con leptocurtosis y sesgo. Esta información se considera importante para los inversionistas que conforman sus portafolios de una manera eficiente.

Received: 2018 November 28; Accepted: 2019 March 1

rmef. 2019 Jun 25; 14(3)
doi: 10.21919/remef.v14i3.409

Keywords: Palabras clave: leptokurtosis, asimetría, varianza gamma, multivariante.
Keywords: Clasificación JEL: C11, C46.
Keywords: Keywords: leptokurtosis, skewness, Variance-Gamma, Multivariate.
Keywords: JEL Classification: C11, C46.

1. Introducción

La tecnología blockchain y las criptomonedas han llamado la atención de una manera importante en la última década. De entre las más sobresalientes están el Bitcoin, que en la literatura es común encontrarla como la primera criptomoneda digital descentralizada, la cual desde su creación utiliza un conjunto de funciones criptográficas; sin embargo, en orden cronológico fue el eCash, el primer sistema criptográfico monetario electrónico (Chaum, 1983) y más tarde, el e-Gold el cual permitió la posibilidad de hacer transferencias instantáneas de valor (GSR, 1998), Henriques, et al. (2018). El Bitcoin creado por Satoshi Nakamoto en 2009, alcanzó la capitalización de mercado en abril del 2018 arriba de los 116 billones de dólares (coinmarketcap.com), y en general la capitalización del mercado de criptomonedas asciende los 295 billones de dólares con más de 1600 (Shuyue et. al. 2018) de ellas que han surgido con diferentes propiedades. Las características distintivas del Bitcoin son que ninguna autoridad central la respalde o controle y que su oferta está limitada por el diseño del protocolo, Bouri et al. (2017) y Dwyer (2015). Por su parte el Ether, la primera moneda basada en tecnología de una cadena de bloques o blockchain; por último, el Ripple distinta del bitcoin debido a que no utiliza el proceso de minería, sino que las transacciones se verifican entre los miembros de la red, son también criptodivisas destacables. Estos novedosos activos tienen nuevas características que los hacen únicos; tienen un control descentralizado en su operación y en su transacción. El uso del cryptocash es distinto al tradicional uso de efectivo, lo cual indica que las transacciones pueden ser realizadas directamente entre agentes económicos sinnecesidad de intermediación financiera y control de cambios, lo cual reduce los costos de transacción y los tiempos. En Sahoo (2017), se discute el futuro del Bitcoin como posible medio de cambio, y se concluye que esta sobrevivencia está sujeta principalmente a que tenga una volatilidad más moderada en el futuro, ya que esto originaría una mayor confianza en los usuarios de dicha criptomoneda. Aunque con datos diarios la volatilidad de los rendimientos del Bitcoin es mayor que los del oro y otras monedas en dólares americanos, es interesante el hallazgo en Dwyer (2015), encuentra que en ciertos periodos de tiempo, las volatilidades mensuales del Bitcoin son menores. De esta manera, no queda claro que de manera uniforme los rendimientos del Bitcoinsean mayores siempre en cualquier circunstancia.

Nuestro trabajo se interesa especialmente en estudiar la dependencia multivariada de los rendimientos de los tipos de cambio de las criptomonedas Bitcoin, Ripple y Ether contra el dólar, y contribuye al estudio y desarrollo de una explicación teórica de la mencionada dependencia.

La investigación se realiza a través de la familia hiperbólica generalizada (GHD) multivariada. En particular, se encuentra que con una distribución Varianza-Gamma de orden tres, es posible describir una estructura de dependencia positiva. El enfoque aplicado corresponde al de distribuciones paramétricas con una relación entre las características de los datos (hechos estilizados) y los momentos de la distribución, que están en función de los parámetros calculados.

Los resultados empíricos proveen la utilidad del trabajo. Por otro lado, planteamos como Corbet et al. (2018b), Glaser et al. (2014), Demir et al. (2018) la naturaleza de las criptomonedas como un medio de cambio e instrumento de inversión especulativo,en este sentido tiene implicaciones en la teoría de portafolios y la teoría de administración de riesgos y reconocemos la naturaleza dual de las criptomonedas, y la determinación de la tecnología y de los factores económicos como una línea de investigación futura, Gandal Halabrurda (2016), Xin, et al. (2016). De manera interesante existen investigaciones que no encuentran que la volatilidad del Bitcoin sea una consecuencia de la especulación, como puede verse en Blau (2017). Los resultados deesta investigación contribuyen a la comprensión del mercado altamente dinámico de las criptomonedas. El trabajo se divide en cinco secciones: En la sección dos se hace la revisión de literatura, en la tres se explica la metodología y los datos. En la sección cuatro exponemos los resultados del estudio, y la sección cinco se refiere a las conclusiones.

2. Revisión de literatura

Los factores que intervienen en la formación de precios del Bitcoin son radicalmente diferentes a los de activos convencionales. Estos factores incluyen noticias en redes, búsquedas en internet, Kristoufek (2013), tendencias en google (García et al., 2014), entre otras. La naturaleza del Bitcoin de acuerdo a Yemack (2013) es especulativa, y aunque en el mercado las principales criptomonedas son el Bitcoin y el Ripple, se sabe que la capitalización y la proporción del mercado de las criptomonedas ha variado fuertemente junto con su volatilidad (White, 2014). También Glaser et al. (2014) encuentra que la mayoría de los inversionistas en Bitcoin lo tratan con un fin especulativo, de lo cual deduce que es más un activo que una moneda. En Baur et al. (2017, 2018) consideran que el Bitcoin es un intermedio entre un metal precioso y una moneda fiduciaria. De manera interesante se establece que tiene una correlación débil o nula con activos tradicionales (acciones, bonos, oro etc.), y por tanto puede utilizarse como opción de diversificación en un portafolio de inversión.

En el sentido anterior, Dyhrberg (2015, 2016) estudia la capacidad de cobertura del Bitcoin contra los tipos de cambio del dólar americano con el Euro y la libra esterlina y concluye que es un valioso activo para períodos de crisis para un perfil de inversionista adverso al riesgo, teniendo estas características intermedias entre el oro y el dólar americano; es decir, entre la característica de una moneda fiduciaria como medio de cambio y la capacidad de guardar valor del oro. Es de notarse la falta de liquidez del Bitcoin con respecto a otros activos, y la presencia de correlación serial en la serie de sus rendimientos. El análisis econométrico se llevó a cabo con un GARCH asimétrico, y el caso del Bitcoin se halló que genera mayor volatilidadcuando hay períodos de crisis y disminuye cuando hay noticias buenas, lo cual es lo contrario de lo que sucede con acciones. Sin embargo, Baur et al. (2017) replica el estudio de Dyhrberg (2016) y encuentra que los rendimientos del Bitcoin siguen un proceso de volatilidad diferente a otros activos como el oro, y demuestra que se encuentra no correlacionado con otros activos como el dólar americano y el mismo oro. Respecto al estudio de las correlaciones, en Bouri et al. (2017) se analizan las propiedades de cobertura y de refugio (como el oro) del Bitcoin, usando un modelo de correlaciones condicionales dinámicas (DCC). Los autores concluyen que el Bitcoin es un activo útil para la diversificación, aunque esta propiedad varía con el tiempo, así como también las propiedades de ser un activo para cobertura y refugio seguro varían dependiendo del periodo de tiempo estudiado.

Como es conocido, la modelación de la volatilidad en finanzas es uno de los aspectos más importantes pues interviene en la administración del riesgo, en la formación de portafolios, en la estructura de formación de burbujas especulativas entre otros aspectos. La modelación de la volatilidad es uno de los aspectos más discutidos recientemente en la literatura. Katsiampa (2017) usó el modelo AR-CGARCH para estimar la volatilidad del Bitcoin y Urquhart (2017) concluyó que los modelos HAR son mejoresen estimar la volatilidad del Bitcoin que los modelos GARCH usuales.

Complementando lo anterior también se ha trabajado en relacionar variables macroeconómicas con la volatilidad del Bitcoin. Cemak (2017) usó un GARCH (1,1) para países con alto volumen de trading. Concluye que el Bitcoin se porta de manera similar a las monedas fiduciarias en China, Estados Unidos y Europa, pero no en Japón. En Chu et al. (2017) modelan la volatilidad de siete criptomonedas con doce modelos tipo GARCH, a decir, SGARCH, EGARCH, GJRGARCH, APARCH, IGARCH, CSGARCH, GARCH, TGARCH, AVGARCH, NGARCH, NAGARCH y ALLGARCH para el Bitcoin, Dash, Dogecoin, Litecoin, Maidsafecoin, Monero y Ripple. Los mejores modelos de acuerdo con estos autores fueron el IGARCH y el GJRGARCH, en donde es de destacar que el crecimiento del Bitcoin, Ether, Litecoiny Riple ha sido mayor que el resto de las criptomonedas. Los autores plantean en sus conclusiones que es necesario investigar y desarrollar el comportamiento multivariado de los tipos de cambio del Bitcoin, Dogecoin, Maidsafecoin, Monero y Ripple. En contraste, en Katsiampa (2017) se explica que la inclusión de las componentes de corto y largo plazo en la modelación de la volatilidad condicional, se encuentra que el modelo AR-CGARCH es el que mejor ajusta los datos del rendimiento del Bitcoin.

Además de las características antes mencionadas del Bitcoin y de las criptomonedas, los enfoques de econofísica han contribuido para ir más a fondo en el estudio de los orígenes de las burbujas y caídas desde la perspectiva monetarista, ver Corsi and Sornette (2014). Las burbujas y caídas para el caso de las criptomonedas tienen un origen diferente y representa un reto importante para la economía y los mercados financieros el descubrir su naturaleza. Las burbujas especulativas son estudiadas en Fry y Chea (2016) desde el enfoque de la econofísica, demostrando primeramente la utilidad de estas herramientas en el estudio del mercado de las criptomonedas y también sugieren su potencial uso en la cobertura y trading en el mercado financiero. Katsiampa (2017) menciona igualmente que el uso del Bitcoin como activo más que como moneda, lo hace especulativo y volátil, así como más susceptible a la generación de burbujas que las monedas. Sin embargo, en Blau (2017) se encuentra que la volatilidad delos rendimientos del Bitcoin no es consecuencia de la especulación con dicha criptomoneda, en el periodo de 2010 a 2013. El 2013 es un punto de referencia importante pues el valor del Bitcoin comenzó con un valor de 13 h a s t a l l e g a r a m s d e 1000. Elanterior estudio está basado en una medida de trading especulativo propuesta por Llorente et al. (2002).

Otra herramienta que es fundamental en diversas áreas de finanzas como la administración de riesgos financieros, la formación de portafolios de inversión, la simulación Monte Carlo, entre otras, es la estimación de la distribución de probabilidad. En Chu et al. (2015) se encuentra que la distribución que mejor se ajusta a los datos del rendimiento del Bitcoin vs. Dólar americano, es la hiperbólica generalizada, la cual tiene un excelente grado de flexibilidad debido a que usan cinco parámetros.

Adicional a lo anterior, Chan et al. (2017) estudian los rendimientos de tipo de cambio contra el dólar americano de las seis principales criptomonedas (de acuerdo a la capitalización del mercado son el Bitcoin, Dash, Dogecoin, Litecoin, MaidSafeCoin y Ripple), y se estiman las distribuciones paramétricas.

Los rendimientos de las criptomonedas exhiben muy alto valor de la curtosis y sesgo negativo a excepción del Bitcoin y el Litecoin. Las distribuciones encontradas son la hiperbólica generalizada, la normal inversa Gaussiana, t Generalizada y Laplace, todas ellas univariadas. Es decir, para este tipo de rendimientos existen colas pesadas que determinan de manera importante valores como la Valor en Riesgo y el VaR condicionado, así como la matriz de varianzas y covarianzas. En otro estudio, Zhang et al. (2018) revisa los hechos estilizados de los rendimientos del Bitcoin, Ether, Ripple, Litecoin, Stellar, Dash, Monero y NEM, donde se destaca el hecho de que todas tienen una curtosis mayor a 3 y la distribución no es normal al 1% de nivel de significancia. Las autocorrelaciones para los rendimientos absolutos decaen lentamente, y se prueban los fenómenos de apalancamiento con el modelo GJR(1,1) y cúmulos de volatilidad con un modelo GARCH(1,1). Es importante estudiar las distribuciones multivariadas de los rendimientos de diferentes criptomonedas, ya que se puede encontrar que existe una dinámica entre ellas, que es significativa, y que por tanto tiene un impacto importante en el diseño de portafolios y medición del riesgo entre otros aspectos.

3. Metodología y Datos

Datos

Los datos seleccionados son los rendimientos diarios de las criptomonedas Bitcoin, Ripple y Ether del 02 de febrero del 2016 al 15 de octubre de 2018, de investing: https://mx.investing.com/. Se calculará la distribución multivariada de dimensión tres.

Distribución de probabilidad y prueba de hipótesis

La distribución hiperbólica generalizada (GHD) tiene parámetros que describen los hechos estilizados de los rendimientos, en este caso lo aplicamos a los rendimientos de los tipos de cambio de criptomonedas contra el dólar (Protassov 2004 y Paolella 2007). La distribución se caracteriza por cinco parámetros en la función de densidad:


<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>;</mml:mo>
<mml:mi>λ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>α</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>β</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>δ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>λ</mml:mi>
<mml:mo>−</mml:mo>
<mml:mn>12</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>−</mml:mo>
<mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>δ</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>,</mml:mo>
<mml:msup>
<mml:mi>α</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msqrt>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>π</mml:mi>
</mml:mrow>
</mml:msqrt>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>λ</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msup>
<mml:mi>δ</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>,</mml:mo>
<mml:msup>
<mml:mi>α</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>−</mml:mo>
<mml:msup>
<mml:mi>β</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mtext>exp</mml:mtext>
<mml:mrow>
<mml:mo> [</mml:mo>
<mml:mrow>
<mml:mi>β</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>−</mml:mo>
<mml:mi>μ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>

con las restricciones α , δ 0 , λ , μ , β y β [ α , α ] y la función de Bessel de tercer orden,


<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>λ</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:munderover>
<mml:mstyle>
<mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mn>0</mml:mn>
<mml:mi>∞</mml:mi>
</mml:munderover>
<mml:msup>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>λ</mml:mi>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mi>x</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:msup>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>.</mml:mo>
</mml:mrow>

donde λ y x > 0 . La función de densidad de una variable aleatoria Gaussiana Inversa Generalizada (GIG), ver Barndorff-Nielsen (1977) es:


<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>;</mml:mo>
<mml:mi>λ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>χ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ψ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>λ</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>χ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ψ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mi>x</mml:mi>
<mml:mrow>
<mml:mi>λ</mml:mi>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mtext>exp</mml:mtext>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>χ</mml:mi>
<mml:msup>
<mml:mi>x</mml:mi>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>ψ</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>


<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>λ</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>χ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ψ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mn>2</mml:mn>
<mml:msup>
<mml:mi>η</mml:mi>
<mml:mrow>
<mml:mfrac>
<mml:mi>λ</mml:mi>
<mml:mn>2</mml:mn>
</mml:mfrac>
</mml:mrow>
</mml:msup>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>λ</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>ω</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>


<mml:mrow>
<mml:mi>η</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>χ</mml:mi>
<mml:mo>/</mml:mo>
<mml:mi>ψ</mml:mi>
</mml:mrow>


<mml:mrow>
<mml:mi>ω</mml:mi>
<mml:mo>=</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mi>χ</mml:mi>
<mml:mi>ψ</mml:mi>
</mml:mrow>
</mml:msqrt>
</mml:mrow>

Los parámetros λ , χ 0 y ψ 0 determinan la localización y la dispersión de x + ,

También con los siguientes cambios,


<mml:mrow>
<mml:mi>χ</mml:mi>
<mml:mo>:</mml:mo>
<mml:mo>=</mml:mo>
<mml:msup>
<mml:mi>δ</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>,</mml:mo>
</mml:mrow>


<mml:mrow>
<mml:mi>ψ</mml:mi>
<mml:mo>:</mml:mo>
<mml:mo>=</mml:mo>
<mml:msup>
<mml:mi>α</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>−</mml:mo>
<mml:msup>
<mml:mi>β</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>

Que el valor esperado


<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mi>X</mml:mi>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>μ</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>β</mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>λ</mml:mi>
<mml:mo>+</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>χ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ψ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>λ</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>χ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ψ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>

Y la varianza


<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mi>X</mml:mi>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>λ</mml:mi>
<mml:mo>+</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>χ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ψ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>λ</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>χ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ψ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>β</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>λ</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>χ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ψ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>λ</mml:mi>
<mml:mo>+</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>χ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ψ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>−</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>λ</mml:mi>
<mml:mo>+</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>χ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ψ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>λ</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>χ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ψ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>

Para el caso de una distribución hiperbólica generalizada multivariada, tomemos un vector aleatorio X de tamaño n × 1 y diremos que sigue una función de densidad hiperbólica generalizada multivariada (GH), si esta función se expresa como


<mml:mi>f</mml:mi>
<mml:mfenced>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>;</mml:mo>
<mml:mi>λ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>γ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>β</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>μ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>δ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>Δ</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msqrt>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi>δ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mfenced>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>μ</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:msup>
<mml:mrow>
<mml:mi>Δ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mfenced>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>μ</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mi>'</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi>γ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>β</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>'</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mi>Δ</mml:mi>
<mml:mi>β</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:msqrt>
<mml:msup>
<mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mi>λ</mml:mi>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mfenced>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>π</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:msup>
<mml:msup>
<mml:mrow>
<mml:mfenced>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>γ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>δ</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mi>λ</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>K</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>λ</mml:mi>
<mml:mo>-</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:msub>
<mml:msqrt>
<mml:msup>
<mml:mrow>
<mml:mi>δ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mfenced>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>μ</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:msup>
<mml:mrow>
<mml:mi>Δ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mfenced>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>μ</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mi>'</mml:mi>
<mml:mo>(</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>γ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>β</mml:mi>
<mml:mi>'</mml:mi>
<mml:mi>Δ</mml:mi>
<mml:mi>β</mml:mi>
<mml:mo>)</mml:mo>
</mml:msqrt>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>K</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>λ</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>(</mml:mo>
<mml:mi>δ</mml:mi>
<mml:mi>γ</mml:mi>
<mml:mo>)</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mi>exp</mml:mi>
</mml:mrow>
<mml:mo>⁡</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>-</mml:mo>
<mml:mo>(</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>μ</mml:mi>
<mml:mo>)</mml:mo>
<mml:mi>β</mml:mi>
<mml:mi>'</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mfrac>

con χ = δ 2 y ψ = γ 2 β ' Δ β , los parámetros de una variable GIG. También


<mml:mrow>
<mml:mi>X</mml:mi>
<mml:mo>|</mml:mo>
<mml:mi>W</mml:mi>
<mml:mo>∼</mml:mo>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>μ</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>w</mml:mi>
<mml:mi>β</mml:mi>
<mml:mtext>Δ</mml:mtext>
<mml:mo>,</mml:mo>
<mml:mi>w</mml:mi>
<mml:mtext>Δ</mml:mtext>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>


<mml:mrow>
<mml:mi>W</mml:mi>
<mml:mo>∼</mml:mo>
<mml:mi>G</mml:mi>
<mml:mi>I</mml:mi>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>λ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>χ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ψ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>

La matriz de covarianzas para el vector aleatorio X es V [ X ] = E [ W ] Δ + V [ W ] β β y el valor esperado es E [ X ] = μ + E [ W ] β (Hu, 2005).

Respecto a la interpretación de los parámetros tenemos que λ está relacionado con qué tan planas son las colas de la distribución, χ y ψ de manera conjunta determinan la forma de la distribución (por ejemplo cuando estos parámetros toman valores grandes en su dominio, la distribución se aproxima a la normal), μ nos dice dónde está centrada la distribución, y finalmente γ es el parámetro que indica el sesgo, Konlack y Wilcox (2014)

Estimación

El procedimiento para la estimación de los parámetros de la GH es a través del algoritmo EM (Expectation-Maximization), ver Protassov (2004) y Breymann (2013). En particular, la distribución Varianza-Gamma también puede verse de la siguiente manera:


<mml:mi>f</mml:mi>
<mml:mfenced>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mo>=</mml:mo>
<mml:mi>c</mml:mi>
<mml:mfenced>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>σ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>θ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ϑ</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:msup>
<mml:mrow>
<mml:mfenced>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:mi>ϑ</mml:mi>
<mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mfenced>
<mml:mrow>
<mml:mi>θ</mml:mi>
<mml:mfenced>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mo>/</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>σ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msup>
<mml:msub>
<mml:mrow>
<mml:mi>K</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:mi>ϑ</mml:mi>
<mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mfenced>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mfenced>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:msqrt>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msup>
<mml:mrow>
<mml:mi>σ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mi>ϑ</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>θ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:msqrt>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi>σ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>

donde K ϑ es la función de Bessel de tercer orden,


<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>σ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>θ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>ϑ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>2</mml:mn>
<mml:mrow>
<mml:mi>σ</mml:mi>
<mml:msqrt>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>π</mml:mi>
</mml:mrow>
</mml:msqrt>
<mml:msup>
<mml:mi>ϑ</mml:mi>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:mi>ϑ</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mi>Γ</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:mi>ϑ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msup>
<mml:mi>σ</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mi>ϑ</mml:mi>
</mml:mfrac>
<mml:mo>+</mml:mo>
<mml:msup>
<mml:mi>θ</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>ϑ</mml:mi>
</mml:mfrac>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>

Estimación de la función de densidad

Si el vector de parámetros está definido como Θ = ( λ , γ , δ , β , μ , Δ ) , el logaritmo de la función de verosimilitud es


<mml:mrow>
<mml:mtext>ln</mml:mtext>
<mml:mo>(</mml:mo>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtext>Θ</mml:mtext>
<mml:mo>;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>…</mml:mo>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:munderover>
<mml:mstyle>
<mml:mo>∑</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>m</mml:mi>
</mml:munderover>
<mml:mtext>ln</mml:mtext>
<mml:mo>(</mml:mo>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>;</mml:mo>
<mml:mtext>Θ</mml:mtext>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>

donde x 1 , , x m es una muestra de vectores de números aleatorios. El algoritmo EM considera variables latentes w 1 , w 2 , , w m con una distribución GIG. La función de verosimilitud expandida es


<mml:mrow>
<mml:mtext>ln</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msup>
<mml:mi>L</mml:mi>
<mml:mo>′</mml:mo>
</mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtext>Θ</mml:mtext>
<mml:mo>;</mml:mo>
<mml:mo> </mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>…</mml:mo>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>w</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>w</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>…</mml:mo>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>w</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:munderover>
<mml:mstyle>
<mml:mo>∑</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>m</mml:mi>
</mml:munderover>
<mml:mtext>ln</mml:mtext>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>X</mml:mi>
<mml:mo>|</mml:mo>
<mml:mi>W</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>|</mml:mo>
<mml:msub>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>;</mml:mo>
<mml:mi>μ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mtext>Δ</mml:mtext>
<mml:mo>,</mml:mo>
<mml:mi>β</mml:mi>
<mml:mo>)</mml:mo>
<mml:mo>)</mml:mo>
<mml:mo>+</mml:mo>
<mml:munderover>
<mml:mstyle>
<mml:mo>∑</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>m</mml:mi>
</mml:munderover>
<mml:mtext>ln</mml:mtext>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>W</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>;</mml:mo>
<mml:mi>λ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>γ</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>δ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>

Después de la selección de un punto inicial Θ 0 , se desarrollan los siguientes pasos:

  1. a) Estimación del valor esperado de la función del logaritmo de la función de máxima verosimilitud aumentada,
  2. h ( Θ ; Θ k ) = E [ ln ( L ( Θ ;   x 1 , , x m , w 1 , w 2 , , w m ) ) | x 1 , , x m ; Θ k ]
  3. b) Maximizar la función objetivo h ( Θ ; Θ k ) para obtener la aproximación
  4. Θ k + 1 hasta que | Θ k + 1 Θ k | < ε .

El anterior algoritmo se programó en el paquete R. Se probaron diversos puntos iniciales del método numérico, entre ellos la media y mediana de cada serie de datos. El mismo programa ghyp en R realiza diversas pruebas con valores iniciales, Breymann Lüthi (2013).

Se puede demostrar que el estimador Θ ^ tiene distribución normal multivariada asintótica, es decir,


<mml:mrow>
<mml:msqrt>
<mml:mi>m</mml:mi>
</mml:msqrt>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover>
<mml:mi>Θ</mml:mi>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>−</mml:mo>
<mml:msub>
<mml:mi>Θ</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>→</mml:mo>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>0,</mml:mn>
<mml:msup>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>

donde


<mml:mrow>
<mml:mi>I</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>C</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>


<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>∂</mml:mo>
<mml:mi>ln</mml:mi>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mover>
<mml:mi>Θ</mml:mi>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mo>∂</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mover>
<mml:mi>Θ</mml:mi>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>

El intervalo para el estimador (Barndorff-Nielsen, 2012), Θ ^ i es


<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover>
<mml:mi>Θ</mml:mi>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>±</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>α</mml:mi>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>m</mml:mi>
</mml:mfrac>
<mml:msub>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msup>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover>
<mml:mi>Θ</mml:mi>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mrow>
<mml:mover>
<mml:mi>Θ</mml:mi>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msqrt>
</mml:mrow>


<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover>
<mml:mi>Θ</mml:mi>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>−</mml:mo>
<mml:mi>Θ</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:msqrt>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover>
<mml:mi>Θ</mml:mi>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mfrac>
</mml:mrow>

Para un nivel de significancia α . La bondad de ajuste entre una muestra de datos y la GH multivariada se desarrolla mediante la prueba de McAssey (2013).

La prueba de McAssey consiste en dos pasos:

  1. I) Simulación de una muestra aleatoria GH, { u ^ 1 , u ^ 2 , , u ^ T } , con los parámetros estimados a partir de algoritmo EM (Ross, 2012).
  2. II) Cálculo de la distancia de Mahalanobis entre la muestra de datos y a i) con la media muestral μ ^ y la matriz de covarianzas Σ ^ .


<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover>
<mml:mi>d</mml:mi>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover>
<mml:mi>u</mml:mi>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>−</mml:mo>
<mml:mover>
<mml:mi>μ</mml:mi>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>'</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mover>
<mml:mi>Σ</mml:mi>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mo>−</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover>
<mml:mi>u</mml:mi>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>−</mml:mo>
<mml:mover>
<mml:mi>μ</mml:mi>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:msqrt>
</mml:mrow>

Intuitivamente, la distancia de Mahalanobis coincide con la distancia Euclidiana si los vectores que estamos utilizando no están correlacionados. Por otro lado, si existe correlación entre los vectores, entonces se pondera por medio de la matriz de covarianzas Σ ^ , para obtener la distancia entre dichos vectores.

Por tanto, se obtiene una muestra aleatoria { A 1 , A 2 , , A p } con


<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:munderover>
<mml:mstyle>
<mml:mo>∑</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>T</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mo>|</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>−</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>O</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>|</mml:mo>
</mml:mrow>
</mml:mrow>

donde E j es la frecuencia esperada de las observaciones en un intervalo fijo, O j es la frecuencia observada en dicho intervalo y el valor p de la prueba es calculado.

4. Resultados

Para las series de tiempo de las criptomonedas Bitcoin, Ether y Ripple (rendimientos), se obtuvieron las estadísticas de la tabla 1. Es de notarse la alta curtosis de los rendimientos del Bitcoin y del Ripple. Existe un sesgo negativo en los tres casos.

Tabla 1.

Estadística Descriptiva


Variable Promedio Desviación Estándar Sesgo Curtosis Mínimo Máximo
BITCOIN 0.0029 0.0460 -0.5880 10.2213 -0.3716 0.3083
ETHER 0.0031 0.0631 -0.0227 2.9306 -0.3134 0.2586
RIPPLE 0.0026 0.1743 -1.3621 85.3564 -2.4079 2.3026

Aplicando las pruebas de raíces unitarias obtenemos que ninguna de las series de rendimientos presenta raíces unitarias (Tabla 2). La gráfica de los rendimientos se muestra también (gráfica 1). Esta gráfica ciertamente confirma lo obtenido las pruebas de Dickey-Fuller, Phillips-Perron y KPSS del paquete EVIEWS. Por tanto, podemos continuar con el cálculo de los parámetros de la distribución de dimensión tres.

Tabla 2.

Pruebas de Raíces Unitarias


Variable

  1. Dickey Fuller
  2. Estadístico t

  1. Phillips-Perron
  2. Estadístico t

  1. KPSS
  2. Estadístico t

BITCOIN -49.33*** -49.34*** 0.21
ETHER -30.69*** -30.75*** 0.35
RIPPLE -17.30*** -17.65*** 0.15

TFN1***: Nivel de significancia de 1%,

TFN2**: Nivel de significancia de 5%,

TFN3*: Nivel de significancia de 10%



[Figure ID: f1] Figura 1.

Para la estimación de los parámetros se usó el algoritmo EM, de tal manera que una vez que se tienen los parámetros, probamos si la distribución es normal multivariada o Varianza-Gamma. De acuerdo con la prueba de McAssey (2013), se rechaza la hipótesis de normalidad de los datos y no se rechaza la hipótesis de una distribución Varianza-Gamma, Tabla 3. Algunas de las ventajas que tiene la distribución Varianza-Gamma son que incluye colas más pesadas que las de la distribución normal y todos losmomentos de la distribución son finitos entre otras propiedades, Madan y Seneta (1990). Las distribuciones varianza-gamma forman una subclase de la distribución hiperbólica propuesto por O.E. Barndorff-Nielsen (1977).

Podemos notar de la matriz de covarianzas de la distribución de orden 3, que existe una dependencia positiva lo cual nos proporciona información importante; por ejemplo, en la formación de portafolios, pues no podemos diversificar en gran medida con las tres criptomonedas, pues existe una relación en el mismo sentido. La diversificación de un posible portafolio se tiene que buscar a través de otros activos. También hay que decir que en este punto no podemos determinar que estas tres criptomonedas sean usadas para especulación solamente o al menos de manera importante.

Tabla 3.

Estimación de los parámetros de la distribución de dimensión tres


Log-verosimilitud 4798.77 λ 0.7510 (1.68)*

<mml:mi>μ</mml:mi>

  1. 0.0026
  2. (-1.82)*

  1. -0.0023
  2. (-1.73)*

  1. -0.0026
  2. (1.93)*


<mml:mtext>Δ</mml:mtext>
BITCOIN ETHER RIPPLE
BITCOIN

  1. 0.0019
  2. (1.96)**

  1. 0.0014
  2. (1.83)*

  1. 0.0012
  2. (2.03)**

ETHER

  1. 0.0014
  2. (1.83)**

  1. 0.0043
  2. (2.05)*

  1. 0.0020
  2. (1.82)*

RIPPLE

  1. 0.0012
  2. (2.03)*

  1. 0.0020
  2. (1.82)*

  1. 0.0051
  2. (1.95)*


<mml:mi>γ</mml:mi>

  1. 0.0004
  2. (1.96)**

  1. 0.0054
  2. (1.88)*

  1. 0.0069
  2. (1.47)

Prueba de McAssey (Normal) 1368.35 Valor-p 0.0000
Prueba de McAssey (VG) 1133.36 Valor-p 0.1021

TFN4***: Nivel de significancia de 1%; ** : Nivel de significancia de 5%; * : Nivel de significancia de 10%.


Como se observa de la matriz de varianzas y covarianzas para la distribución de dimensión tres, resaltan las covarianzas positivas y pequeñas, lo cual tiene consecuencias importantes en términos de inversión. Por un lado, nos indica que habría limitaciones para invertir en las tres criptomonedas en un mismo portafolio, pues hay limitaciones para generar un portafolio diversificado. Por otro lado, las características de volatilidad de las tres criptomonedas puede deberse a una fuerte especulación en este mercado, aunque falta desarrollar estudios más profundos para investigar este último punto, ya que se ha visto para ciertos periodos de tiempo que la alta volatilidad característica no es necesariamente consecuencia de la especulación. La prueba de bondad de ajuste de MacAssey nos dice que los datos se ajustan a la distribución varianza gamma de una manera razonable. La normalidad de los datos es rechazada como era de esperarse para datos con colas pesadas, debido a la volatilidad existente.

5. Conclusiones

El surgimiento de las criptomonedas ha significado la aparición de lo que llamaríamos un activo y que tiene características de moneda. Sin embargo, no está clara su naturaleza económica y no se ha definido lo que es. A pesar de lo anterior, existe unmercado importante de las criptomonedas,- sobre todo el Bitcoin-, y es un gran reto analizar y modelar estadísticamente la dependencia de este con el Ripple y el Ether. Hemos encontrado una dependencia positiva y pequeña, a través de una distribuciónde colas pesadas. La matriz de correlaciones es un ingrediente determinante de la formación de portafolios a través de la diversificación. Las distribuciones Varianza-Gamma forman una subclase de la familia hiperbólica generalizada. En la literatura,puede ser que éstos sean utilizados con fines de diversificación, cobertura o como refugio seguro (mantiene valor). De acuerdo con Baur y Lucey (2010), y con Ratner y Chiu (2013) un activo se puede usar para diversificación si mantiene una correlación positiva débil con otro activo en promedio. Recordemos que la volatilidad de estos activos financieros es muy alta y esa es una oportunidad para los especuladores o tomadores de riesgo. Sin embargo, dependiendo del periodo de tiempo, no necesariamente la alta volatilidad es debida a la especulación con éstos activos. Lo anterior está publicado en algunos artículos de la literatura. Falta esperar la evolución de este interesante mercado que combina tecnología, finanzas y características que no han permitido definir qué tipo de objeto es. Una extensión del trabajo es investigar si la volatilidad alta que caracteriza a este tipo de activos es debido realmente a especulación, o bajo qué condiciones es consecuencia de la especulación. Análogamente, si también la existencia de burbujas en este mercado nos indicaría que la naturaleza de estos activos es especulativa, y si eso es lo que mantiene el mercado activo.

Otros estudios han publicado que existe correlación cero de los rendimientos del Bitcoin con otros activos comunes como las acciones, bonos, commodities etc. Esta es otra ventana de oportunidad, puesto que se puede estudiar la correlación de otras criptomonedas con éstos activo más comunes, a partir de una distribución multivariada.


Referencias
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
18.
19.
20.
21.
22.
23.
24.
25.
26.
27.
28.
29.
30.
31.
32.
33.
34.
35.
36.
37.

Enlaces refback

  • No hay ningún enlace refback.


Licencia de Creative Commons
Este obra está bajo una licencia de Creative Commons Reconocimiento-NoComercial 4.0 Internacional.
Métricas de artículo
Cargando métricas ...

Metrics powered by PLOS ALM