Financial Time Series Volatility Analysis Using Gaussian Process State-Space Models
|
|
- Melissa Deirdre Russell
- 5 years ago
- Views:
Transcription
1 15 IEEE Global Conference on Signal and Information Processing (GlobalSIP) Financial Time Series Volatility Analysis Using Gaussian Process State-Space Models Jianan Han, Xiao-Ping Zhang Department of Electrical and Computer Engineering, Ryerson University 35 Victoria Street, Toronto, ON, Canada M5B K3 jhan, Abstract In this paper, we propose a novel nonparametric modeling framework for financial time series data analysis, and apply it to the problem of time varying volatility modeling. Existing parametric models have a rigid-form transition function and they often have over-fitting problems when model parameters are estimated using maximum likelihood methods. These drawbacks effect the models prediction performance. To solve this problem, we take Bayesian nonparametric modeling approach. By adding Gaussian process prior to the hidden state transition process, we extend the standard state-space model to a Gaussian process state-space model. We introduce the Gaussian process regression stochastic volatility (GPRSV) model and instead of using maximum likelihood methods, we use Monte Carlo inference algorithms. We demonstrate performance of our model and inference methods with both simulated and empirical financial data. I. INTRODUCTION The problem of analysing financial time series data is an important task for both research and investment. In the last decades, many researchers took the modeling approach to describe financial data. Modeling give us a way of discovering knowledge from data, and making prediction. From this point, modeling financial time series data is very similar to modeling signals in engineering application. For example, in presence of noise, filtering methods like Kalman filter and particle filter can be applied to financial data too. With the recent development of Bayesian nonparametric modeling in signal processing community, we can model financial data with more flexible tools, like Gaussian Process and Copula Process. In this work, we apply Gaussian process regression to the transition function of volatility modeling. For inference methods, we apply Bayesian inference framework to compute the jointly distribution. Compared with the traditional parametric models, our new method is more flexible and free of over fitting problem. Volatility can be expressed as the standard deviation of an asset return and it is widely used to describe the variability of financial time series data. In some degree, volatility is a measurement of the financial time series magnitude and fluctuation s speed. Volatility modeling has been one of the most active financial time series research areas in the last decade, and it is of great importance for both finance market practitioners and academic researchers. Volatility is a forwardlooking concept, we often model the financial time series return variance conditioned on all the relevant information I t 1 which is defined as follow: σ t = var(r t I t 1 )=E((r t μ t ) I t 1 ), (1) where r t is the asset return and μ t is expected value of r t. II. REVIEW OF EXISTING VOLATILITY MODELS There are two main classes of modeling time changing variance: the generalized autoregressive conditional heteroskedasticity (GARCH) and the stochastic volatility (SV) model. For the GARCH typed models, ARCH [1] was first introduced in 198 by Engle. GARCH [] was generalized form of ARCH. The other type is Stochastic Volatility models, see [3]. There are many variants for both these two types of models. Numerous researchers developed extensions of these models applying for empirical analyses in financial economics and macroeconomics. A lot of variants and extensions have been studied. A. GARCH models Autoregressive Conditional Heteroscedastictiy (ARCH) model was first introduced by Engle(198). Bollerslev extended the model, and give a form of the Generalised Autoregressive Conditional Heteroscedastictiy models. Standard GARCH(1,1) model assumes the return of the asset follows a Gaussian distribution. In this work, we simplify its mean μ to be zero and the variance of the distribution is time-varying: x t N(,σt ) σt = α + α 1 x t 1 + βσt 1 (a) (b) Some GARCH extensions are reviewed here too. The Exponential GARCH (EGARCH) [4] can be summarized as: log(σt )=α + α 1 g(x t 1)+βlog(αt 1) g(x t )=θx t + λ x t (3a) (3b) Another popular GARCH extension is the GJR-GARCH model [5]: σ t = α + βα t 1 + γx t 1I t 1 (4a) /15/$ IEEE 358
2 15 IEEE Global Conference on Signal and Information Processing (GlobalSIP) {, if x t 1 I t 1 = (4b) 1, if x t 1 < Parameters of GARCH class models can be learned using maximum likelihood methods. There are a lot of papers on the topic and many software environment provide the program implementing the algorithms. Among these, see Kevin Shepperd s Matlab code of Oxford MFE Toolbox. B. Stochastic Volatility models The first discrete time-varying stochastic volatility model was introduced by Taylor [3], [6]. The logarithm of variance was modeled by a latent AR(1) process. Taylor s stochastic model can be presented as: r t = μ t + a t = μ t + σ t ɛ t (5a) log(σt )=α + α 1 log(σt 1)+σ n η t (5b) where α 1 is a parameter which controls the persistence of logarithm variance and the value of α 1 is between (-1,1). There are two independent and identically distributed random variables ɛ t and η t. The original idea of SV model assume these two noise parts to be i.i.d. normal distributed. The inference of a SV model parameters is not as straightforward as the corresponding simple GARCH typed model. In [7], Shephard reviewed SV models and inference methods like methods of moments (MM) and quasi-maximum likelihood (QML). Simulation-based methods to learn SV models become more and more popular because of their accuracy and flexibility of handling complicated models. III. GAUSSIAN PROCESS STATE-SPACE MODELS A. Gaussian Process A Gaussian Process can be viewed as an extension of multivariate Gaussian distribution to infinite dimensions. Any finite subset of sampled from the process follows a multiple Gaussian distribution. Also, Gaussian Process can be considered as a normal distribution over function, and it is determined by mean function and covariance function. f(x) =GP(m(x),k(x, x )) (6) All values of f(x) at any location x are jointly Gaussian distributed. B. State-Space Models State-space model (SSMs) are one of the most widely used models for effectively modeling time series and describing dynamical system.in Finance area, State-space models can generalize other popular time series models such as ARMA, (G) ARCH and SV models as we mentioned before. A State- Space Model can be a general framework for descripting dynamical phenomena. The model consists of two parts: an hidden state x t and observation variable y t. The general form of State-Space model can be summarized as: x t = f(x t 1 )+ɛ, x t R M (7a) y t = g(x t )+ν, y t R D (7b) where, ɛ and are normal distributed noise. The unknown function f describe the system dynamics and function g links the observation and the system hidden states. Both the f and g functions can be linear or non-lienar. We can combine Gaussian process and state-space model together. The way of combining the two is to use the statespace model s structure and apply Gaussian process to describe the hidden state transition function. The essence of the GP- SSMs is to change the rigid form of states transition function of traditional state-space models with a Gaussian process prior. Financial data exhibits many dynamics because the market is changing all the time and a lot of small change of the involved factors can result in significant fluctuation. As more and more data is available to process, the rigid form of state transition function in traditional state-space models becomes the bottleneck of improving the models forecast performance. We assume the hidden state transition function f is sampled from a Gaussian process, we extend the standard state-space model to a Gaussian process state-space model (GP-SSM). Compared with standard SSM, GP-SSM is a more flexible and powerful tool to model time series data. The linear SSMs with Gaussian noise can be inferenced using Kalman Filter [8]. But linear Gaussian SSMs can only model a limited set of phenomena. To learn a GP-SSM is more difficult than a standard SSM. In our work, we take the Bayesian approach to solve the GP-SSMs inference problem. Bayesian filtering is a type of technique used to estimate the hidden states in dynamic systems. The goal of Bayesian filtering is to recursively compute the posterior distribution of the current hidden state given the whole history of observations. IV. GAUSSIAN PROCESS REGRESSION STOCHASTIC VOLATILITY MODEL Our Gaussian process regression stochastic volatility (GPRSV) model can be viewed as an instant of GP-SSM. We treat volatility as unobserved latent variable, like the standard stochastic volatility models. Instead of assuming that the latent volatility follows a p-step autoregressive process. We place a Gaussian process prior to the unknown transition function f on the volatility process. A GPRSV model can be presented with the following equations: a t = r t μ = σ t ɛ t v t = log(σt )=f(v t 1 )+τη t f GP(m(x),k(x, x )) (8a) (8b) (8c) where r t is the asset return at time t and μ is the mean of asset return series, and a t is the innovation of the return series. v t is the logarithm of variance at time t, ɛ t and η t are i.i.d. Gaussian (or student s t) distributed noises. τ is unknown scaling parameters to be estimated. The function f is the hidden state transition function. Here we assume this function f follows a Gaussian process, which is defined by the mean function m(x) and covariance function k(x, x ). We show the graphical model representation of GPRSV model in Fig
3 15 IEEE Global Conference on Signal and Information Processing (GlobalSIP) f 1 f f * mdoel with RAPCF algorithm. For particle filter review, see [13] and [14]. v v 1 v a a 1 a... Fig. 1. Graphical model representation of a Gaussian process regression stochastic volatility model, where a t is the observation variable at time t, and v t are the hidden variable (logarithm of volatility) at time t. f t is the Gaussian process sampled function value at time t, and the thick horizontal line represent fully connected nodes. Hyper-parameters of the Gaussian process are omitted in the figure. In GPRSV model we use logarithm of variance instead of standard deviation directly in our model. The advantage of using logarithm form is explained in [3] and [4]. We build a GPRSV model with the following four steps: Step 1: specify mean equation. First we need to test the serial dependence in the return series. If the series are linear dependent, we should use an econometric model to remove the linear dependence in return series. Step : test ARCH effect. The residuals of the asset return a t expressed in (8a) are often used to test the series conditional heteroscedasticity. This conditional heteroscedasticity is also known as the ARCH effects [9]. Step 3: specify volatility equation. We need to specify both the mean and covariance functions for Gaussian process. Besides these functions, the initial value of hyperparameters associated with them need to be specified as well. Step 4: estimate model parameters and check model fitness. We use training data to estimate the unknown parameters. Once we get our estimated parameters we can use testing data to test learned model, and it is necessary to check the fitness of model. V. INFERENCE FOR THE GPRSV MODEL It is much more challenging to learn a GPRSV model than its parametric competitors. We target at jointly learning the hidden states trajectory and the unknown Gaussian process regression function values and hyper-parameters. Our approach is marginalizing out the Gaussian process regression function values, and then jointly learning the hidden volatility states and hyper-parameters using Monte Carlo methods. Most of previous work on inference GP-SSMs focused on filtering and smoothing hidden variables without jointly learning the Gaussian process hyper-parameters. In [1], a novel particle Markov chain Monte Carlo (particle MCMC) algorithm, particle Gibbs with ancestor sampling (PGAS) was proposed. In [11], Frigola et al. apply the algorithm to learn a GP-SSM s hidden states and Gaussian process dynamics jointly. In [1], a regularized auxiliary particle filter which the authors named as Regularized Auxiliary Particle Chain Filter (RAPCF) was introduced. In Algorithm 1 we show how to learn a GPRSV v * a * Algorithm 1 RAPCF for GPRSV Model 1: Input: return data r 1:T, number of particles N, shrinkage parameter <λ<1, prior p(θ). : Remove linear dependence from r 1:T to get the residuals a 1:T. 3: Sample N parameter particles from the prior, and set initial importance weights, W i =1/N 4: for t =1to T do 5: Shrink parameter particles towards empirical means θ t 1 =Σ N i=1wt 1θ i t 1 i (9a) θ t i = λθt 1 i +(1 λ) θ t 1 (9b) 6: Compute expected states: 7: Compute important weights μ i t = E(v t θ i t,v i 1:t 1) (1) g i t W i t 1p(a t μ i t, θ i t) (11) 8: 9: Resample N auxiliary indices {j} according to {gt}. i Propagate the chains of v t forward, {v j 1:t 1 } j J. 1: Add jitter: θt 1 i N(θ j t, (1 λ )V t 1 ), and V t 1 is empirical covariance of θ t 1. 11: Propose new states v j t p(v t θ j t,v j 1:t 1,a 1:t 1) 1: Adjust weights with newly proposed states: W i t p(a t v j t, θ j t )/p(a t μ j t, θ j t ) (1) 13: end for 14: Output: particles of v j 1:T, particles of θj t and particle weights W j t. VI. EXPERIMENTS We apply both the simulated and empirical financial data to demonstrate performance of our GPRSV model and inference method. For simulated data, we generated one synthetic dataset of length T = according to equation 8. We sample our hidden state transition function f from a Gaussian process prior. We specified that the mean function m(x t ) and the covariance function k(y, z) as follow: m(x t )=ax t 1 k(y, z) =γexp(.5 y z /l ) (13a) (13b) We show the learned hidden states density for the simulated data in Fig. and the predictive Log-likelihood value compared with the true value in Fig. 3. For real financial data, we chose IBM stock daily closing price data. The data period is from January 1, 1988 to September 14, 3. We compare our GPRSV model with standard parametric volatility models: GARCH, GJR-GARCH and standard SV. For GPRSV model, the Gaussian process dynamics are specified as follows: the mean function is m(x t )=ax t 1 36
4 15 IEEE Global Conference on Signal and Information Processing (GlobalSIP) Evolution of the state density of the most widely-used loss functions in financial research literature. In our experiment, we adopt a class of statistical loss functions instead of a particular one. We take the following loss functions: MAD : L(ˆσ t+m,σ t+m )=n 1 ˆσ t+m σ t+m (14) Fig.. Estimated hidden states densities of simulated Data. There are iteration steps for the simulated data, and we plot every 5 densities in this figure. The distribution densities are generated using particles and weights as described in Algorithm 1. Value Predictive Log-likelihood Values True Value GPRSV MLAE : L(ˆσ t+m,σ t+m )=n 1 QLIKE : L(ˆσ t+m,σ t+m) =n 1 HMSE : L(ˆσ t+m,σ t+m) =n 1 n log( ˆσ t+m σ t+m ) (15) (ˆσ t+m/σ t+m+logσ t+m) (ˆσ t+m/σ t+m 1) (16) (17) In our experiment, we apply the high frequency data to calculate the realized volatility [16] and make it as the proxy for true volatility. Compared with the squared return, realized volatility is considered to be more precise proxy for volatility forecast evaluation. The loss function result is presented in Table II. TABLE II RESULTS OF IBM VOLATILITY FORECAST USING LOSS FUNCTIONS Time(t) Fig. 3. RAPCF algorithm learned predictive Log-likelihood value are compared with true value calculate from equation 8. We discard the first 5 burn in iterations. and the covariance function is the squared exponential covariance function k(y, z) =γexp(.5 y z /l ). The hyperparameters include a, γ, l and likelihood function parameter log(sn). The learned parameters are presented in Table I. There is many metrics to evaluate different forecast models. One of the most popular approach is using a particular statistical loss function, the model which achieved a minimized loss function value is the best forecasting model [15]. There are extensive choices of lost functions. QLIKE function is one TABLE I ESTIMATED GPRSV MODEL HYPER-PARAMETERS RESULTS FOR IBM DAILY RETURN DATA a γ l log(sn) Model MAD MLAE HMSE QLIKE GARCH GJR-GARCH SV GPRSV Note: The lowest loss function value is marked using the bold text. GPRSV model obtained the lowest loss function values for all loss functions. The volatility proxy is the 65-minutes sampled realized volatility. VII. CONCLUSIONS In this paper, we proposed a Gaussian process regression based volatility model for the problem of predicting the time varying volatility of financial time series data. After we introduced the GPRSV model, we gave a solution to jointly learning the hidden volatility states and the Gaussian process dynamics. Based on our experiment results, we can successfully learn the Gaussian process regression stochastic volatility model s hidden states and hyper-parameters. Also the Gaussian process regression based stochastic volatility models can achieve better performance compared with the standard economic parametric models. 361
5 15 IEEE Global Conference on Signal and Information Processing (GlobalSIP) REFERENCES [1] R. F. Engle, Autoregressive conditional heteroscedasticity with estimates of the variance of united kingdom inflation, Econometrica: Journal of the Econometric Society, pp , 198. [] T. Bollerslev, Generalized autoregressive conditional heteroskedasticity, Journal of Econometrics, vol. 31, pp , [3] S. J. Taylor, Modelling financial time series, [4] D. Nelson, Conditional Heteroskedasticity in Asset Returns: A New Approach, Econometrica, vol. 59, no., pp , March [5] L. Glosten, R. Jagannathan, and D. Runkle, On the relation between the expected value and the volatility of the nominal excess return on stocks, The Journal of Finance, vol. 48, no. 5, pp , [6] S. J. Taylor, Financial returns modelled by the product of two stochastic processes, a study of daily sugar prices. Oxford University Press, 5. [7] S. Neil and T. Andersen, Stochastic Volatility: Origins and Overview, University of Oxford, Department of Economics, Economics Series Working Papers 389, 8. [8] R. E. Kalman, A new approach to linear filtering and prediction problems, Transactions of the ASME Journal of Basic Engineering, vol. 8, no. Series D, pp , 196. [9] R. Tsay, Analysis of financial time series. Wiley, 1. [1] F. Lindsten, M. Jordan, and T. Schön, Particle gibbs with ancestor sampling, Journal of Machine Learning Research, vol. 15, pp , 14. [11] R. Frigola, F. Lindsten, T. B. Schön, and C. E. Rasmussen, Bayesian inference and learning in gaussian process state-space models with particle MCMC, in Advances in Neural Information Processing Systems 6, 13, pp [1] Y. Wu, J. M. Hernández-Lobato, and Z. Ghahramani, Gaussian process volatility model, in Advances in Neural Information Processing Systems, 14, pp [13] M. S. Arulampalam, S. Maskell, N. Gordon, and T. Clapp, A tutorial on particle filters for online nonlinear/non-gaussian bayesian tracking, Signal Processing, IEEE Transactions on, vol. 5, no., pp ,. [14] C. Hue, J.-P. Le Cadre, and P. Perez, Sequential monte carlo methods for multiple target tracking and data fusion, Signal Processing, IEEE Transactions on, vol. 5, no., pp , Feb. [15] C. Brownlees, R. Engle, and B. Kelly, A practical guide to volatility forecasting through calm and storm, Journal of Risk, vol. 14(), pp. 1, 11. [16] A. Torben, T. Bollerslev, F. Diebold, and P. Labys, Modeling and forecasting realized volatility, Econometrica, vol. 71, no., pp , 3. 36
Chapter 6 Forecasting Volatility using Stochastic Volatility Model
Chapter 6 Forecasting Volatility using Stochastic Volatility Model Chapter 6 Forecasting Volatility using SV Model In this chapter, the empirical performance of GARCH(1,1), GARCH-KF and SV models from
More informationUniversity of Toronto Financial Econometrics, ECO2411. Course Outline
University of Toronto Financial Econometrics, ECO2411 Course Outline John M. Maheu 2006 Office: 5024 (100 St. George St.), K244 (UTM) Office Hours: T2-4, or by appointment Phone: 416-978-1495 (100 St.
More informationConditional Heteroscedasticity
1 Conditional Heteroscedasticity May 30, 2010 Junhui Qian 1 Introduction ARMA(p,q) models dictate that the conditional mean of a time series depends on past observations of the time series and the past
More informationMachine Learning for Quantitative Finance
Machine Learning for Quantitative Finance Fast derivative pricing Sofie Reyners Joint work with Jan De Spiegeleer, Dilip Madan and Wim Schoutens Derivative pricing is time-consuming... Vanilla option pricing
More informationState Switching in US Equity Index Returns based on SETAR Model with Kalman Filter Tracking
State Switching in US Equity Index Returns based on SETAR Model with Kalman Filter Tracking Timothy Little, Xiao-Ping Zhang Dept. of Electrical and Computer Engineering Ryerson University 350 Victoria
More informationIndian Institute of Management Calcutta. Working Paper Series. WPS No. 797 March Implied Volatility and Predictability of GARCH Models
Indian Institute of Management Calcutta Working Paper Series WPS No. 797 March 2017 Implied Volatility and Predictability of GARCH Models Vivek Rajvanshi Assistant Professor, Indian Institute of Management
More informationResearch Article The Volatility of the Index of Shanghai Stock Market Research Based on ARCH and Its Extended Forms
Discrete Dynamics in Nature and Society Volume 2009, Article ID 743685, 9 pages doi:10.1155/2009/743685 Research Article The Volatility of the Index of Shanghai Stock Market Research Based on ARCH and
More informationStatistical Inference and Methods
Department of Mathematics Imperial College London d.stephens@imperial.ac.uk http://stats.ma.ic.ac.uk/ das01/ 14th February 2006 Part VII Session 7: Volatility Modelling Session 7: Volatility Modelling
More informationCourse information FN3142 Quantitative finance
Course information 015 16 FN314 Quantitative finance This course is aimed at students interested in obtaining a thorough grounding in market finance and related empirical methods. Prerequisite If taken
More informationLecture Note 9 of Bus 41914, Spring Multivariate Volatility Models ChicagoBooth
Lecture Note 9 of Bus 41914, Spring 2017. Multivariate Volatility Models ChicagoBooth Reference: Chapter 7 of the textbook Estimation: use the MTS package with commands: EWMAvol, marchtest, BEKK11, dccpre,
More informationARCH and GARCH models
ARCH and GARCH models Fulvio Corsi SNS Pisa 5 Dic 2011 Fulvio Corsi ARCH and () GARCH models SNS Pisa 5 Dic 2011 1 / 21 Asset prices S&P 500 index from 1982 to 2009 1600 1400 1200 1000 800 600 400 200
More informationStatistical Models and Methods for Financial Markets
Tze Leung Lai/ Haipeng Xing Statistical Models and Methods for Financial Markets B 374756 4Q Springer Preface \ vii Part I Basic Statistical Methods and Financial Applications 1 Linear Regression Models
More informationVolatility Models and Their Applications
HANDBOOK OF Volatility Models and Their Applications Edited by Luc BAUWENS CHRISTIAN HAFNER SEBASTIEN LAURENT WILEY A John Wiley & Sons, Inc., Publication PREFACE CONTRIBUTORS XVII XIX [JQ VOLATILITY MODELS
More informationUsing MCMC and particle filters to forecast stochastic volatility and jumps in financial time series
Using MCMC and particle filters to forecast stochastic volatility and jumps in financial time series Ing. Milan Fičura DYME (Dynamical Methods in Economics) University of Economics, Prague 15.6.2016 Outline
More informationFinancial Econometrics
Financial Econometrics Volatility Gerald P. Dwyer Trinity College, Dublin January 2013 GPD (TCD) Volatility 01/13 1 / 37 Squared log returns for CRSP daily GPD (TCD) Volatility 01/13 2 / 37 Absolute value
More informationIdiosyncratic risk, insurance, and aggregate consumption dynamics: a likelihood perspective
Idiosyncratic risk, insurance, and aggregate consumption dynamics: a likelihood perspective Alisdair McKay Boston University June 2013 Microeconomic evidence on insurance - Consumption responds to idiosyncratic
More informationShort-selling constraints and stock-return volatility: empirical evidence from the German stock market
Short-selling constraints and stock-return volatility: empirical evidence from the German stock market Martin Bohl, Gerrit Reher, Bernd Wilfling Westfälische Wilhelms-Universität Münster Contents 1. Introduction
More informationIntroduction to Sequential Monte Carlo Methods
Introduction to Sequential Monte Carlo Methods Arnaud Doucet NCSU, October 2008 Arnaud Doucet () Introduction to SMC NCSU, October 2008 1 / 36 Preliminary Remarks Sequential Monte Carlo (SMC) are a set
More informationDynamic conditional score volatility models Szabolcs Blazsek GESG seminar 30 January 2015 Universidad Francisco Marroquín, Guatemala
Dynamic conditional score volatility models Szabolcs Blazsek GESG seminar 30 January 2015 Universidad Francisco Marroquín, Guatemala From GARCH(1,1) to Dynamic Conditional Score volatility models GESG
More informationTime series: Variance modelling
Time series: Variance modelling Bernt Arne Ødegaard 5 October 018 Contents 1 Motivation 1 1.1 Variance clustering.......................... 1 1. Relation to heteroskedasticity.................... 3 1.3
More informationFinancial Data Mining Using Flexible ICA-GARCH Models
55 Chapter 11 Financial Data Mining Using Flexible ICA-GARCH Models Philip L.H. Yu The University of Hong Kong, Hong Kong Edmond H.C. Wu The Hong Kong Polytechnic University, Hong Kong W.K. Li The University
More informationFE570 Financial Markets and Trading. Stevens Institute of Technology
FE570 Financial Markets and Trading Lecture 6. Volatility Models and (Ref. Joel Hasbrouck - Empirical Market Microstructure ) Steve Yang Stevens Institute of Technology 10/02/2012 Outline 1 Volatility
More informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2017, Mr. Ruey S. Tsay. Solutions to Final Exam
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2017, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (40 points) Answer briefly the following questions. 1. Describe
More informationDoes Volatility Proxy Matter in Evaluating Volatility Forecasting Models? An Empirical Study
Does Volatility Proxy Matter in Evaluating Volatility Forecasting Models? An Empirical Study Zhixin Kang 1 Rami Cooper Maysami 1 First Draft: August 2008 Abstract In this paper, by using Microsoft stock
More informationMeasuring DAX Market Risk: A Neural Network Volatility Mixture Approach
Measuring DAX Market Risk: A Neural Network Volatility Mixture Approach Kai Bartlmae, Folke A. Rauscher DaimlerChrysler AG, Research and Technology FT3/KL, P. O. Box 2360, D-8903 Ulm, Germany E mail: fkai.bartlmae,
More informationAbsolute Return Volatility. JOHN COTTER* University College Dublin
Absolute Return Volatility JOHN COTTER* University College Dublin Address for Correspondence: Dr. John Cotter, Director of the Centre for Financial Markets, Department of Banking and Finance, University
More informationOn modelling of electricity spot price
, Rüdiger Kiesel and Fred Espen Benth Institute of Energy Trading and Financial Services University of Duisburg-Essen Centre of Mathematics for Applications, University of Oslo 25. August 2010 Introduction
More informationVolatility Analysis of Nepalese Stock Market
The Journal of Nepalese Business Studies Vol. V No. 1 Dec. 008 Volatility Analysis of Nepalese Stock Market Surya Bahadur G.C. Abstract Modeling and forecasting volatility of capital markets has been important
More informationHigh-Frequency Data Analysis and Market Microstructure [Tsay (2005), chapter 5]
1 High-Frequency Data Analysis and Market Microstructure [Tsay (2005), chapter 5] High-frequency data have some unique characteristics that do not appear in lower frequencies. At this class we have: Nonsynchronous
More informationMODELING VOLATILITY OF US CONSUMER CREDIT SERIES
MODELING VOLATILITY OF US CONSUMER CREDIT SERIES Ellis Heath Harley Langdale, Jr. College of Business Administration Valdosta State University 1500 N. Patterson Street Valdosta, GA 31698 ABSTRACT Consumer
More informationTechnical Appendix: Policy Uncertainty and Aggregate Fluctuations.
Technical Appendix: Policy Uncertainty and Aggregate Fluctuations. Haroon Mumtaz Paolo Surico July 18, 2017 1 The Gibbs sampling algorithm Prior Distributions and starting values Consider the model to
More information12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006.
12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006. References for this Lecture: Robert F. Engle. Autoregressive Conditional Heteroscedasticity with Estimates of Variance
More informationFinancial Time Series Analysis (FTSA)
Financial Time Series Analysis (FTSA) Lecture 6: Conditional Heteroscedastic Models Few models are capable of generating the type of ARCH one sees in the data.... Most of these studies are best summarized
More informationModelling Returns: the CER and the CAPM
Modelling Returns: the CER and the CAPM Carlo Favero Favero () Modelling Returns: the CER and the CAPM 1 / 20 Econometric Modelling of Financial Returns Financial data are mostly observational data: they
More informationInference of the Structural Credit Risk Model
Inference of the Structural Credit Risk Model using Yuxi Li, Li Cheng and Dale Schuurmans Abstract Credit risk analysis is not only an important research topic in finance, but also of interest in everyday
More informationDiscussion Paper No. DP 07/05
SCHOOL OF ACCOUNTING, FINANCE AND MANAGEMENT Essex Finance Centre A Stochastic Variance Factor Model for Large Datasets and an Application to S&P data A. Cipollini University of Essex G. Kapetanios Queen
More informationStochastic Volatility (SV) Models
1 Motivations Stochastic Volatility (SV) Models Jun Yu Some stylised facts about financial asset return distributions: 1. Distribution is leptokurtic 2. Volatility clustering 3. Volatility responds to
More informationBayesian Estimation of the Markov-Switching GARCH(1,1) Model with Student-t Innovations
Bayesian Estimation of the Markov-Switching GARCH(1,1) Model with Student-t Innovations Department of Quantitative Economics, Switzerland david.ardia@unifr.ch R/Rmetrics User and Developer Workshop, Meielisalp,
More informationList of tables List of boxes List of screenshots Preface to the third edition Acknowledgements
Table of List of figures List of tables List of boxes List of screenshots Preface to the third edition Acknowledgements page xii xv xvii xix xxi xxv 1 Introduction 1 1.1 What is econometrics? 2 1.2 Is
More informationApplication of Bayesian Network to stock price prediction
ORIGINAL RESEARCH Application of Bayesian Network to stock price prediction Eisuke Kita, Yi Zuo, Masaaki Harada, Takao Mizuno Graduate School of Information Science, Nagoya University, Japan Correspondence:
More informationEstimation of Long Memory in Volatility
1 Estimation of Long Memory in Volatility Rohit S. Deo and C. M. Hurvich New York University Abstract We discuss some of the issues pertaining to modelling and estimating long memory in volatility. The
More informationApplication of MCMC Algorithm in Interest Rate Modeling
Application of MCMC Algorithm in Interest Rate Modeling Xiaoxia Feng and Dejun Xie Abstract Interest rate modeling is a challenging but important problem in financial econometrics. This work is concerned
More informationIntroductory Econometrics for Finance
Introductory Econometrics for Finance SECOND EDITION Chris Brooks The ICMA Centre, University of Reading CAMBRIDGE UNIVERSITY PRESS List of figures List of tables List of boxes List of screenshots Preface
More informationThe Great Moderation Flattens Fat Tails: Disappearing Leptokurtosis
The Great Moderation Flattens Fat Tails: Disappearing Leptokurtosis WenShwo Fang Department of Economics Feng Chia University 100 WenHwa Road, Taichung, TAIWAN Stephen M. Miller* College of Business University
More informationIdentifying Long-Run Risks: A Bayesian Mixed-Frequency Approach
Identifying : A Bayesian Mixed-Frequency Approach Frank Schorfheide University of Pennsylvania CEPR and NBER Dongho Song University of Pennsylvania Amir Yaron University of Pennsylvania NBER February 12,
More informationFin285a:Computer Simulations and Risk Assessment Section 7.1 Modeling Volatility: basic models Daníelson, ,
Fin285a:Computer Simulations and Risk Assessment Section 7.1 Modeling Volatility: basic models Daníelson, 2.1-2.3, 2.7-2.8 Overview Moving average model Exponentially weighted moving average (EWMA) GARCH
More informationChapter 5 Univariate time-series analysis. () Chapter 5 Univariate time-series analysis 1 / 29
Chapter 5 Univariate time-series analysis () Chapter 5 Univariate time-series analysis 1 / 29 Time-Series Time-series is a sequence fx 1, x 2,..., x T g or fx t g, t = 1,..., T, where t is an index denoting
More informationFiltering Stochastic Volatility Models with Intractable Likelihoods
Filtering Stochastic Volatility Models with Intractable Likelihoods Katherine B. Ensor Professor of Statistics and Director Center for Computational Finance and Economic Systems Rice University ensor@rice.edu
More informationEstimating Bivariate GARCH-Jump Model Based on High Frequency Data : the case of revaluation of Chinese Yuan in July 2005
Estimating Bivariate GARCH-Jump Model Based on High Frequency Data : the case of revaluation of Chinese Yuan in July 2005 Xinhong Lu, Koichi Maekawa, Ken-ichi Kawai July 2006 Abstract This paper attempts
More informationA STUDY ON ROBUST ESTIMATORS FOR GENERALIZED AUTOREGRESSIVE CONDITIONAL HETEROSCEDASTIC MODELS
A STUDY ON ROBUST ESTIMATORS FOR GENERALIZED AUTOREGRESSIVE CONDITIONAL HETEROSCEDASTIC MODELS Nazish Noor and Farhat Iqbal * Department of Statistics, University of Balochistan, Quetta. Abstract Financial
More informationA potentially useful approach to model nonlinearities in time series is to assume different behavior (structural break) in different subsamples
1.3 Regime switching models A potentially useful approach to model nonlinearities in time series is to assume different behavior (structural break) in different subsamples (or regimes). If the dates, the
More informationQuarterly Storage Model of U.S. Cotton Market: Estimation of the Basis under Rational Expectations. Oleksiy Tokovenko 1 Lewell F.
Quarterly Storage Model of U.S. Cotton Market: Estimation of the Basis under Rational Expectations Oleksiy Tokovenko 1 Lewell F. Gunter Selected Paper prepared for presentation at the American Agricultural
More informationStudy on Dynamic Risk Measurement Based on ARMA-GJR-AL Model
Applied and Computational Mathematics 5; 4(3): 6- Published online April 3, 5 (http://www.sciencepublishinggroup.com/j/acm) doi:.648/j.acm.543.3 ISSN: 38-565 (Print); ISSN: 38-563 (Online) Study on Dynamic
More informationForecasting Volatility of USD/MUR Exchange Rate using a GARCH (1,1) model with GED and Student s-t errors
UNIVERSITY OF MAURITIUS RESEARCH JOURNAL Volume 17 2011 University of Mauritius, Réduit, Mauritius Research Week 2009/2010 Forecasting Volatility of USD/MUR Exchange Rate using a GARCH (1,1) model with
More informationPosterior Inference. , where should we start? Consider the following computational procedure: 1. draw samples. 2. convert. 3. compute properties
Posterior Inference Example. Consider a binomial model where we have a posterior distribution for the probability term, θ. Suppose we want to make inferences about the log-odds γ = log ( θ 1 θ), where
More informationEstimating a Dynamic Oligopolistic Game with Serially Correlated Unobserved Production Costs. SS223B-Empirical IO
Estimating a Dynamic Oligopolistic Game with Serially Correlated Unobserved Production Costs SS223B-Empirical IO Motivation There have been substantial recent developments in the empirical literature on
More informationModels with Time-varying Mean and Variance: A Robust Analysis of U.S. Industrial Production
Models with Time-varying Mean and Variance: A Robust Analysis of U.S. Industrial Production Charles S. Bos and Siem Jan Koopman Department of Econometrics, VU University Amsterdam, & Tinbergen Institute,
More informationWeek 7 Quantitative Analysis of Financial Markets Simulation Methods
Week 7 Quantitative Analysis of Financial Markets Simulation Methods Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 November
More informationInternational Journal of Computer Engineering and Applications, Volume XII, Issue II, Feb. 18, ISSN
Volume XII, Issue II, Feb. 18, www.ijcea.com ISSN 31-3469 AN INVESTIGATION OF FINANCIAL TIME SERIES PREDICTION USING BACK PROPAGATION NEURAL NETWORKS K. Jayanthi, Dr. K. Suresh 1 Department of Computer
More informationAmath 546/Econ 589 Univariate GARCH Models: Advanced Topics
Amath 546/Econ 589 Univariate GARCH Models: Advanced Topics Eric Zivot April 29, 2013 Lecture Outline The Leverage Effect Asymmetric GARCH Models Forecasts from Asymmetric GARCH Models GARCH Models with
More informationAn Implementation of Markov Regime Switching GARCH Models in Matlab
An Implementation of Markov Regime Switching GARCH Models in Matlab Thomas Chuffart Aix-Marseille University (Aix-Marseille School of Economics), CNRS & EHESS Abstract MSGtool is a MATLAB toolbox which
More informationRunning head: IMPROVING REVENUE VOLATILITY ESTIMATES 1. Improving Revenue Volatility Estimates Using Time-Series Decomposition Methods
Running head: IMPROVING REVENUE VOLATILITY ESTIMATES 1 Improving Revenue Volatility Estimates Using Time-Series Decomposition Methods Kenneth A. Kriz Wichita State University Author Note The author wishes
More informationAnalyzing Oil Futures with a Dynamic Nelson-Siegel Model
Analyzing Oil Futures with a Dynamic Nelson-Siegel Model NIELS STRANGE HANSEN & ASGER LUNDE DEPARTMENT OF ECONOMICS AND BUSINESS, BUSINESS AND SOCIAL SCIENCES, AARHUS UNIVERSITY AND CENTER FOR RESEARCH
More informationSELECTION OF VARIABLES INFLUENCING IRAQI BANKS DEPOSITS BY USING NEW BAYESIAN LASSO QUANTILE REGRESSION
Vol. 6, No. 1, Summer 2017 2012 Published by JSES. SELECTION OF VARIABLES INFLUENCING IRAQI BANKS DEPOSITS BY USING NEW BAYESIAN Fadel Hamid Hadi ALHUSSEINI a Abstract The main focus of the paper is modelling
More informationPortfolio construction by volatility forecasts: Does the covariance structure matter?
Portfolio construction by volatility forecasts: Does the covariance structure matter? Momtchil Pojarliev and Wolfgang Polasek INVESCO Asset Management, Bleichstrasse 60-62, D-60313 Frankfurt email: momtchil
More informationNonlinear Filtering of Asymmetric Stochastic Volatility Models and VaR Estimation
Nonlinear Filtering of Asymmetric Stochastic Volatility Models and VaR Estimation Nikolay Nikolaev Goldsmiths College, University of London, UK n.nikolaev@gold.ac.uk Lilian M. de Menezes Cass Business
More informationA Macro-Finance Model of the Term Structure: the Case for a Quadratic Yield Model
Title page Outline A Macro-Finance Model of the Term Structure: the Case for a 21, June Czech National Bank Structure of the presentation Title page Outline Structure of the presentation: Model Formulation
More informationModeling the volatility of FTSE All Share Index Returns
MPRA Munich Personal RePEc Archive Modeling the volatility of FTSE All Share Index Returns Bayraci, Selcuk University of Exeter, Yeditepe University 27. April 2007 Online at http://mpra.ub.uni-muenchen.de/28095/
More informationThe Time-Varying Effects of Monetary Aggregates on Inflation and Unemployment
経営情報学論集第 23 号 2017.3 The Time-Varying Effects of Monetary Aggregates on Inflation and Unemployment An Application of the Bayesian Vector Autoregression with Time-Varying Parameters and Stochastic Volatility
More informationVolatility Clustering of Fine Wine Prices assuming Different Distributions
Volatility Clustering of Fine Wine Prices assuming Different Distributions Cynthia Royal Tori, PhD Valdosta State University Langdale College of Business 1500 N. Patterson Street, Valdosta, GA USA 31698
More informationFINANCIAL ECONOMETRICS PROF. MASSIMO GUIDOLIN
Massimo Guidolin Massimo.Guidolin@unibocconi.it Dept. of Finance FINANCIAL ECONOMETRICS PROF. MASSIMO GUIDOLIN SECOND PART, LECTURE 1: VOLATILITY MODELS ARCH AND GARCH OVERVIEW 1) Stepwise Distribution
More informationFinancial Econometrics Notes. Kevin Sheppard University of Oxford
Financial Econometrics Notes Kevin Sheppard University of Oxford Monday 15 th January, 2018 2 This version: 22:52, Monday 15 th January, 2018 2018 Kevin Sheppard ii Contents 1 Probability, Random Variables
More informationModelling Stock Market Return Volatility: Evidence from India
Modelling Stock Market Return Volatility: Evidence from India Saurabh Singh Assistant Professor, Graduate School of Business,Devi Ahilya Vishwavidyalaya, Indore 452001 (M.P.) India Dr. L.K Tripathi Dean,
More informationCross-Sectional Distribution of GARCH Coefficients across S&P 500 Constituents : Time-Variation over the Period
Cahier de recherche/working Paper 13-13 Cross-Sectional Distribution of GARCH Coefficients across S&P 500 Constituents : Time-Variation over the Period 2000-2012 David Ardia Lennart F. Hoogerheide Mai/May
More informationBayesian analysis of GARCH and stochastic volatility: modeling leverage, jumps and heavy-tails for financial time series
Bayesian analysis of GARCH and stochastic volatility: modeling leverage, jumps and heavy-tails for financial time series Jouchi Nakajima Department of Statistical Science, Duke University, Durham 2775,
More informationIntroduction Dickey-Fuller Test Option Pricing Bootstrapping. Simulation Methods. Chapter 13 of Chris Brook s Book.
Simulation Methods Chapter 13 of Chris Brook s Book Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 April 26, 2017 Christopher
More informationInternational Journal of Computer Engineering and Applications, Volume XII, Issue II, Feb. 18, ISSN
International Journal of Computer Engineering and Applications, Volume XII, Issue II, Feb. 18, www.ijcea.com ISSN 31-3469 AN INVESTIGATION OF FINANCIAL TIME SERIES PREDICTION USING BACK PROPAGATION NEURAL
More informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2009, Mr. Ruey S. Tsay. Solutions to Final Exam
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2009, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (42 pts) Answer briefly the following questions. 1. Questions
More informationEstimation of Stochastic Volatility Models : An Approximation to the Nonlinear State Space Representation
Estimation of Stochastic Volatility Models : An Approximation to the Nonlinear State Space Representation Junji Shimada and Yoshihiko Tsukuda March, 2004 Keywords : Stochastic volatility, Nonlinear state
More informationRegime-dependent Characteristics of KOSPI Return
Communications for Statistical Applications and Methods 014, Vol. 1, No. 6, 501 51 DOI: http://dx.doi.org/10.5351/csam.014.1.6.501 Print ISSN 87-7843 / Online ISSN 383-4757 Regime-dependent Characteristics
More informationEC316a: Advanced Scientific Computation, Fall Discrete time, continuous state dynamic models: solution methods
EC316a: Advanced Scientific Computation, Fall 2003 Notes Section 4 Discrete time, continuous state dynamic models: solution methods We consider now solution methods for discrete time models in which decisions
More informationINFORMATION EFFICIENCY HYPOTHESIS THE FINANCIAL VOLATILITY IN THE CZECH REPUBLIC CASE
INFORMATION EFFICIENCY HYPOTHESIS THE FINANCIAL VOLATILITY IN THE CZECH REPUBLIC CASE Abstract Petr Makovský If there is any market which is said to be effective, this is the the FOREX market. Here we
More informationImplied Volatility v/s Realized Volatility: A Forecasting Dimension
4 Implied Volatility v/s Realized Volatility: A Forecasting Dimension 4.1 Introduction Modelling and predicting financial market volatility has played an important role for market participants as it enables
More information. Large-dimensional and multi-scale effects in stocks volatility m
Large-dimensional and multi-scale effects in stocks volatility modeling Swissquote bank, Quant Asset Management work done at: Chaire de finance quantitative, École Centrale Paris Capital Fund Management,
More informationEstimation of Monthly Volatility: An Empirical Comparison of Realized Volatility, GARCH and ACD-ICV Methods
Estimation of Monthly Volatility: An Empirical Comparison of Realized Volatility, GARCH and ACD-ICV Methods Shouwei Liu School of Economics, Singapore Management University Yiu-Kuen Tse School of Economics,
More informationThe Complexity of GARCH Option Pricing Models
JOURNAL OF INFORMATION SCIENCE AND ENGINEERING 8, 689-704 (01) The Complexity of GARCH Option Pricing Models YING-CHIE CHEN +, YUH-DAUH LYUU AND KUO-WEI WEN + Department of Finance Department of Computer
More informationUsing Agent Belief to Model Stock Returns
Using Agent Belief to Model Stock Returns America Holloway Department of Computer Science University of California, Irvine, Irvine, CA ahollowa@ics.uci.edu Introduction It is clear that movements in stock
More informationTop-down particle filtering for Bayesian decision trees
Top-down particle filtering for Bayesian decision trees Balaji Lakshminarayanan 1, Daniel M. Roy 2 and Yee Whye Teh 3 1. Gatsby Unit, UCL, 2. University of Cambridge and 3. University of Oxford Outline
More informationEstimation of dynamic term structure models
Estimation of dynamic term structure models Greg Duffee Haas School of Business, UC-Berkeley Joint with Richard Stanton, Haas School Presentation at IMA Workshop, May 2004 (full paper at http://faculty.haas.berkeley.edu/duffee)
More informationOption Pricing Using Bayesian Neural Networks
Option Pricing Using Bayesian Neural Networks Michael Maio Pires, Tshilidzi Marwala School of Electrical and Information Engineering, University of the Witwatersrand, 2050, South Africa m.pires@ee.wits.ac.za,
More informationInstitute of Actuaries of India Subject CT6 Statistical Methods
Institute of Actuaries of India Subject CT6 Statistical Methods For 2014 Examinations Aim The aim of the Statistical Methods subject is to provide a further grounding in mathematical and statistical techniques
More informationForecasting Volatility of Wind Power Production
Forecasting Volatility of Wind Power Production Zhiwei Shen and Matthias Ritter Department of Agricultural Economics Humboldt-Universität zu Berlin July 18, 2015 Zhiwei Shen Forecasting Volatility of Wind
More informationFORECASTING PERFORMANCE OF MARKOV-SWITCHING GARCH MODELS: A LARGE-SCALE EMPIRICAL STUDY
FORECASTING PERFORMANCE OF MARKOV-SWITCHING GARCH MODELS: A LARGE-SCALE EMPIRICAL STUDY Latest version available on SSRN https://ssrn.com/abstract=2918413 Keven Bluteau Kris Boudt Leopoldo Catania R/Finance
More informationForecasting jumps in conditional volatility The GARCH-IE model
Forecasting jumps in conditional volatility The GARCH-IE model Philip Hans Franses and Marco van der Leij Econometric Institute Erasmus University Rotterdam e-mail: franses@few.eur.nl 1 Outline of presentation
More informationModel Estimation. Liuren Wu. Fall, Zicklin School of Business, Baruch College. Liuren Wu Model Estimation Option Pricing, Fall, / 16
Model Estimation Liuren Wu Zicklin School of Business, Baruch College Fall, 2007 Liuren Wu Model Estimation Option Pricing, Fall, 2007 1 / 16 Outline 1 Statistical dynamics 2 Risk-neutral dynamics 3 Joint
More informationAnalysis of the Bitcoin Exchange Using Particle MCMC Methods
Analysis of the Bitcoin Exchange Using Particle MCMC Methods by Michael Johnson M.Sc., University of British Columbia, 2013 B.Sc., University of Winnipeg, 2011 Project Submitted in Partial Fulfillment
More informationForecasting Volatility in the Chinese Stock Market under Model Uncertainty 1
Forecasting Volatility in the Chinese Stock Market under Model Uncertainty 1 Yong Li 1, Wei-Ping Huang, Jie Zhang 3 (1,. Sun Yat-Sen University Business, Sun Yat-Sen University, Guangzhou, 51075,China)
More informationLecture Note of Bus 41202, Spring 2008: More Volatility Models. Mr. Ruey Tsay
Lecture Note of Bus 41202, Spring 2008: More Volatility Models. Mr. Ruey Tsay The EGARCH model Asymmetry in responses to + & returns: g(ɛ t ) = θɛ t + γ[ ɛ t E( ɛ t )], with E[g(ɛ t )] = 0. To see asymmetry
More informationRelevant parameter changes in structural break models
Relevant parameter changes in structural break models A. Dufays J. Rombouts Forecasting from Complexity April 27 th, 2018 1 Outline Sparse Change-Point models 1. Motivation 2. Model specification Shrinkage
More informationBooth School of Business, University of Chicago Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay. Solutions to Midterm
Booth School of Business, University of Chicago Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay Solutions to Midterm Problem A: (34 pts) Answer briefly the following questions. Each question has
More information