poster-hmm
0 likes102 views
This document presents a comparison of Gaussian and non-Gaussian stochastic volatility models for modeling financial asset returns. It estimates parameters for these models using a hidden Markov model approach on index fund daily return data from 2006 to 2016. The results show that non-Gaussian models generally perform better in terms of goodness-of-fit measures. Specifically, indexes for stocks, emerging markets and the Pacific performed better with a non-Gaussian assumption, while a bond index was nearly normally distributed. The document also discusses model specifications and concludes it would be interesting to relax independence assumptions between error terms.
![Hidden Markov Model for Stochastic Volatility
Vasin Suntayodom, Mengyuan Wu
University of Massachusetts, Amherst
Introduction
Recently, Bayesian estimation of stochastic volatility models via Markov chain Monte Carlo (MCMC) have gained some popularity because Bayesian
estimators have become much easier to compute. However, the purpose of Hidden Markov model for stochastic volatility project is
to show that it provides alternative approach to estimate parameters in stochastic volatility model.
to compare the goodness of 鍖t between Gaussian Stochastic Volatility model and Non-Gaussian Stochastic Volatility
The key idea is the use of iterated numerical integration. This method involves an approximation to the stochastic volatility likelihood that can be made
arbitrarily accurate.
Numerical Results
The Gaussian stochastic volatility model with-
out leverage was 鍖tted to the daily continuously
compounded returns on the index funds for the
period from the end of 1 May 2006 to the end of
27 May 2016. The starting values are 袖0 = 0.7,
0 = 0.95 and 0 = 0.12. Using m = 200 and gt
values from -2.5 to 2.5, we obtain the result as
follows :
硫 AIC BIC
v鍖nx 0.01490 0.99755 0.19384 -16219.70590 -16202.18969
veurx 0.01535 0.99254 0.14358 -14976.17353 -14958.65732
veiex 0.02027 0.99785 0.15170 -14976.24524 -14958.72903
vbltx 0.00616 0.99167 0.07627 -18365.26568 -18347.74947
vbisx 0.00205 0.99905 0.06691 -26627.52034 -26610.00413
vpacx 0.01316 0.98901 0.15359 -15761.37733 -15743.86112
The result obtained in table above lead us to
conclude that the parameter estimated obtained
by HMM stabilizes for m somewhere between
100 and 200. Next we examine Non-Gaussian
Stochastic Volatility with gt values from -7.5 to
7.5, we obtain the result as follows:
硫 僚 AIC BIC
v鍖nx 0.00859 0.98457 0.17785 12.91058 -16228.02737 -16204.67242
veurx 0.01118 0.98899 0.13716 14.01139 -14989.97749 -14966.62254
veiex 0.01162 0.98766 0.14411 27.46222 -14981.03284 -14957.67789
vbltx 0.00616 0.99168 0.07627 31058.90450 -18363.26517 -18339.91022
vbisx 0.00106 0.99782 0.04748 7.73066 -26688.34194 -26664.98699
vpacx 0.00990 0.98141 0.16345 28.24942 -15772.48712 -15749.13217
According to the table above, S&P 500 index,
European stock index, Emerging markets index,
Short term bond index and Paci鍖c stock index
under Non-Gaussian Stochastic Volatility model
perform better than the Gaussian Stochastic
Volatility.
On the contrary, the long term bond index is
almost normally distributed and the degree of
freedom we got also support this evidence (僚 =
31058.90450). Both AIC and BIC are almost
identical in two models, therefore the models
are indistinguishable. In the case of Short term
bond index, kernal density also indicates that it
deviates from normality assumption.
As a result, the Non-Gaussian stochastic volatil-
ity model performs much better than Gaussian
stochastic volatility model.
Data Descriptions
We will use various major index funds obtained
from Bloomberg from May 2006 until May 2016
on daily basis. The following is the brief descrip-
tion of the data :
1. S&P 500 index (v鍖nx)
2. European stock index (veurx)
3. Emerging markets stock index (veiex)
4. Long term bond fund (vbltx)
5. Short term bond fund (vbisx)
6. Paci鍖c stock index (vpacx)
The data set covers daily closing price data from
the end of 1 May 2006 to the end of 27 May 2016.
Model Speci鍖cation
g1 g2 g3 gt gT
r1 r2 r3 rt rT
. . . . . .
The 鍖rst form of the model which we consider
here and the best known, is Gaussian Stochas-
tic Volatility or Heston Model where the asset
returns rt on the observation equation satisfy
rt = t硫 exp (gt/2)
gt+1 = gt + 侶t,
where || < 1 and { t} is Gaussian white noise
sequence with mean 0 and variance 1, {侶t} is also
Gaussian white noise sequence with mean 0 and
variance 2
僚. Suppose { t} and {僚t} are inde-
pendent for Gaussian stochastic volatility model
without leverage e鍖ect. Then we improve the ob-
servation equation in the model by
rt = t (硫 exp (gt/2) + 両) .
The additional parameter 両 0 is persuasive on
the fact that some baseline volatility is always
presented. We also relax normality assumption
and assume that t has t-distribution with 僚 de-
grees of freedom. The basic model which assume
standard normal distribution for t is a special
case when 僚 .
Conclusion
In this work, we introduced the problem of volatility estimation in prices of 鍖nancial assets. We have
discussed why this is an interesting endeavor to use stochastic volatility models and have introduced
a popular stochastic volatility model to capture some of the statistical features found in real-world
data. It would be interesting to experiment by relaxing assumption that { t} and {僚t} have to
be independent for Gaussian stochastic volatility model with leverage e鍖ect. Speci鍖cally, for all t
t
侶t
N (0, 裡) with 裡 =
1
2 .
Reference
[1] J. Hull and A. White. The pricing of options on assets with stochastic volatility the Journal of Finance, pp. 42(2):281-300, 1987.
[2] E. Jacquier, N. G. Polson, and P. E. Rossi. Bayesian analysis of stochastic volatility models with fat-tails and correlated errors Journal of Econometric,
122(1):185-212, 2004.
[3] P.Glasserman Monte Carlo methods in 鍖nancial engineering (2004), Volume 53, Springer
[4] Bartolucci, F.,De Luca, G. Maximum likelihood estimation of a latent variable time-series model.Applied Stochastic Models in Business and Industry, pp. 17:
5-17, 2001
[5] Fridman, M., Harris, L. A maximum likelihood approach for non-Gaussian stochastic volatility models.Journal of Business and Economic Statistics, pp. 16:
284-291, 1998
Acknowledgment
We wish to thank particularly Panit Arunanondchai and Rene Cabrera for his proof reading and very useful comments during the preparation of this paper.](https://image.slidesharecdn.com/ed26e66b-e6e0-46e4-a6be-fac84d7884a0-170104154329/85/poster-hmm-1-320.jpg)
poster-hmm
- 1. Hidden Markov Model for Stochastic Volatility Vasin Suntayodom, Mengyuan Wu University of Massachusetts, Amherst Introduction Recently, Bayesian estimation of stochastic volatility models via Markov chain Monte Carlo (MCMC) have gained some popularity because Bayesian estimators have become much easier to compute. However, the purpose of Hidden Markov model for stochastic volatility project is to show that it provides alternative approach to estimate parameters in stochastic volatility model. to compare the goodness of 鍖t between Gaussian Stochastic Volatility model and Non-Gaussian Stochastic Volatility The key idea is the use of iterated numerical integration. This method involves an approximation to the stochastic volatility likelihood that can be made arbitrarily accurate. Numerical Results The Gaussian stochastic volatility model with- out leverage was 鍖tted to the daily continuously compounded returns on the index funds for the period from the end of 1 May 2006 to the end of 27 May 2016. The starting values are 袖0 = 0.7, 0 = 0.95 and 0 = 0.12. Using m = 200 and gt values from -2.5 to 2.5, we obtain the result as follows : 硫 AIC BIC v鍖nx 0.01490 0.99755 0.19384 -16219.70590 -16202.18969 veurx 0.01535 0.99254 0.14358 -14976.17353 -14958.65732 veiex 0.02027 0.99785 0.15170 -14976.24524 -14958.72903 vbltx 0.00616 0.99167 0.07627 -18365.26568 -18347.74947 vbisx 0.00205 0.99905 0.06691 -26627.52034 -26610.00413 vpacx 0.01316 0.98901 0.15359 -15761.37733 -15743.86112 The result obtained in table above lead us to conclude that the parameter estimated obtained by HMM stabilizes for m somewhere between 100 and 200. Next we examine Non-Gaussian Stochastic Volatility with gt values from -7.5 to 7.5, we obtain the result as follows: 硫 僚 AIC BIC v鍖nx 0.00859 0.98457 0.17785 12.91058 -16228.02737 -16204.67242 veurx 0.01118 0.98899 0.13716 14.01139 -14989.97749 -14966.62254 veiex 0.01162 0.98766 0.14411 27.46222 -14981.03284 -14957.67789 vbltx 0.00616 0.99168 0.07627 31058.90450 -18363.26517 -18339.91022 vbisx 0.00106 0.99782 0.04748 7.73066 -26688.34194 -26664.98699 vpacx 0.00990 0.98141 0.16345 28.24942 -15772.48712 -15749.13217 According to the table above, S&P 500 index, European stock index, Emerging markets index, Short term bond index and Paci鍖c stock index under Non-Gaussian Stochastic Volatility model perform better than the Gaussian Stochastic Volatility. On the contrary, the long term bond index is almost normally distributed and the degree of freedom we got also support this evidence (僚 = 31058.90450). Both AIC and BIC are almost identical in two models, therefore the models are indistinguishable. In the case of Short term bond index, kernal density also indicates that it deviates from normality assumption. As a result, the Non-Gaussian stochastic volatil- ity model performs much better than Gaussian stochastic volatility model. Data Descriptions We will use various major index funds obtained from Bloomberg from May 2006 until May 2016 on daily basis. The following is the brief descrip- tion of the data : 1. S&P 500 index (v鍖nx) 2. European stock index (veurx) 3. Emerging markets stock index (veiex) 4. Long term bond fund (vbltx) 5. Short term bond fund (vbisx) 6. Paci鍖c stock index (vpacx) The data set covers daily closing price data from the end of 1 May 2006 to the end of 27 May 2016. Model Speci鍖cation g1 g2 g3 gt gT r1 r2 r3 rt rT . . . . . . The 鍖rst form of the model which we consider here and the best known, is Gaussian Stochas- tic Volatility or Heston Model where the asset returns rt on the observation equation satisfy rt = t硫 exp (gt/2) gt+1 = gt + 侶t, where || < 1 and { t} is Gaussian white noise sequence with mean 0 and variance 1, {侶t} is also Gaussian white noise sequence with mean 0 and variance 2 僚. Suppose { t} and {僚t} are inde- pendent for Gaussian stochastic volatility model without leverage e鍖ect. Then we improve the ob- servation equation in the model by rt = t (硫 exp (gt/2) + 両) . The additional parameter 両 0 is persuasive on the fact that some baseline volatility is always presented. We also relax normality assumption and assume that t has t-distribution with 僚 de- grees of freedom. The basic model which assume standard normal distribution for t is a special case when 僚 . Conclusion In this work, we introduced the problem of volatility estimation in prices of 鍖nancial assets. We have discussed why this is an interesting endeavor to use stochastic volatility models and have introduced a popular stochastic volatility model to capture some of the statistical features found in real-world data. It would be interesting to experiment by relaxing assumption that { t} and {僚t} have to be independent for Gaussian stochastic volatility model with leverage e鍖ect. Speci鍖cally, for all t t 侶t N (0, 裡) with 裡 = 1 2 . Reference [1] J. Hull and A. White. The pricing of options on assets with stochastic volatility the Journal of Finance, pp. 42(2):281-300, 1987. [2] E. Jacquier, N. G. Polson, and P. E. Rossi. Bayesian analysis of stochastic volatility models with fat-tails and correlated errors Journal of Econometric, 122(1):185-212, 2004. [3] P.Glasserman Monte Carlo methods in 鍖nancial engineering (2004), Volume 53, Springer [4] Bartolucci, F.,De Luca, G. Maximum likelihood estimation of a latent variable time-series model.Applied Stochastic Models in Business and Industry, pp. 17: 5-17, 2001 [5] Fridman, M., Harris, L. A maximum likelihood approach for non-Gaussian stochastic volatility models.Journal of Business and Economic Statistics, pp. 16: 284-291, 1998 Acknowledgment We wish to thank particularly Panit Arunanondchai and Rene Cabrera for his proof reading and very useful comments during the preparation of this paper.