VWO ARIMA
Sandra Gonzalez Arzate
07 Abril de 2020
# VWO Vanguard FTSE Emerging Markets ETF
VWO
Vanguard FTSE Emerging Markets ETF Es una de las compañías gestoras más grandes del mundo. Fundada el 1 de mayo de 1975. El primer fondo que estableció fue el Wellington Fund, cuya fecha de incepción fue el primero de julio de 1929. Actualmente administra más de 2 trillones de dólares en fondos internacionales y estadounidenses, según cifras al 30 de septiembre de 2013, con más de 170 fondos establecidos. Cuenta con oficinas en Estados Unidos, Holanda, Australia, Hong Kong, Francia, Inglaterra, Singapur, Japón, Canadá y Suiza.
VWO se refiere a Vanguard FTSE Emerging Markets ETF El fondo busca replicar el rendimiento de un índice dereferencia que mide el rendimiento de la inversión de acciones emitidas por empresas ubicadas en países de mercados emergentes..
ETF de Vanguard® Las acciones no son canjeables con el fondo emisor que no sea en agregaciones muy grandes que valen millones de dólares. En su lugar, los inversores deben comprar y vender Vanguard ETF Shares en el mercado secundario y mantener esas acciones en una cuenta de corretaje. Al hacerlo, el inversionista puede incurrir en comisiones de corretaje y puede pagar más que el valor neto del activo al comprar y recibir menos del valor neto del activo al vender.
Composición de VWO
Comportamiento del precio de VWO
Apartir del 02 de enero de 2015 al 30 de agosto de 2019. La tendencia que presenta la emisora de enero de 2015 a enero de 2018 es a la alza llegando a registrar un máximo de 51.25 USD por acción.
A continuación, se presentan el comportamiento del precio de cierre de VWO a partir del 1 de enero de 2015 al 30 de marzo de 2020. EDZ presentó su mayor crecimiento a finales de enero de 2018, esto está relacionado con la incertidumbre que se generó en torno a las elecciones presidenciales de Estados Unidos llegando a los $51usd. Posterior a las elecciones, los índices de las economías emergentes comenzaron a recuperar en incluso llegaron a registrar máximos históricos, como fue el caso de Índice de Precios y Cotizaciones en México (IPC) en agosto de 2017 cuando superó los 50 mil puntos base. VWO, al ser un instrumento apalancado que paga el inverso de los índices, presentó un fuerte declive, tal es el caso de VWO que presentó una caída hasta los 30usd por acción incluso libernado un gran volumen de acciónes a lo largo de su ultimo trimestre. La Figura 1 muestra el comportamiento del precio de cierre de EDZ.
Figura 1. Precio de Cierre de VWO y S&P500: enero 2015 - marzo 2020
Fuente: elaboración propia con datos de Yahoo Finance.
Si se compara este comportamiento con la tendencia registrada por el índice S&P500, que representa las 500 emisoras más representativas que cotizan en Estados Unidos, se puede ver un comportamiento inverso a VWO, esto tiene que ver con la fuerza de arrastre y contagio que tiene el S&P500 sobre los índices de economías emergentes sobre todo en VWO ya que su valor se encuentra únicamente en Emerging Markets, lo cual se puede corroborar en la Figura 2.
En lo que refiere a los rendimientos, el S&P500 en promedio osciló entre 2% en tanto que EDZ estuvo fluctuando entre 17% debido a que es un instrumento que presenta mayor volatilidad, es decir, presenta más dispersión en sus precios.
A partir de la contingencia sanitaria COVID 19 que se originó en la ciudad de Wuhan, China, que se dispersó en diferentes países ocasionando contagios y situaciones de emergencia internacional cerrando fronteras.En nuestro caso a partir del inicios de 2020, no se vió beneficiado el tipo de ETF porque los mercados emergentes fueron los primeros en resentir la caida del mercado bursátil.
Figura 2. Rendimientos de S&P500 y VWO: enero de 2015 a marzo 2020 
Fuente: elaboración propia con salida de R.
Correlación entre VWO y S&P500
La correlación que hay entre los precios de cierre a niveles y entre los rendimientos de VWO y el S&P500, muestra una alta relación inversa del 65% respectivamente.
Tabla 2. Correlación de VWO y S&P500 a niveles
SP500 VWO
SP500 1.0000000 0.6519263
VWO 0.6519263 1.0000000Fuente: elaboración propia con salida de R.
Tabla 3. Correlación de WVO y S&P500 en rendimientos
SP500 VWO
SP500 1.0000000 0.8105438
VWO 0.8105438 1.0000000Fuente: elaboración propia con salida de R.
Histogramas y gráficos Q-Q
Los histogramas son gráficos que representan frecuencia de un fenómeno o de una variable mediante una distribución de los datos. En el caso de VWO y del S&P500, a partir de los intervalos o marcas de clase que se hacen sobre ellos, se puede identificar el número de veces (frecuencia) que los precios caen en dicho intervalo.
En la figura 3 se presentan los histogramas a niveles del S&P500 y del VWO; el eje vertical representan las frecuencias y en el eje horizontal los valores de las variables (puntos base y precios respectivamente).
El histograma del S&P500 a niveles indica que, en el periodo de muestra, el índice tuvo mayor número de repeticiones en los 2100 puntos (48 veces). Sin embargo, la mayor parte de la distribución se centra entre los 2000 y los 3000 puntos. Para el caso de VWO, su mayor número de repeticiones es 400 para el precio de 42usd. Los valores más extremos y con pocas repeticiones se sitúan en más de los 300usd.
Figura 3. Histogramas a niveles S&P500 y VWO: enero de 2015 a marzo 2020
Fuente: elaboración propia con salida de R.En lo que refiere a los rendimientos, en promedio, los rendimientos presentan un proceso de reversión a la media (0), sin embargo, la distribución de los rendimientos del S&P500 oscila entre 2% y los de VWO en 10%, tal cual como se observó en la figura 2. La figura 4 presenta el histograma de los dos instrumentos que se han estado analizando en rendimientos.
Figura 4. Histogramas en rendimientos S&P500 y VWO: enero de 2015 a marzo 2020
Fuente: elaboración propia con salida de R.Los gráficos Cuantil-Cuantil (también referidos como q-q plots) es la representación gráfica de la distribución de un conjunto de datos a alguna distribución ideal o a priori que se asume como dada. La finalidad de estos gráficos es comparar la distribución teórica (la que suponemos o queremos como se comporte) contra la distribución empírica (la que realmente presenta la variable).
El siguiente gráfico muestra Q-Q del S&P y de VWO; los cuantiles teóricos o la distribución contra la que se están comparando los precios es contra una distribución normal; si la distribución empírica fuera así, entonces los puntos de dispersión deberían de distribuirse en torno a la recta.
Lo que se observa es que sí hay una parte de la distribución que se asocia a la línea recta, sin embargo, son más los datos, sobre todo en los extremos o en las colas, donde la distribución se “despega” de la normalidad.
Figura 5. Q-Q plot a niveles S&P500 y VWO: enero de 2015 a marzo 2020
Fuente: elaboración propia con salida de R.Lo mismo se observa en el caso del gráfico Q-Q de los rendimientos, sin embargo, en este ejemplo, nótese que los datos, al menos en la parte central de la distribución, están más pegados a la recta, esto tiene que ver con la propiedad que cumplen los rendimientos (media cero o constante que es uno de los supuestos que se debe de cumplir para la estacionariedad de las series), sin embargo, ambos instrumentos tuvieron días que presentaron rendimientos que rebasaron su media, provocando mayor dispersión en sus datos.
Con esta representación, no se puede garantizar la normalidad en los datos, y en lo que respecta a los instrumentos financieros, lo más normal es que no sean normales.
Figura 5. Q-Q plot en rendimientos S&P500 y VWO: enero de 2015 a marzo 2020
Fuente: elaboración propia con salida de R.Estacionariedad y pruebas de raices unitarias
El concepto de estacionariedad es importante para la estimación y para la elaboración de pronósticos, el no garantizar esta condición implicaría que las series, no serían independientes e idénticamente distribuidas, ocasionado problemas de sesgo en las estimaciones, regresiones espurias o el mal cálculo de las bandas de confianza a partir de datos que se encuentran correlacionados.
Las pruebas de raíces unitarias permiten identificar si la serie es estacionaria o no, verificando si la serie tiene alguna estructura de dependencia con los datos anteriores. Al pronosticar series de tiempo, se asumen que estas son aleatorias, por lo tanto:
Ecuación 1
\[E\left ( Y_t \right| \phi_t )=0\]
Donde \(Y_t\) es el valor esperado de la variable condicionado a \(\phi_t\), que refiere a la información pasada o registrada de la misma variable. Si esta variable es aleatoria, entonces su valor esperado es 0. La ecuación 1 también se le conoce como un proceso estocástico y en este caso, los precios se comportan de manera aleatoria, es decir:
Ecuación 2
\[f\left ( Y_t \right| Y_{t-1} )=f ( Y_{t} )\]
Cuando llega nueva información, los precios de las acciones fluctuarán aleatoriamente, al menos así lo dice la teoría.
Adicional al supuesto de la ecuación 1, las condiciones de estacionariedad también implican que las series sean homocedásticas, es decir, que su varianza sea constante. Este supuesto es difícil de cumplir para las series financieras debido a la dispersión o volatilidad que presentan los datos, sin embargo, de este supuesto nos encargaremos después.
Lo primero que se requiere garantizar es que la serie no tenga problemas de raíces unitarias, para que al menos se pueda garantizar el primer supuesto (valor esperado = 0).
Pruebas de raíces unitarias
Las pruebas que se utilzian para detectar raíces unitarias en este análisis son: Dickey Fuller Aumentada (DFA), Phillips Perron y la prueba Kwiatkowski - Phillips - Schmidt - Shin (KPSS). La tabla 4 muestra los resultados de S&P500 y de EDZ a niveles y rendimientos.
Tabla 4. Pruebas de raíces unitarias
| Variable | \(DFA^{a/}\)(Valor p) | \(Phillips-Perron^{b/}\)(Valor p) | \(KPSS^{c/}\)(Valor p) |
|---|---|---|---|
| S&P500 (a niveles) | 0.4464 | 0.5515 | 0.01 |
| S&P500 (rendimientos) | 0.01 | 0.01 | 0.1 |
| VWO (a niveles) | 0.6177 | 0.488 | 0.01 |
| VWO (rendimientos) | 0.01 | 0.01 | 0.1 |
\(^{a/}H0\): La serie tiene raíz unitaria
\(^{b/}H0\): La serie tiene raíz unitaria
\(^{c/}H0\): La serie es estacionaria
Fuente. Elaboración propia con salida de R.
NUNCA OLVIDAR:
Si valor p mayor a 0.05 No rechazo (acepto) H0.
Si valor p menor a 0.05 Rechazo H0.
¿Por qué la serie en rendimientos no tiene raíz unitaria?
Se debe a lo siguiente:
\[Y_t = \alpha + \beta Y_{t-1} + e_t\]
Suponga \(\beta=1\).
\[Y_t = \alpha + Y_{t-1} + e_t\]
Donde \(Y_t\) depende del valor pasado \(Y_{t-1}\), si esto es cierto, entonces la serie no es aleatoria, hay dependencia con el dato anterior y no podemos cumplir con el primer supuesto (ecuación 1).
A este proceso se le conoce también como: “caminata aleatoria”.
Se aplican primeras diferencias en ambas partes de la ecuación.
\[Y_t - Y_{t-1} = \alpha + \beta Y_{t-1} - Y_{t-1} + e_t\]
\[\Delta Y_t= \alpha + Y_{t-1}(\beta -1) + e_t\]
Recordemos que \(\beta=1\).
\[\Delta Y_t= \alpha + Y_{t-1}(1 -1) + e_t\]
\[\Delta Y_t= \alpha + Y_{t-1}(0) + e_t\] \[\therefore \]
\[ \Delta Y_t= \alpha + e_t\]
La serie, en primeras diferencias, no tiene raíz unitaria, solo depende del error y del intercepto, pero no de los valores pasados o registrados del precio, por lo tanto, es estacionaria.
A este proceso también se le conoce como “ruido blanco”.
Modelos ARIMA (VWO)
Ahora, se va a calcular el primer modelo ARIMA para hacer los pronósticos, utilizando la metodología de Box & Jenkins.
Se obtiene la Función de Autocorrelación (MA) y Función de Autocorrelación parcial (AR). Ambas series requieren ser integrada de orden I, es decir, se les tiene que aplicar una primera diferencia para que al menos puedan ser estacionarias en media. La aplicación de la primera diferencia es congruente con los resultados de las pruebas unitarias, en donde es necesario que las series se transformen en rendimientos.
Figura 6. Componentes de autocorrelación ACF y PACF 
Fuente: elaboración propia con salida de R.
Al revisar el correlograma (a pesar de diferenciar una vez la serie), se identifican componentes de autocorrelación tanto en el procero Autorregresivo (PACF) y en el proceso de media móvil (ACF).
El primer ajuste que se hace para el pronóstico de VWO es utilizando la función auto.arima de R, que propone una combinación de ARIMA(4,1,0) para corregir los problemas de autocorrelación.
Tabla 5. Resultados del ARIMA(4,1,0) para VWO
Series: VWO
ARIMA(5,1,2)
Coefficients:
ar1 ar2 ar3 ar4 ar5 ma1 ma2
-1.6203 -0.8461 0.0403 0.0688 0.0466 1.5403 0.7698
s.e. 0.1482 0.1623 0.0580 0.0559 0.0343 0.1458 0.1433
sigma^2 estimated as 0.2487: log likelihood=-944.58
AIC=1905.16 AICc=1905.27 BIC=1946.59Fuente: elaboración propia con salida de R.
Figura 7. Resultados del ARIMA(4,1,0) para VWO 
Ljung-Box test
data: Residuals from ARIMA(5,1,2)
Q* = 9.7193, df = 23, p-value = 0.9929
Model df: 7. Total lags used: 30Fuente: elaboración propia con salida de R.
El resultado muestra que no se han terminado de corregir los problemas de autocorrelación. Aplicando la prueba de Ljung-Box, donde la H0 es: los datos se distribuyen de forma independiente o dicho de otra forma, los residuales del ARIMA no están correlacionados. Para el ARIMA(4,1,0) la H0 se rechaza. Si bien se puede realizar un pronóstico con estos resultados, se cae el riesgo de obtener resultados sesgados (debido a los problemas de autocorrelación).
A continuación, se muestra la estabilidad del modelo a partir del gráfico de raíces uniarias, tanto en el proceso AR como en el de MA.
Figura 8. Prueba de racíces unitarias ARIMA(4,1,0) - círculo unitario 
Fuente: elaboración propia con salida de R.
Si bien el modelo es estable, los problemas de autocorrelación no han sido solventados en su totalidad, por lo que se propone el siguiente modelo.
Propuesta de modelo ARIMA(5,1,21) para VWO
Este modelo mejora significativamente los resultados propuestos por el ARIMA, se corrigen los problemas (en su mayoría) de autocorrelación en los residuales de acuerdo a los resultados de la prueba de Ljung-Box.
Tabla 6. Resultados del ARIMA(5,1,21) para VWO
Call:
arima(x = VWO, order = c(5, 1, 21))
Coefficients:
ar1 ar2 ar3 ar4 ar5 ma1 ma2 ma3
-0.4127 -0.5439 -0.3526 -0.6373 -0.8027 0.3287 0.5711 0.3519
s.e. NaN NaN NaN NaN NaN NaN NaN NaN
ma4 ma5 ma6 ma7 ma8 ma9 ma10 ma11 ma12
0.6432 0.8511 -0.0539 0.1260 -0.0392 0.0632 0.0114 -0.0385 0.0409
s.e. NaN NaN 0.0476 0.0467 0.0415 0.0422 0.0430 0.0467 0.0517
ma13 ma14 ma15 ma16 ma17 ma18 ma19 ma20 ma21
-0.0516 0.0649 -0.0286 0.0244 0.0634 0.0067 0.0431 0.0124 0.0521
s.e. 0.0516 0.0512 0.0507 0.0505 0.0447 0.0397 0.0345 0.0323 0.0280
sigma^2 estimated as 0.2429: log likelihood = -934.99, aic = 1923.99Fuente: elaboración propia con salida de R.
Figura 9. Resultados del ARIMA(5,1,21) para VWO 
Ljung-Box test
data: Residuals from ARIMA(5,1,21)
Q* = 4.7069, df = 4, p-value = 0.3187
Model df: 26. Total lags used: 30Fuente: elaboración propia con salida de R.
El modelo, a pesar de ser menos parsimonioso, sigue siendo estable.
Figura 10. Prueba de racíces unitarias ARIMA(5,1,21) - círculo unitario
Fuente: elaboración propia con salida de R.
Se presentan los pronósticos obtenidos por ambos modelos.
Figura 11. Pronóstico a 20 días de VWO con ARIMA(4,1,0) 
Fuente: elaboración propia con salida de R.
Los problemas de una mala especificación o el no corregir los problemas de autocorrelación del modelo, implica que no se obtengan resultados confiables. La propuesta del ARIMA(5,1,21) mejora significativamente el pronóstico.
Figura 12. Pronóstico a 20 días de VWO con ARIMA(5,1,21) 
Fuente: elaboración propia con salida de R.
Los pronósticos para VWO los días 23 y 24 de marzo son:
| Fecha | Dato real | Pronosticado ARIMA (4,1,0) | Pronosticado ARIMA (5,1,21) |
|---|---|---|---|
| 23-mar-20 | 74.17 | 69.4100 | 70.0552 |
| 24-mar-20 | 57.57 | 70.3100 | 71.8474 |
| Criterio de información | AIC | 8194.642 | 8146.5 |
Finalmente, el Criterio de Información de Akaike muestra un mejor ajuste para el ARIMA(5,1,21).
Si bien el modelo ARIMA(5,1,21) presenta una mejor aproximación al precio real inmediato pronosticado (23 de marzo), la volatilidad del precio de EDZ provocó que el 24 de marzo cerrará en los 57usd (caída del más del 20%).
¿Cómo ajustar el pronóstico para instrumentos volátiles? ¿Y el pronóstico del S&P500?
Se presentará en un siguiente trabajo con otro tipo de modelación.
Conclusiones
En este trabajo se analizó el comportamiento de VWO, se comparó contra el S&P500 revisando su comportamiento a niveles y en rendimientos. Posteriormente, se hicieron histogramas y gráficos Q-Q que permitieron visualizar la distribución que siguen las series y la mayor parte de la concentración tanto en precios como en rendimientos.
Después, se realizaron pruebas de raíces unitarias para identificar la estacionariedad de las series en donde los resultados indicaros que las series, para que cumplan con este supuesto (al menos en media o un sentido débil), deben de ser integradas de orden I, es decir, se les tiene que aplicar una primera diferencia.
Consecutivamente, se obtuvieron los correlogramas para identificar los procesos de autocorrelación de las series, en el caso de VWO, se plantearon dos modelos, un ARIMA(4,1,0) propuesto por R y un ARIMA(5,1,21). La primera especificación, a pesar de que era un modelo estable, no corregía en su totalidad los problemas de autocorrelación provocando distorsiones en los pronósticos en tanto que el un ARIMA(5,1,21) muestra mejores resultados, corrige autocorrelación y mejora los pronósticos.
Si bien esta propuesta mejora la implementada por R, no significa que este sea el mejor modelo o que no se puedan hacer otras especificaciones, pero cumple con los supuestos y permite que tengamos más herramientas para la toma de decisiones.
Finalmente, con base al pronóstico realizado, se sugiere estimar el ARIMA(5,1,21) incorporando siempre la nueva información para ajustar el modelo. Con una perspectiva no tan clara respecto a las economías emergentes, la incertidumbre que se vive en el mercado con la propagación del COVID-19 aunada a las tensiones del sector petrolero, se recomienda una posición de hold o mantener.
Pero, si los estímulos fiscales de Europa, Asia y América Latina no tienen el efecto esperado de alentar a los agentes económicos a través de una inyección de mayor liquidez y la economía china no muestra signos de recuperación a pesar de reducir considerablemente el número de contagio del COVID-19, así como las economías emergentes, entonces se recomienda comprar VWO.
Referencias
[1] Vanguard. (2020). https://investor.vanguard.com/etf/profile/VWO